Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,596,878:1 +G 5.9% intergenic (‑141/+27) ypfJ ← / ← purC putative neutral zinc metallopeptidase/phosphoribosylaminoimidazole‑succinocarboxamide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,596,8781.G 5.9% 350.3 / 15.5 85intergenic (‑141/+27)ypfJ/purCputative neutral zinc metallopeptidase/phosphoribosylaminoimidazole‑succinocarboxamide synthetase
Reads supporting (aligned to +/- strand):  ref base . (40/40);  new base G (3/2);  total (43/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.91e-01

CAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACA‑CCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTCCCACAGG  >  NC_000913/2596730‑2597022
                                                                                                                                                     |                                                                                                                                                 
cAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCt                                                                                                                                                    >  1:1128960/1‑149 (MQ=255)
cAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCt                                                                                                                                                    >  2:650075/1‑149 (MQ=255)
 aaCGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑g                                                                                                                                                  <  2:1013239/149‑1 (MQ=255)
    gctttgATACCTCAACAATATATTAATTATAGAAATGCTCGAAGGTGGTTCATAGGAAGATTACCACAGGAAAGCTGGATGAAATGGCGGAGGAATGAGGGTTCAATGCTTATATTCCGAGACCGGTTGTCTTGCCCTGCGAGCT‑GCAc                                                                                                                                               <  2:516122/144‑1 (MQ=255)
      aTAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGGGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACgg                                                                                                                                             <  2:1411985/149‑1 (MQ=255)
           tACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAga                                                                                                                                        <  2:589777/149‑1 (MQ=255)
             ccTCAACAATATATTCATTATAGAAATGATCGTTGGTGGTTCAGAGGAAGATTACCAAAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTTCGAGGT‑GCGCTGCAAgatt                                                                                                                                      >  2:785286/1‑148 (MQ=255)
             ccTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAgatt                                                                                                                                      >  2:424804/1‑148 (MQ=255)
              cTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAgatta                                                                                                                                     >  1:61994/1‑147 (MQ=255)
              cTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAgatga                                                                                                                                     <  2:59014/149‑1 (MQ=255)
              cTCAACAATATATTCATAATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAgatga                                                                                                                                     <  2:119567/149‑1 (MQ=255)
               tCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTATCTTGCCGTG‑CAGCTGGCACGGCAAgacaac                                                                                                                                    <  2:1484600/149‑5 (MQ=255)
                 aacaaTATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGAtgatga                                                                                                                                  >  2:888337/1‑149 (MQ=255)
                 aacaaTATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGAtgatga                                                                                                                                  >  2:1347505/1‑149 (MQ=255)
                  acaaTATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGaa                                                                                                                                 <  1:180753/149‑1 (MQ=255)
                   caaTATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGaaa                                                                                                                                >  1:596179/1‑148 (MQ=255)
                   caaTATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAAc                                                                                                                                >  1:1340876/1‑149 (MQ=255)
                            cATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAAta                                                                                                                       >  2:338413/1‑148 (MQ=255)
                            cATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATc                                                                                                                       <  2:924720/149‑1 (MQ=255)
                             aTTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCa                                                                                                                      >  2:1531765/1‑149 (MQ=255)
                              ttATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAg                                                                                                                     <  1:762820/149‑1 (MQ=255)
                              ttATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAg                                                                                                                     <  2:1492034/149‑1 (MQ=255)
                              ttATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAg                                                                                                                     >  2:1333004/1‑149 (MQ=255)
                               tataGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCATTGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCATGATGATGAACAGAAAAATCAGt                                                                                                                    <  1:911024/149‑1 (MQ=255)
                               tataGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGACCAGAAAAATCAGt                                                                                                                    >  2:1042798/1‑149 (MQ=255)
                                ataGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTc                                                                                                                   <  2:1538246/149‑1 (MQ=255)
                                ataGAAATGATCGTAGGTGGTTCAGAGGAAGAATACCACAGGAAAGCGGGATGAAATGGAGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTTCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTc                                                                                                                   >  1:1133071/1‑149 (MQ=255)
                                   gAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATTCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAACATGATGAACAGAAAAATCAGTCCAg                                                                                                                >  2:197904/1‑149 (MQ=255)
                                    aaaTGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGc                                                                                                               >  2:285751/1‑149 (MQ=255)
                                       tGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGACAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGt                                                                                                            <  1:251606/149‑1 (MQ=255)
                                           cGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTaca‑c                                                                                                       >  1:478855/1‑149 (MQ=255)
                                            gTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACA‑cc                                                                                                      <  1:1212524/149‑1 (MQ=255)
                                             tAGGTGGTTCAGAGGAAGATTACCAAAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACA‑ccc                                                                                                     >  2:1006183/1‑149 (MQ=255)
                                                    ttCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACA‑CCCAggcggc                                                                                              >  1:935936/1‑149 (MQ=255)
                                                         agGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACA‑CCCAGGCggcgggcg                                                                                         >  1:533269/1‑149 (MQ=255)
                                                         agGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACA‑CCCAGGCggcgggcg                                                                                         >  2:613426/1‑149 (MQ=255)
                                                          ggAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACA‑CCCAGGCGGCGGGCGa                                                                                        <  1:557969/149‑1 (MQ=255)
                                                          ggAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGACGTGGCAGAA‑TCACGGCAAAATGGAAAACAAAAAAATAAATCCACCGAGACACCCCAGGcgccggcga                                                                                         >  2:1212197/1‑146 (MQ=255)
                                                           gAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACA‑CCCAGGCGGCGGGCGAc                                                                                       >  1:683309/1‑149 (MQ=255)
                                                           gAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACA‑CCCAGGCGGCGGGCGAc                                                                                       <  2:1340876/149‑1 (MQ=255)
                                                                   ccACAGGCAAGCGGGAGGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGACTTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACA‑CCCAGGCGGCGGGCGACGGCTTCAt                                                                               <  2:95647/149‑1 (MQ=255)
                                                                    cacaGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACA‑CCCAGGCGGCGGGCGACGGCTTCATa                                                                              >  2:919545/1‑149 (MQ=255)
                                                                      caGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACA‑CCCAGGCGGCGGGCGACGGCTTCATAgg                                                                            >  1:1282911/1‑149 (MQ=255)
                                                                       aGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACA‑CCCAGGCGGCGGGCGACGGCTTCATAGGc                                                                           <  1:1449717/149‑1 (MQ=255)
                                                                       aGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACA‑CCCAGGCGGCGGGCGACGGCTTCATAGGc                                                                           <  2:1128960/149‑1 (MQ=255)
                                                                         gAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACA‑CCCAGGCGGCGGGCGACGGCTTCATAGGCtt                                                                         >  1:1099120/1‑149 (MQ=255)
                                                                         gAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACA‑CCCAGGCGGCGGGCGACGGCTTCATAGGCtt                                                                         <  2:1474582/149‑1 (MQ=255)
                                                                             gagggCTGAAATGGCGGAGGACTGAGCGTTCAATGCATATAATCCGCGAGCGGGTGGCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACC‑CCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGAt                                                                     <  2:65683/147‑1 (MQ=255)
                                                                             gCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACC‑CCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGAt                                                                     >  1:649481/1‑149 (MQ=255)
                                                                               gggATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTTCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTACAGCTTTAAA‑CCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCa                                                                   >  2:504311/1‑149 (MQ=255)
                                                                               gggATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACA‑CCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCa                                                                   >  2:1049091/1‑149 (MQ=255)
                                                                                  aTGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACA‑CCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGc                                                                <  2:935936/149‑1 (MQ=255)
                                                                                      aaTGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACA‑CCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCAcc                                                            <  2:1282911/149‑1 (MQ=255)
                                                                                                 aaTGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACA‑CCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGc                                                 <  2:61994/149‑1 (MQ=255)
                                                                                                  aTGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTTCCGTGCCAGCT‑GCCCGGCAAGATGATGAACAGAAAAATCAGTACAGCTGTACA‑CCCCCGCGGCGGGCGCCGGGTTCATAGGCTTCGATCAGGCCACCCACGCTCTGGCg                                                >  2:1319193/1‑149 (MQ=255)
                                                                                                  aTGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACA‑CCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCg                                                >  2:300393/1‑149 (MQ=255)
                                                                                                  aTGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACA‑CCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCg                                                >  1:1595093/1‑149 (MQ=255)
                                                                                                        gTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACA‑CCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACg                                          <  2:1043359/149‑1 (MQ=255)
                                                                                                        gTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACA‑CCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACg                                          >  2:1042873/1‑149 (MQ=255)
                                                                                                           aaaTGCATATATTGCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GGACGGAAAGATGATGAACAGAAAAATCAGTCCATCTGTACA‑CCCAGGGGGCGGGCGACGGCTTCATAGTCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTc                                       <  2:1508161/148‑1 (MQ=255)
                                                                                                               gCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACA‑CCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTg                                   <  1:1515940/149‑1 (MQ=255)
                                                                                                                cATATATTCCGAGAGCGGTTGTCTTGCCTTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACA‑CCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGt                                  <  1:650075/149‑1 (MQ=255)
                                                                                                                cATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACTGCAAGATGATGAACAGAAAATTCAGTCCAGCTGTACA‑CCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGTCTCTGGCGGAAACGGTCTTTGt                                  <  2:608255/149‑1 (MQ=255)
                                                                                                                   atatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACA‑CCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCa                               >  1:924286/1‑149 (MQ=255)
                                                                                                                   atatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACA‑CCCAGGCGGCGGGCGACGGCGTCATATGCTTCGATCAGGCCCCCGAGGCTCTGGCGGAAACGGTCTTTGTCCa                               >  2:480417/1‑149 (MQ=255)
                                                                                                                    aatTCCGAGAGCGGTGGTCTTGCCGTGCGAGCT‑GCACGGCAAGATGATGCACAGAAAAATCCGTCATGCTGTACT‑ACCAGGCGGCTGGCGACGGCTTCATAGGCTTCGATCAGTCCACCGAGGCGCTGGCGGAAACTGTCTTTGTCCAt                              <  1:237657/148‑1 (MQ=255)
                                                                                                                     atTCCGAGAGCGGTTGTCTTGCCATGCCAGCA‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAACTGTACA‑CCCAGGCGGCGGGCGACGGCTTCATAGGCTACGATCAAGCCACCGAGGCTCTGGCGGAAACCGTCTTAGTCCaat                             >  2:448949/1‑147 (MQ=255)
                                                                                                                      ttCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACA‑CCCAGGCGGCGGGGGACGGCGTCATAGGCTTGGGTCAGGCCACCGAGGCTCTGGGGGAAACGGTCTTTGTCCAttt                            >  1:609694/1‑149 (MQ=255)
                                                                                                                        ccGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACA‑CCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTtc                          >  1:198274/1‑149 (MQ=255)
                                                                                                                        ccGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACA‑CCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTtc                          <  2:478855/149‑1 (MQ=255)
                                                                                                                        ccGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACA‑CCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGGCTTTGTCCATTTtc                          <  1:285751/149‑1 (MQ=255)
                                                                                                                         cGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCCCGGGAAGATGATGCACAGAAAAATCAGTCCAGCTGTTCA‑CCCAGGCGGCGGGCGCCGGGGTCATAGGGGTCGAGCAGGCCCCCGCGGCGCTGGCGGAAACGGTCTTTGTGCATTTtct                         >  1:647373/1‑149 (MQ=255)
                                                                                                                         cGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATAAGTCCAGCTGTACA‑CCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGACGCTCTGGCGGAAACGGTCTTTGTACAttttat                         >  2:1376112/1‑147 (MQ=255)
                                                                                                                             cgCGGTTGTCTTGCCGTGCCAGCG‑GCACGGCAAGATGATGCACAGAAAAATCAGTCCCGCTGTACA‑CCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAg                     <  2:1161405/148‑1 (MQ=255)
                                                                                                                             agCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACA‑CCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAg                     <  2:1167752/149‑1 (MQ=255)
                                                                                                                              gCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACA‑CCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGc                    <  1:1347505/149‑1 (MQ=255)
                                                                                                                                ggTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACA‑CCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGt                  <  1:47273/149‑1 (MQ=255)
                                                                                                                                 gTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACA‑CCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGtt                 <  1:1319193/149‑1 (MQ=255)
                                                                                                                                  tcatCTTGCCGTG‑CAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACA‑CCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGttt                >  2:322493/4‑149 (MQ=255)
                                                                                                                                    gTCTTGCCGGTCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTAGA‑CCCATGCGGCGGGCGACGGCTTCTTAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCt              <  1:1212197/149‑1 (MQ=255)
                                                                                                                                        tGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACA‑CCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGt          >  2:1141801/1‑149 (MQ=255)
                                                                                                                                        tGCCGGGCCAGCT‑GCAAGGCAAGATGATGAACAGAAAAATCAGTCCCGCTGTCCA‑CCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGt          <  1:1006183/149‑1 (MQ=255)
                                                                                                                                         gCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACA‑CCCAGGCGGAGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTc         <  2:830284/149‑1 (MQ=255)
                                                                                                                                          ccGTG‑CAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACA‑CCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCCCCGCGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTcc        >  1:461559/1‑149 (MQ=255)
                                                                                                                                          ccGTG‑CAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACA‑CCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGaaa                                            >  1:1412899/1‑113 (MQ=255)
                                                                                                                                          ccGTG‑CAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACA‑CCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGaaa                                            <  2:1412899/113‑1 (MQ=255)
                                                                                                                                               ccAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACA‑CCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTCCCACAg   >  2:679726/1‑149 (MQ=255)
                                                                                                                                                cAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACA‑CCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTCCCACAgg  <  1:1438202/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
CAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACA‑CCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTCCCACAGG  >  NC_000913/2596730‑2597022

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑