Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,596,873 Δ1 bp 5.7% intergenic (‑136/+32) 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,8730C. 5.7% 256.6 / 11.0 87intergenic (‑136/+32)ypfJ/purCputative neutral zinc metallopeptidase/phosphoribosylaminoimidazole‑succinocarboxamide synthetase
Reads supporting (aligned to +/- strand):  ref base C (40/40);  new base . (3/2);  total (44/43)
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.59e-01

CTTGCCAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACA‑CCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTCCCACAG  >  NC_000913/2596725‑2597021
                                                                                                                                                    |                                                                                                                                                      
cTTGCCAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGc                                                                                                                                                        <  1:1267607/149‑1 (MQ=255)
  tGCCAACGCATACATACCTCAACAATATAATCCTTATAGAAATGATCCTAGGTGGTGCAGAGGAAGATTACCACAGGAACGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCa                                                                                                                                                      <  2:1222774/149‑1 (MQ=255)
     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:1347505/1‑149 (MQ=255)
                      aacaaTATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGAtgatga                                                                                                                                 >  2:888337/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                                                                                                                    <  2:1492034/149‑1 (MQ=255)
                                   ttATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAg                                                                                                                    <  1:762820/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                                                                          <  2:1128960/149‑1 (MQ=255)
                                                                            aGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACA‑CCCAGGCGGCGGGCGACGGCTTCATAGGc                                                                          <  1:1449717/149‑1 (MQ=255)
                                                                              gAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACA‑CCCAGGCGGCGGGCGACGGCTTCATAGGCtt                                                                        <  2:1474582/149‑1 (MQ=255)
                                                                              gAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACA‑CCCAGGCGGCGGGCGACGGCTTCATAGGCtt                                                                        >  1:1099120/1‑149 (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                                               >  1:1595093/1‑149 (MQ=255)
                                                                                                       aTGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACA‑CCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCg                                               >  2:300393/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                         <  2:478855/149‑1 (MQ=255)
                                                                                                                             ccGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACA‑CCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTtc                         >  1:198274/1‑149 (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)
                                                                                                                                                    |                                                                                                                                                      
CTTGCCAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACA‑CCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTCCCACAG  >  NC_000913/2596725‑2597021

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: ‑