Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 505,155 T→C 5.8% F179F (TTT→TTC PP_0416 → phosphoglycolate phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947505,1550TC 5.8% 202.0 / 2.2 72F179F (TTT→TTCPP_0416phosphoglycolate phosphatase
Reads supporting (aligned to +/- strand):  ref base T (34/32);  new base C (2/2);  total (37/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01

GGCCCGCTGCTGGACCAGATGAAGATCGGTAATTTCTTCCGCTGGATCATTGGTGGCGACACCTTGCCGCAGAAAAAGCCCGACCCGGCGGCGTTGCTGTTCGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTTGTCGGTGATTCGCGCAGTGACGTGCAGGCGGCCAAGGCGGCTGGCGTGCAATGCGTGGGCCTGACCTACGGCTATAACCATGGCCGGCCGATCGACGCCGAATCCCCCAGCCTGGTAATCGACGACCTGCGCGCATTGCTGCC  >  NC_002947/505012‑505298
                                                                                                                                               |                                                                                                                                               
ggCCCGCTGCTGGACCAGATGAAGATCGGTAATTTCTTCCGCTGGATCATTGGTGGCGACACCTTGCCGCAGAAAAAGCCCGACCCGGCGGCGTTGCTGTTCGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTTGTCgg                                                                                                                                            >  6:125593/1‑149 (MQ=255)
ggCCCGCTGCTGGACCAGATGAAGATCGGTAATTTCTTCCGCTGGATCATTGGTGGCGACACCTTGCCGCAGAAAAAGCCCGACCCGGCGGCGTTGCTGTTCGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTTGTCg                                                                                                                                             <  6:269594/148‑1 (MQ=255)
 gCCCGCTGCTGGACCAGATGAAGATCGGTAATTTCTTCCGCTGGATCATTGGTGGCGACACCTTGCCGCAGAAAAAGCCCGACCCGGCGGCGTTGCTGTTCGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTTGTCGGt                                                                                                                                           <  1:33479/149‑1 (MQ=255)
    cGCTGCTGGACCAGATGAAGATCGGTAATTTCTTCCGCTGGATCATTGGTGGCGACACCTTGCCGCAGAAAAAGCCCGACCCGGCGGCGTTGCTGTTCGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTTGTCGGTGAt                                                                                                                                        >  3:248489/1‑149 (MQ=255)
        gctgGACCAGATGAAGATCGGTAATTTCATCCGCTGGATCATTGGTGGCGACACCTTGCCGCAGAAAAAGCCCGACCCGGCGGCGTTGCTGTTCGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTTGTCGGTGATTcgc                                                                                                                                    <  6:273295/149‑1 (MQ=255)
         ctgGACCAGATGAAGATCGGTAATTTCTTCCGCTGGATCATTGGTGGCGACACCTTGCCGCAGAAAAAGCCCGACCCGGCGGCGTTGCTGTTCGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTTGTCGGTGATTcgcg                                                                                                                                   <  5:67882/149‑1 (MQ=255)
            gACCAGATGAAGATCGGTAATTTCTTCCGCTGGATCATTGGTGGCGACACCTTGCCGCAGAAAAAGCCCGACCCGGCGGCGTTGCTGTTCGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTTGTCGGTGATTCgcgc                                                                                                                                  <  6:87616/147‑1 (MQ=255)
            gACCAGATGAAGATCGGTAATTTCTTCCGCTGGATCATTGGTGGCGACACCTTGCCGCAGAAAAAGCCCGACCCGGCGGCGTTGCTGTTCGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTTGTCGGTGATTAGCGCAg                                                                                                                                >  5:198250/1‑149 (MQ=255)
             aCCAGATGAAGATCGGTAATTTCTTCCGCTGGATCATTGGTGGCGACACCTTGCCGCAGAAAAAGCCCGACCCGGCGGCGTTGCTGTTGGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTTGTCGGTGATTCGCGCAGt                                                                                                                               >  2:158931/1‑149 (MQ=255)
                aGATGAAGATCGGTAATTTCTTCCGCTGGATCATTGGTGGCGACACCTTGCCGCAGAAAAAGCCCGACCCGGCGGCGTTGCTGTTCGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTTGTCGGTGATTCGCGCAGTGAc                                                                                                                            <  4:49678/149‑1 (MQ=255)
                aGATGAAGATCGGTAATTTCTTCCGCTGGATCATTGGTGGCGACACCTTGCCGCAGAAAAAGCCCGACCCGGCGGCGTTGCTGTTCGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTTGTCGGTGATTCGCGCAGTGAc                                                                                                                            <  2:324296/149‑1 (MQ=255)
                aGATGAAGATCGGTAATTTCTTCCGCTGGATCATTGGTGGCGACACCTTGCCGCAGAAAAAGCCCGACCCGGCGGCGTTGCTGTTCGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTTGTCGGTGATTCGCGCAGTGAc                                                                                                                            <  4:248489/149‑1 (MQ=255)
                 gATGAAGATCGGTAATTTCTTCCGCTGGATCATTGGTGGCGACACCTTGCCGCAGAAAAAGCCCGACCCGGCGGCGTTGCTGTTCGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTTGTCGGTGATTCGCGCAGTGACg                                                                                                                           >  2:303454/1‑149 (MQ=255)
                 gATGAAGATCGGTAATTTCTTCCGCTGGATCATTGGTGGCGACACCTTGCCGCAGAAAAAGCCCGACCCGGCGGCGTTGCTGTTCGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTCGTCGGTGATTCGCGCAGTGACg                                                                                                                           >  4:297580/1‑149 (MQ=255)
                 gATGAAGATCGGTAATTTCTTCCGCTGGATCATTGGTGGCGACACCTTGCCGCAGAAAAAGCCCGACCCGGCGGCGTTGCTGTTCGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTCGTCGGTGATTCGCGCAGTGACg                                                                                                                           >  4:292651/1‑149 (MQ=255)
                        aTCGGTAATTTCTTCCGCTGGATCATTGGTGGCGACACCTTGCCGCAGAAAAAGCCCGACCCGGCGGCGTTGCTGTTCGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTTGTCGGTGATTCGCGCAGTGACGTGCAggc                                                                                                                    <  6:293356/149‑1 (MQ=255)
                        aTCGGTAATTTCTTCCGCTGGATCATTGGTGGCGACACCTTGCCGCAGAAAAAGCCCGACCCGGCGGCGTTGCTGTTCGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTTGTCGGTGATTCGCGCAGTGACGTGCAggc                                                                                                                    >  4:298658/1‑149 (MQ=255)
                           ggTAATTTCTTCCGCTGGATCATTGGTGGCGACACCTTGCCGCAGAAAAAGCCCGACCCGGCGGCGTTGCTGTTCGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTTGTCGGTGATTCGCGCAGTGACGTGCAggcgg                                                                                                                  <  8:90186/148‑1 (MQ=255)
                            gTAATTTCTTCCGCTGGATCATTGGTGGCGACACCTTGCCGCAGAAAAAGCCCGACCCGGCGGCGTTGCTGTTCGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTTGTCGGTGATTCGCGCAGTGACGTGCAggcggc                                                                                                                 >  2:59597/1‑148 (MQ=255)
                            gTAATTTCTTCCGCTGGATCATTGGTGGCGACACCTTGCCGCAGAAAAAGCCCGACCCGGCGGCGTTGCTGTTCGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTTGTCGGTGATTCGCGCAGTGACGTGCAggcggc                                                                                                                 >  4:289882/1‑148 (MQ=255)
                             tAATTTCTTCCGCTGGATCATTGGTGGCGACACCTTGCCGCAGAAAAAGCCCGACCCGGCGGCGTTGCTGTTCGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTTGTCGGTGATTCGCGCAGTGACGTGCAGGCGGcc                                                                                                                >  3:28094/1‑148 (MQ=255)
                               aTTTCTTCCGCTGGATGATTGGTGGCGACACCTTGGAGCAGAAAAAGCCCGACCCGGCGGAGTTGCTGTTCGTCATGCAGATGGCCGGTGTCACCCCGCAGGAGTCCCTGTTTGTCGGTGATTCGCGCAGTGACGTg                                                                                                                         <  4:99538/137‑1 (MQ=255)
                               aTTTCTTCCGCTGGATCATTGGTGGCGACACCTTGCCGCAGAAAAAGCCCGACCCGGCGGCGTTGCTGTTCGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTTGTCGGTGATTCGCGCAGTGACGTGCAGGCGGCCAAg                                                                                                             >  8:191572/1‑149 (MQ=255)
                               aTTTCTTCCGCTGGATCATTGGTGGCGACACCTTGCCGCAGAAAAAGCCCGACCCGGCGGCGTTGCTGTTCGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTAGTCGGTGATTCGCGCAGTGACGTg                                                                                                                         >  3:99538/1‑137 (MQ=255)
                                tttCTTCCGCTGGATCATTGGTGGCGACACCTTGCCGCAGAAAAAGCCCGACCCGGCGGCGTTGCTGTTCGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTTGTCGGTGATTCGCGCAGTGACGTGCAGGCGGCCAAgg                                                                                                            <  4:137121/149‑1 (MQ=255)
                                 ttcttcCGCTGGATCATTGGTGGCGACACCTTGCCGCAGAAAAAGCCCGACCCGGCGGCGTTGCTGTTCGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTTGTCGGTGATTCGCGCAGTGACGTGCAGGCGGCCAAggc                                                                                                           >  2:259692/1‑149 (MQ=255)
                                   cttcCGCTGGATCATTGGTGGCGACACCTTGCCGCAGAAAAAGCCCGACCCGGCGGCGTTGCTGTTCGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTTGTCGGTGATTCGCGCAGTGACGTGCAGGCGGCCAAgg                                                                                                            <  8:47489/146‑1 (MQ=255)
                                   cttcCGCTGGATCATTGGTGGCGACACCTTGCCGCAGAAAAAGCCCGACCCGGCGGCGTTGCTGTTCGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTTGTCGGTGATTCGCGCAGTGACGTGCAGGCGGCCAAgg                                                                                                            >  7:47489/1‑146 (MQ=255)
                                       cGCTGGATCATTGGTGGCGACACCTTGCCGCAGAAAAAGCCCGACCCGGCGGCGTTGCTGTTCGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTTGTCGGTGATTCGCGCAGTGACGTGCAGGCGGCCAAGGCGGCTg                                                                                                      <  2:34899/148‑1 (MQ=255)
                                        gCTGGATCATTGGTGGCGACACCTTGCCGCAGAAAAAGCCCGACCCGGCGGCGTTGCGGTTCGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTTGTCGGTGATTCGCGCAGTGACGTGCAGGCGGCCAAGGCGGCTGGc                                                                                                    <  2:261352/149‑1 (MQ=255)
                                             aTCATTGGTGGCGACACCTTGCCGCAGAAAAAGCCCGACCCGGCGGCGTTGCTGTTCGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTTGTCGGTGATTcgc                                                                                                                                    >  3:66565/1‑112 (MQ=255)
                                             aTCATTGGTGGCGACACCTTGCCGCAGAAAAAGCCCGACCCGGCGGCGTTGCTGTTCGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTTGTCGGTGATTcgc                                                                                                                                    <  4:66565/112‑1 (MQ=255)
                                              tCATTGGTGGCGACACCTTGCCGCAGAAAAAGCCCGACCCGGCGGCGTTGCTGTTCGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTCGTCGGTGATTCGCGCAGTGACGTGCAGGCGGCCAAGGCGGCTGGCGTGCaa                                                                                              <  3:297580/149‑1 (MQ=255)
                                              tCATTGGTGGCGACACCTTGCCGCAGAAAAAGCCCGACCCGGCGGCGTTGCTGTTCGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTCGTCGGTGATTCGCGCAGTGACGTGCAGGCGGCCAAGGCGGCTGGCGTGCaa                                                                                              <  3:292651/149‑1 (MQ=255)
                                                aTTGGTGGCGACACCTTGCCGCAGAAAAAGCCCGACCCGGCGGCGTTGCTGTTCGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTTGTCGGTGATTCGCGCAGTGACGTGCAGGCGGCCAAGGCGGCTGGCGTGCAAt                                                                                             <  5:125593/148‑1 (MQ=255)
                                                aTTGGTGGCGACACCTTGCCGCAGAAAAAGCCCGACCCGGCGGCGTTGCTGTTCGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTTGTCGGTGATTCGCGCAGTGACGTGCAGGCGGCCAAGGCGGCTGGCGTGCAAt                                                                                             <  1:59597/148‑1 (MQ=255)
                                                 ttGGTGGCGACACCTTGCCGCAGAAAAAGCCCGACCCGGCGGCGTTGCTGTTCGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTTGTCGGTGATTCGCGCAGTGACGTGCAGGCGGCCAAGGCGGCTGGCGTGCAATGc                                                                                           >  5:79972/1‑149 (MQ=255)
                                                 ttGGTGGCGACACCTTGCCGCAGAAAAAGCCCGACCCGGCGGCGTTGCTGTTCGTCATGCAGATGGCCGGTGTCACCCCACAGCAGTCCCTGTTTGTCGGTGATTCGCGCAGTGACGTGCAGGCGGCCAAGGCGGCTGGCGTGCAATGc                                                                                           >  4:54066/1‑149 (MQ=255)
                                                     tggCGACACCTTGCCGCAGAAAAAGCCCGACCCGGCGGCGTTGCTGTTCGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTTGTCGGTGATTCGCGCAGTGACGTGCAg                                                                                                                      >  2:249723/1‑118 (MQ=255)
                                                     tggCGACACCTTGCCGCAGAAAAAGCCCGACCCGGCGGCGTTGCTGTTCGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTTGTCGGTGATTCGCGCAGTGACGTGCAg                                                                                                                      <  1:249723/118‑1 (MQ=255)
                                                     tggCGACACCTTGCCGCAGAAAAAGCCCGACCCGGCGGCGTTGCTGTTCGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTTGTCGGTGATTCGCGCAGTGACGTGCAGGCGGCCAAGGCGGCTGGCg                                                                                                   <  3:268926/137‑1 (MQ=255)
                                                     tggCGACACCTTGCCGCAGAAAAAGCCCGACCCGGCGGCGTTGCTGTTCGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTTGTCGGTGATTCGCGCAGTGACGTGCAGGCGGCCAAGGCGGCTGGCg                                                                                                   >  4:268926/1‑137 (MQ=255)
                                                       gCGACACCTTGCCGCAGAAAAAGCCCGACCCGGCGGCGTTGCTGTTCGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTTGTCGGTGATTCGCGCAGTGACGTGCAGGCGGCCAAGGCGGCTGGCGTGCAATGCGTGGGc                                                                                     <  7:191572/149‑1 (MQ=255)
                                                          acacCTTGCCGCAGAAAAAGCCCGACCCGGCGGCGTTGCTGTTCGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTTGTCGGTGATTCGCGCAGTGACGTGCAGGCGGCCAAGGCGGCTGGCGTGCAATGCGTGGGCCt                                                                                   >  2:229668/1‑148 (MQ=255)
                                                                 gCCGCAGAAAAAGCCCGACCCGGCGGCGTTGCTGTTCGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTTGTCGGTGATTCGCGCAGTGACGTGCAGGCGGCCAAGGCGGCTGGCGTGCAATGCGTGGGCCTGACCTAc                                                                            <  6:198250/148‑1 (MQ=255)
                                                                 gCCGCAGAAAAAGCCCGACCCGGCGGCGTTGCTGTTCGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTTGTCGGTGATTCGCGCAGTGACGTGCAGGCGGCCAAGGCGGCTGGCGTGCAATGCGTGGGCCTGACCTAc                                                                            <  6:79972/148‑1 (MQ=255)
                                                                    gCAGAAAAAGCCCGACCCGGCGGCGTTGCTGTTCGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTTGTCGGTGATTCGCGCAGTGACGTGCAGGCGGCCAAGGCGGCTGGCGTGCAATGCGTGGGCCTGACCTACGGc                                                                         >  7:272650/1‑148 (MQ=255)
                                                                     cAGAAAAAGCCCGACCCGGCGGCGTTGCTGTTCGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTTGTCGGTGATTCGCGCAGTGACGTGCAGGCGGCCAAGGCGGCTGGCGTGCAATGCGTGGGCCTGACCTACGGCta                                                                       >  6:79091/1‑149 (MQ=255)
                                                                      aGAAAAAGCCCGACCCGGCGGCGTTGCTGTTGGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTTGTCGGTGATTCGCGCAGTGACGTGCAGGCGGCCAAGGCGGCTGGCGTGCAATGCGTGGGCCTGACCTACGGCtat                                                                      <  1:158931/149‑1 (MQ=255)
                                                                           aaGCCCGACCCGGCGGCGTTGCTGTTCGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTTGTCGGTGATTCGCGCAGTGACGTGCAg                                                                                                                      >  6:292961/1‑96 (MQ=255)
                                                                             gCCCGACCCGGCGGCGTTGCTGTTCGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTTGTCGGTGATTCGCGCAGTGACGTGCAGGCGGCCAAGGCGGCTGGCGTGCAATGCGTGGGCCTGACCTACGGCTATAACCATg                                                               <  4:28094/149‑1 (MQ=255)
                                                                                   cccGGCGGCGTTGCTGTTCGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTTGTCGGTGATTCGCGCAGTGACGTGCAGGCGGCCAAGGCGGCTGGCGTGCAATGCGTGGGCCTGACCTACGGCTATAACCATggccggc                                                         <  8:272650/149‑1 (MQ=255)
                                                                                   cccGGCGGCGTTGCTGTTCGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTTGTCGGTGATTCGCGCAGTGACGTGCAGGCGGCCAAGGCGGCTGGCGTGCAATGCGTGGGCCTGACCTACGGCTATAACCATggccggc                                                         <  1:229668/149‑1 (MQ=255)
                                                                                      ggcggcGTTGCTGTTCGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTTGTCGGTGATTCGCGCAGTGACGTGCAGGCGGCCAAGGCGGCTGGCGTGCAATGCGTGGGCCTGACCTACGGCTATAACCATGGCCGGCCGa                                                      <  3:289882/149‑1 (MQ=255)
                                                                                                     cGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTTGTCGGTGATTCGCGCAGTGACGTGCAGGCGGCCAAGGCGGCTGGCGTGCAATGCGTGGGCCTGACCTACGGCTATAACCATGGCCGGCCGATCGACGCCGAATccc                                       >  5:205509/1‑149 (MQ=255)
                                                                                                         aTGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTTGTCGGTGATTCGCGCAGTGACGTGCAGGCGGCCAAGGCGGCTGGCGTGCAATGCGTGGGCCTGACCTACGGCTATAACCAt                                                                >  6:175343/1‑120 (MQ=255)
                                                                                                             aGATGGCCGGTGTCACCCCACAGCAGTCCCTGTTTGTCGGTGATTCGCGCAGTGACGTGCAGGCGGCCAAGGCGGCTGGCGTGCAATGCGTGGGCCTGACCTACGGCTATAACCATGGCCGGCCGATCGACGCCGAATCCCCCAGCCTg                               <  3:54066/149‑1 (MQ=255)
                                                                                                             aGATGGCCGGTGGCACCCCGCAGCAGGCCCTGGTGGGCGGTGATGCGCGCAGTGACGTGCAGCCGGCCAAGGCGGCGGGCGTGCATTGCGTGGGTCTGACCTACGGCTATAACCATGGCCGGCCGATCGACGCCGAATCCCCCAGCCTg                               <  1:259692/149‑1 (MQ=255)
                                                                                                                 ggCCGGTGTCACCCCGCAGCAGTCCCTGTTTGTCGGTGATTCGCGCAGTGACGTGCAGGCGGCCAAGGCGGCTGGCGTGCAATGCGTGGGCCTGACCTACGGCTATAACCATGGCCGGCCGATCGACGCCGAATCCCCCAGCCTGGTaa                           >  6:89526/1‑149 (MQ=255)
                                                                                                                  gCCGGTGTCACCCCGCAGCAGTCCCTGTTTGTCGGTGATTCGCGCAGTGACGTGCAGGCGGCCAAGGCGGCTGGCGTGCAATGCGTGGGCCTGACCTACGGCTATAACCATGGCCGGCCGATCGACGCCGAATCCCCCAGCCTGGTAAt                          <  5:89526/149‑1 (MQ=255)
                                                                                                                   ccGGTGTCACCCCGCAGCAGTCCCTGTTTGTCGGTGATTCGCGCAGTGACGTGCAGGCGGCCAAGGCGGCTGGCGTGCAATGCGTGGGCCTGACCTACGGCTATAACCATGGCCGGCCGATCGACGCCGAATCCCCCAGCCTGGTAATc                         >  7:154939/1‑149 (MQ=255)
                                                                                                                           aCCCCGCAGCAGTCCCTGTTTGTCGGTGATTCGCGCAGTGACGTGCAGGCGGCCAAGGCGGCTGGCGTGCAATGCGTGGGCCTGACCTACGGCTATAACCCTGGCCGGCCGATCGACGCCGAATCCCCCAGCCTGGTAATCGCCGACCt                 >  8:278011/1‑149 (MQ=255)
                                                                                                                           aCCCCGCAGCAGTCCCTGTTTGTCGGTGATTCGCGCAGTGACGTGCAGGCGGCCAAGGCGGCTGGCGTGCAATGCGTGGGCCTGACCTACGGCTATAACCATGGCCGGCCGATCGACGCCGAATCCCCCAGCCTGGTAAt                          >  7:147888/1‑140 (MQ=255)
                                                                                                                           aCCCCGCAGCAGTCCCTGTTTGTCGGTGATTCGCGCAGTGACGTGCAGGCGGCCAAGGCGGCTGGCGTGCAATGCGTGGGCCTGACCTACGGCTATAACCATGGCCGGCCGATCGACGCCGAATCCCCCAGCCTGGTAAt                          <  8:147888/140‑1 (MQ=255)
                                                                                                                           aCCCCGCAGCAGTCCCTGTTTGTCGGTGATTCGCGCAGTGACGTGCAGGCGGCCAAGGCGGCTGGCGTGCAATGCGTGGGCCTGACCTACGGCTATAACCATGGCCGGCCGATCGACGCCGAATCCCCCAGCCTGGTAATCGACGACCt                 >  1:179926/1‑149 (MQ=255)
                                                                                                                             cccGCAGCAGTCCCTGTTTGTCGGTGATTCGCGCAGTGACGTGCAGGCGGCCAAGGCGGCTGGCGTGCAATGCGTGGGCCTGACCTACGGCTATAACCATGGCCGGCCGATCGACGCCGAATCCCCCAGCCTGGTAATCGACGACCTgc               >  4:2841/1‑149 (MQ=255)
                                                                                                                                gcagcaGTCCCTGTTTGTCGGTGATTCGCGCAGTGACGTGCAGGCGGCCAAGGCGGCTGGCGTGCAATGCGTGGGCCTGACCTACGGCTATAACCATGGCCGGCCGATCGACGCCGAATCCCCCAGCCTGGTAATCGACGACCTgcgcg            >  7:52620/1‑149 (MQ=255)
                                                                                                                                  agcagTCCCTGTTTGTCGGTGATTCGCGCAGTGACGTGCAGGCGGCCAAGGCGGCTGGCGTGCAATGCGTGGGCCTGACCTACGGCTATAACCATGGCCGGCCGATCGACGCCGAATCCCCCAGCCTGGTAATCGACGACCTGCGCGCa          >  7:62624/1‑149 (MQ=255)
                                                                                                                                    cagTCCCTGTTTGTCGGTGATTCGCGCAGTGACGTGCAGGCGGCCAAGGCGGCTGGCGTGCAATGCGTGGGCCTGACCTACGGCTATAACCATGGCCGGCCGATCGACGCCGAATCCCCCAGCCTGGTAATCGACGACCTGCGCGCAtt        >  3:218356/1‑149 (MQ=255)
                                                                                                                                         ccTGTTTGTCGGTGATTCGCGCAGTGACGTGCAGGCGGCCAAGGCGGCTGGCGTGCAATGCGTTGGCCTGACCTACGGCTATATCCATGGCCGGCCGATCGACGCCGAATCCCCCAGCCTGGTAATCGACGACCTGCGCGCATtgctgc   <  5:235103/149‑1 (MQ=255)
                                                                                                                                         ccTGTTTGTCGGTGATTCGCGCAGTGACGTGCAGGCGGCCAAGGCGGCTGGCGTGCAATGCGTGGGCCTGACCTACGGCTATAACCATGGCCGGCCGATCGCCGCCGAATCCCCCAGCCTGGTAATCGACGACCTGCGCGCATtgctgc   >  8:291517/1‑149 (MQ=255)
                                                                                                                                          cTGTTTGTCGGTGATTCGCGCAGTGACGTGCAGGCGGCCAAGGCGGCTGGCGTGCAATGCGTGGGCCTGACCTACGGCTATAACCATGGCCGGCCGATCGACGCCGAATCCCCCAGCCTGGTAATCGACGACCTGCGCGCATTGCTGcc  >  3:133270/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                               
GGCCCGCTGCTGGACCAGATGAAGATCGGTAATTTCTTCCGCTGGATCATTGGTGGCGACACCTTGCCGCAGAAAAAGCCCGACCCGGCGGCGTTGCTGTTCGTCATGCAGATGGCCGGTGTCACCCCGCAGCAGTCCCTGTTTGTCGGTGATTCGCGCAGTGACGTGCAGGCGGCCAAGGCGGCTGGCGTGCAATGCGTGGGCCTGACCTACGGCTATAACCATGGCCGGCCGATCGACGCCGAATCCCCCAGCCTGGTAATCGACGACCTGCGCGCATTGCTGCC  >  NC_002947/505012‑505298

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