Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 74,046 A→G 18.0% L159P (CTC→CCC)  thiP ← thiamine/thiamine pyrophosphate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091374,0460AG18.0% 173.1 / 2.9 72L159P (CTC→CCC) thiPthiamine/thiamine pyrophosphate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (13/43);  new base G (12/2);  total (27/45)
Fisher's exact test for biased strand distribution p-value = 2.79e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCG  >  NC_000913/73810‑74283
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                               
aGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCGAGTTGGCGCAGGTCGCCGGGGATGTTTTCCAGGGCCTGGGGt                                                                                                                                                                                                                                               >  1:454307/1‑239 (MQ=255)
          cAGCGCCTGAAAGATTGCCAGCGCGAAAGTGGTCGCCTGCGGACCGCCCCGCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGGGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGc                                                                                                                                                                                                                                     <  1:68699/239‑1 (MQ=255)
           aGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTCACCACCGCCAATCGCCGACACCGAAAAAATGCAAGCTACGCCTCCCAACCTAGGCGGCAAGTTGACGCTGTTCGCCGCGGGAGTTTTTCCAGGCCTGGGGgaataaaaggc                                                                                                                                                                                                                                     >  2:541477/1‑233 (MQ=255)
           aGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGACGGGATTTGTCGCCGCAACCACGGCCATTCGACGAAGCGGAAAAAATACCAGCTCCGCATCCCAAGCAGGGCGGCAAGGTGACGCCATTCGCCGGGGGTGGTTTCCAGGGCAGGGAGtacaaaggggat                                                                                                                                                                                                                                    >  2:68699/1‑229 (MQ=255)
                   aTAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACCCGGTTCCCCGGGGAGGTTTTCCAGTGCCTGGGGTAATAAACGGCCCCCCATCg                                                                                                                                                                                                                            >  1:336906/1‑239 (MQ=255)
                     aGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCgg                                                                                                                                                                                                                           <  2:919930/238‑1 (MQ=255)
                      gATTGCCAGCTCGATAGGGGTCCCCGGGCGACCGCCCCCCAGCGATAGCCCGGTGGCGAAGCTGGCAAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCa                                                                                                                                                                                                                         <  1:52812/239‑1 (MQ=255)
                      agggacaggCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCa                                                                                                                                                                                                                         <  1:841594/231‑1 (MQ=255)
                              gCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCCTAACCACGGCCATTCGCCGAAGCGGCAAAAATGCCAGCTACGCATCCCCAGCTGGGCGGCAAGTTGACGCTGTTTCCCCGGGGTGTTTTCCCGTGCCCGGGGAAAAAAGCGGCGGCCCACCGGCAGAGGaaaa                                                                                                                                                                                                                >  2:52812/1‑239 (MQ=255)
                              gCTCGATAGGGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGGGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTaaaa                                                                                                                                                                                                                <  1:1085449/239‑1 (MQ=255)
                               cTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCCAGTTGACGCTGTTCGCCGGGGAGGTTTTCCAGTGCCTGGGGTAATAAGCGGCTCGCCATCGGCAGATTaaaaa                                                                                                                                                                                                               >  1:179273/1‑239 (MQ=255)
                                     tggggcgcccTGCGGGCCCCCCCCCCGCGATGGCCCGGTGGCGAAGCCGGGTAAAAGGAGCAAAAAGTATAGCGCAGCAACCGGCGGGATTTGGCGCCGTAACCCCGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAAcaca                                                                                                                                                                                                          <  2:142454/230‑1 (MQ=255)
                                     aggggccggcTGCGGCCCCCCCCCCAGCGATAGCACGGTGGCGAAGGTGGCGAAACAGAGCATAAAGATAAGCCCAGCAACCGGCGGGATTTGTCGCCGTAACCCCGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAAcaca                                                                                                                                                                                                          <  1:1092837/229‑1 (MQ=255)
                                     aGTGGCCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAAcaca                                                                                                                                                                                                          <  1:406512/238‑1 (MQ=255)
                                         gtccccgcgggACCGCCCCCCGGCAGTAGCCCAGTGGCGAAGTTGGCGAACCAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGTCCATTCGACGAAGCGGAAAAAATGCCCGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGc                                                                                                                                                                                                     <  2:372471/230‑1 (MQ=255)
                                         gTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGGTTTCCAGTGCCTGGTGTCATAAGCGGCCCGCCATCGGCCGCTTTAAAAAaacactggc                                                                                                                                                                                                     >  1:435122/1‑234 (MQ=255)
                                         gTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGGAACCACGGCCATTCGACGAAGCGGAAAAAAAGCCAGCTACGCATCCCCAGCTGGGCGGCAAGTTGACGCTGTTTCCCGGGGAGGTTTTCCAGTGCCTGGGGAGAAAAACGGCCCCCCACCGGCAGGTTAAAAAaaacgtctct                                                                                                                                                                                                     >  2:841594/1‑233 (MQ=255)
                                          tCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGc                                                                                                                                                                                                     <  1:493083/238‑1 (MQ=255)
                                          gcGCCGGCGGCCCGCCCCCCAGCGATAACACGGTGGCGAAGGCGGTGAAACAGACGATAAAGAGAAGCGCAGCACACGGGGGGATTTGTCGCCGTAACCACGGCCATGTGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGTCGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGc                                                                                                                                                                                                     <  2:611590/237‑1 (MQ=255)
                                                 ccgACCGCCCCCCAGGGATAGCCCGGTGGCGAAGCTGGGGAAACAGAGACAAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCaaaa                                                                                                                                                                                            <  2:59781/237‑1 (MQ=255)
                                                 cGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCCGCTACGCATCCCAAGCTGGGCGGCAAGTTGACCCTGGTCGCCGGGGAGGTTTTCCCGAGCATGGGGTAATAAGCGGCCCCCCCTCCGCAGGTTTAAAAAACACAGGGCCCGAaaaa                                                                                                                                                                                            >  2:1035525/1‑239 (MQ=255)
                                                      cGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGCTGTTTTCCAGTGCCTGGAGGAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCt                                                                                                                                                                                       >  1:967334/1‑239 (MQ=255)
                                                         ccccgcgCGAGAGAACGGTGGCAAAGTGGGCTACACATAGCATAAAGAGACGCGCAGCACCGGGCGGGATTTGCCGCCGTAACCCGGGCCAGTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCGGGGGGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATGCCTTGc                                                                                                                                                                                    <  2:937788/233‑1 (MQ=255)
                                                                   tAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGGGTAATAAGCGGCTCGCCATCGGCAGCTTAAAAAACACATGGGCCAGCAAAAAACCTTTCAGGCCGTAgg                                                                                                                                                                          >  1:238295/1‑239 (MQ=255)
                                                                    tccgcGGGGGGGAAGGTGGCGAACCAGAGCAAAAAGATAAGCGCAGCGCCGGGGGGGGTTTGGAGCCCTACCCACTGGCAATCGACAAAGCGGAAAAAATGGCAGCGACGCATCCCAAGCTGGGGGGCAAGTGGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCATCAAAATACCGTGCAGGCCGTAgg                                                                                                                                                                          <  2:224511/234‑1 (MQ=255)
                                                                     gCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGGGTAATAAGCGGCTCGCCATCGGCAGGTTAAAAAAAACACGGGCCAGCAAAAAACCTTGCAGGCCGTAaggg                                                                                                                                                                        >  1:522274/1‑238 (MQ=255)
                                                                              ggaggCTGGGAAAACAGGGCATAAAGGAAAGCGCAGGAACCGGGGGGTTTTGTCGGCGTAACCCCGGCCATTCGACGAAGGGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGGCGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGt                                                                                                                                                                <  2:546229/234‑1 (MQ=255)
                                                                              ggAAGCTGGCGAAACAGGGAAAAAAGATAAGCGAAGAACCCGGCGGGATTTGTCGCCGTAACCCCGGCCATGCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGt                                                                                                                                                                <  2:1060563/237‑1 (MQ=255)
                                                                              ctaaGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGt                                                                                                                                                                <  1:192428/236‑1 (MQ=255)
                                                                              cGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGTATTTGTCGCCTTAACCACGGCCTTTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGt                                                                                                                                                                <  1:90226/238‑1 (MQ=255)
                                                                              cGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGt                                                                                                                                                                <  1:819234/238‑1 (MQ=255)
                                                                                aaGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAAGGCCAGCTACGCATCCCCAGCTGGGCGGCAAGTTGACGCTTTTCGCCGGGGATGTTTTCCCGGGCCTGGCGGGATAAGCGGCCCGCCATCGGGCGATTAAAAAAAACATCGGCGAGCAAAATACCTTGCAGGCCGAGGGGGGCAAAgaaca                                                                                                                                                              >  2:192428/1‑234 (MQ=255)
                                                                                 aGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCACGCAGGGCGGCCAGTTGACGCCGTTCGCCGGGGATGTTTTTCAGTGCCTGGGGTAATAACCGGCTCGCCACCGGCCGGTTAAAAAACACCAGGGCCAGCCAAAAAACATTCCGCGCCTAGGGGGAAAAggcacac                                                                                                                                                            >  2:1038925/1‑234 (MQ=255)
                                                                                     ggCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTTCAGGCCGGAAGGCGAAAAAGTCCACCCCa                                                                                                                                                        >  1:1056249/1‑239 (MQ=255)
                                                                                           aCAGAGCATAAAGATAAGCGAAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGa                                                                                                                                                  <  1:638188/239‑1 (MQ=255)
                                                                                                gCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGGGTAATAAGCGGCTCGCCATCGGCCGATTAAAAAACACATGGGCCAGCAAAATACCTTGCAGGCCGTAGGGGGAAAAGGTCCACTCCAGAccagccgat                                                                                                                                             >  1:305805/1‑232 (MQ=255)
                                                                                                gCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAAcaca                                                                                                                                                                                                          >  1:576557/1‑179 (MQ=255)
                                                                                                gCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAAcaca                                                                                                                                                                                                          <  2:576557/179‑1 (MQ=255)
                                                                                                gCATAAAGATAAGCGCAGAAGCCGGCGGGATTGGTCGCCGTACCCACGGCCATTCGACGAAGCGGAAAAAGAGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGc                                                                                                                                                <  1:1036678/236‑1 (MQ=255)
                                                                                                    aaaGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGc                                                                                                                                         <  1:1004320/239‑1 (MQ=255)
                                                                                                           aGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAAAgggt                                                                                                                                  >  1:606889/1‑239 (MQ=255)
                                                                                                                    aCCGGGGGGATTTGGCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTgccagcca                                                                                                                         <  1:284452/239‑1 (MQ=255)
                                                                                                                         cggGTTTTCTCGCCGTAACCACGGCCATTCGCCGAAGCGGAAAAAATGCCAGCTGCGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCt                                                                                                                    <  2:1106928/239‑1 (MQ=255)
                                                                                                                          gggATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTGCAGGCCGAAGGGCGAAAAGGTCCACCCCAAACCGAGCGATTGGCAGAGTGTTGCCAGCCAGcccgg                                                                                                                   >  1:745890/1‑237 (MQ=255)
                                                                                                                              tttGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCg                                                                                                                <  1:435599/238‑1 (MQ=255)
                                                                                                                              tttGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCg                                                                                                                <  1:217155/238‑1 (MQ=255)
                                                                                                                              tttGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGGTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGGGTAATAAGCGGCTCGCCATCGGCCGATTAAAAAACACATGGGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGGACACTCCAAACCACACGACTGGCAGAGGGTTTCCACCCAACCCCgggc                                                                                                                >  2:45298/1‑236 (MQ=255)
                                                                                                                                   cGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATa                                                                                                          <  1:630110/239‑1 (MQ=255)
                                                                                                                                   cGCCCTAACCCCGGCCTTTCGGCGAGGCGGAAAAAATGCCAGCTACGCTTCCCAAGCGGGGGGGCGAGTGGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACCCATGGGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATa                                                                                                          <  2:454307/239‑1 (MQ=255)
                                                                                                                                                gCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGc                                                                                                                                                                                                     >  2:302483/1‑136 (MQ=255)
                                                                                                                                                gCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGc                                                                                                                                                                                                     <  1:302483/136‑1 (MQ=255)
                                                                                                                                                    ttCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGa                                                                                         >  1:27558/1‑239 (MQ=255)
                                                                                                                                                     tCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTTCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAAACCATCGGCCAGCAAAATACCTTGACGGCCGCAGGGCGAAAAGGGCCACTCCAGACCGAGCGATTGGACGAGGGTTGTCCCCCCGCCCTGGCGACCACAGACACCAAGaagaccaaa                                                                                        >  2:585493/1‑232 (MQ=255)
                                                                                                                                                      cgacgaAGCGAAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGaaa                                                                                       >  1:776337/1‑239 (MQ=255)
                                                                                                                                                         cgaAGGGGAAAAAATGCCGGCTACGATCCCCAAGGGGGGGGGAAAGTTTGCGGTGTTCGCCGGGGAAGTTTACCAGTGCCTGGAGTAATAAGCGGCTCGCCAGCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACa                                                                                    <  2:435122/239‑1 (MQ=255)
                                                                                                                                                               ggAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGAc                                                                               <  1:45298/238‑1 (MQ=255)
                                                                                                                                                               ggAAAAAAGGCCAGCTACCAATCCCAAGTGGGGCGGGAAGTTGACGGTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGAc                                                                               <  2:606889/238‑1 (MQ=255)
                                                                                                                                                                  aaaaaGGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGAGGATGTTCGCCGGGGATGTTTTCCAGTGCCGGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCaa                                                                            <  1:1038925/238‑1 (MQ=255)
                                                                                                                                                                        gCCAGTTACGCATCCCAAGTGGGGCGGAAAGTTGACGTGGTTCGCCGGGGATGTTTTCCAGGGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCggg                                                                     <  2:776337/239‑1 (MQ=255)
                                                                                                                                                                          ctggtACGCATCCCAAGCTGGGCGGCAAGTGGACGCTGTTGGCCGGGAATGCTTTCCAGTGCCTGGATTAATAAGCGGCTCGCCATCGGCAGATTAAAAGACACATGGGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGgga                                                                    <  1:1064436/234‑1 (MQ=255)
                                                                                                                                                                          cAGCTCCGCATCCAAAACTGGGTGTCAAGGTGACCCTGTTCGCCGGGGATGTTTTCCAGTGACTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACTTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGgga                                                                    <  1:913933/238‑1 (MQ=255)
                                                                                                                                                                              tACGCCTCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATTCCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGAt                                                                <  1:1103626/238‑1 (MQ=255)
                                                                                                                                                                                            tGGGGGGCAAGTGAGCGCTGTTCGCCGGGGATGTTTTCACGGGCCGGGGAGAATAAGGGGCTCGCCATCGGGAGGTTAAAAAACAAATGGGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGcac                                                  <  2:336906/238‑1 (MQ=255)
                                                                                                                                                                                              gggggcAAGTGGGCGCTGTTCGCCGGGGTTTTTTTACAGTGCCTGGAGTAAGAACCTGCCCGCCAGCGGCAGATTAAAAAACCCATGGGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGcacac                                                <  2:305805/235‑1 (MQ=255)
                                                                                                                                                                                                                ttCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAAcaca                                                                                                                                                                                                          >  1:905321/1‑67 (MQ=255)
                                                                                                                                                                                                                ttCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAAcaca                                                                                                                                                                                                          <  2:905321/67‑1 (MQ=255)
                                                                                                                                                                                                                ttCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGc                                                                                                                 <  2:497201/155‑1 (MQ=255)
                                                                                                                                                                                                                ttCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGc                                                                                                                 >  1:497201/1‑155 (MQ=255)
                                                                                                                                                                                                                               tttCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAgcgc               <  1:895024/239‑1 (MQ=255)
                                                                                                                                                                                                                               tttCCAGTGCCGGGGGTAAAAAACGGCTCGCCATGGGAAAATTAAAAAACCCAGGGGCCAGCAAAATTCCTGGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAgcgc               <  2:522274/239‑1 (MQ=255)
                                                                                                                                                                                                                                ttcaggggCCTGGAGTAAAAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGGGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGcc              <  1:585493/232‑1 (MQ=255)
                                                                                                                                                                                                                                           ggAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAAACGCGACCAACACCGGGAGGAGCAAGGTCACTGCACCCAAACCAAAAAGCGCCACCCGCCCCCGGAAGCGCCTGCGaacaagc   >  2:304624/1‑233 (MQ=255)
                                                                                                                                                                                                                                             aGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATagaacc  >  1:725313/1‑235 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                               
AGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCG  >  NC_000913/73810‑74283

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: