Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 72,361 G→A S189S (AGC→AGT thiQ ← thiamine ABC transporter ATP binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091372,3610GA98.3% 182.7 / NA 59S189S (AGC→AGTthiQthiamine ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (39/19);  minor base C (1/0);  total (40/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCACGCTGTGCGACACCATCAATAGCGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAACATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGTACCAGACAACGCGCTAACGCCACTCGCTGTCGCTGACCGCCGGAAAGCTCGCCCGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGA  >  NC_000913/72071‑72653
                                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                                    
aTGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGac                                                                                                                                                                                                                                                                                             <  1:121815/300‑1 (MQ=255)
  gCCGTTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGTGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGacac                                                                                                                                                                                                                                                                                           <  2:327063/300‑1 (MQ=255)
  gCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTATGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGacac                                                                                                                                                                                                                                                                                           <  1:517084/300‑1 (MQ=255)
  gCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGacac                                                                                                                                                                                                                                                                                           <  1:429736/300‑1 (MQ=255)
  gCCGGTGGATATCAATATTTTGCGTAAAGGGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGTTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGacac                                                                                                                                                                                                                                                                                           <  1:517093/300‑1 (MQ=255)
  gCCGGTATATATCAATATTTTTCGTAAATTGGTTGTGGTTTAGTGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTNACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTCCCAGTCGATGCGCCCATCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGacac                                                                                                                                                                                                                                                                                           <  2:222386/300‑1 (MQ=255)
                       gCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCAtctttt                                                                                                                                                                                                                                                                       >  2:640463/1‑299 (MQ=255)
                                             gggattgTCTTGTGCTGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCAcc                                                                                                                                                                                                                                                <  1:202351/294‑1 (MQ=255)
                                             gggattgTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGTCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCAcc                                                                                                                                                                                                                                                <  1:202130/294‑1 (MQ=255)
                                                   gTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGacac                                                                                                                                                                                                                                                                                           >  2:464620/1‑251 (MQ=255)
                                                   gTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGacac                                                                                                                                                                                                                                                                                           <  1:464620/251‑1 (MQ=255)
                                                   gTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTc                                                                                                                                                                                                                                          >  1:555664/1‑300 (MQ=255)
                                                                    tGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGt                                                                                                                                                                                                                           <  1:625268/298‑1 (MQ=255)
                                                                       tGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGc                                                                                                                                                                                                                      <  1:79377/300‑1 (MQ=255)
                                                                                     cgTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTACTGACGTAACGCCGGATCGAgcgaa                                                                                                                                                                                                         >  2:267130/1‑297 (MQ=255)
                                                                                          tATCAGGCCTCCAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCACAATAGTGCTGAAGCACTCGCCTTACCGCTCAACACCTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTATGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCATCGTCagca                                                                                                                                                                                                                                      >  2:272352/1‑262 (MQ=255)
                                                                                          tATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTACTGACGTAACGCCGGATCGAGCGCAGAGAAc                                                                                                                                                                                                   >  2:271928/1‑300 (MQ=255)
                                                                                               ggCCTACAAAGCGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCACAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTAACCTGCCAGGCGATGCGCCCGTCGGCGACCACCAACGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTACTGACGTAACGCAGGATCGAGCGCAGAGAACGGTTc                                                                                                                                                                                              >  2:343824/1‑300 (MQ=255)
                                                                                               ggCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGTGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTc                                                                                                                                                                                              >  1:430546/1‑300 (MQ=255)
                                                                                               ggCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTACTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTc                                                                                                                                                                                              >  1:430551/1‑300 (MQ=255)
                                                                                                                gtccGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTACTGACGTAACGCAGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATc                                                                                                                                                                              >  2:246201/3‑299 (MQ=255)
                                                                                                                                             gcATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTACTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAAc                                                                                                                                                >  2:552149/1‑300 (MQ=255)
                                                                                                                                               agctgACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCACTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGc                                                                                                                                              <  2:616651/296‑1 (MQ=255)
                                                                                                                                                  tgACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTANCGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCAc                                                                                                                                           >  1:503796/2‑300 (MQ=255)
                                                                                                                                                       gTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACgcgct                                                                                                                                      >  2:578712/1‑300 (MQ=255)
                                                                                                                                                       gTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACgcgct                                                                                                                                      >  1:599132/1‑300 (MQ=255)
                                                                                                                                                       gTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTACTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACgcgct                                                                                                                                      >  2:483914/1‑300 (MQ=255)
                                                                                                                                                                     gTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAGTGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGAccgc                                                                                                                          <  2:405792/298‑1 (MQ=255)
                                                                                                                                                                        aTCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAANGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGaa                                                                                                                     >  1:145352/1‑300 (MQ=255)
                                                                                                                                                                          ccccAATAGTGCTGAANCACTCTCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCCCACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTAGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACTAGACAACNCGCTAACGCCACGCGCTGNCGCTGACCGCCGGaaa                                                                                                                    >  1:581160/1‑299 (MQ=255)
                                                                                                                                                                          ccccAATAGTGCTGAANCACTCTCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACNCGCTAACGCCACGCGCTGNCGCTGACCGCCGGaaa                                                                                                                    >  1:581092/1‑299 (MQ=255)
                                                                                                                                                                          ccccAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCAACAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGaaa                                                                                                                    >  1:240087/1‑299 (MQ=255)
                                                                                                                                                                           gcgagatgtgtataagagacagCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGCCGACGACCACCAACGAGCGAGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGTTGCTGGCAACTCGTGCTCACCAGCGTCAGCATCCCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATAGGGCAATaaaa                                                                                                                                                                                >  2:399618/22‑238 (MQ=255)
                                                                                                                                                                                  gTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGNATCTTCCACACTGTGTGAAACCATCAATAACGNCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTAAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGAAAACGAGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGcca                                                                                                            >  2:513315/1‑298 (MQ=255)
                                                                                                                                                                                  gTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCAAACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATTGGCTGTTCGCGCACCAGCCAACGCGCTAACGCAACGTGCTGTCGCTGACCGCCGGAAAGCTCGcca                                                                                                            >  2:60019/1‑298 (MQ=255)
                                                                                                                                                                                  gTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTAACCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTACTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAACAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGcca                                                                                                            >  2:59394/1‑298 (MQ=255)
                                                                                                                                                                                     ctgAAGCACTCGTCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCAGCGACCACCACCGAGCGCGTGGCGATCCGCGTGGCATCTTCCACACTGTGTGACGCCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGCTGATCGAGCAATaa                                                                                                                                                                                  >  2:59489/3‑226 (MQ=255)
                                                                                                                                                                                                 ccTTACCGCTCAACAACTCATTGGTCTTACCCTGCCGGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGAAGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGAGCTAACGCCACGCGCTGTAGCTGACCGCAGGAAAGCTCGCCAGGTAACCGTGCCATTa                                                                                            >  2:401310/1‑300 (MQ=255)
                                                                                                                                                                                                   ttACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTaaa                                                                                          >  1:423995/1‑300 (MQ=255)
                                                                                                                                                                                                            aacaacTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAAt                                                                                  <  1:267130/299‑1 (MQ=255)
                                                                                                                                                                                                                 cTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCTCGCCATTAAATTATCAATCCCCAt                                                                            >  1:230412/1‑300 (MQ=255)
                                                                                                                                                                                                                  tCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATc                                                                           >  1:394166/1‑300 (MQ=255)
                                                                                                                                                                                                                        gtcttACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTggc                                                                       >  1:194273/5‑298 (MQ=255)
                                                                                                                                                                                                                        gtcttACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTACTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTggc                                                                       >  1:194311/5‑298 (MQ=255)
                                                                                                                                                                                                                                       gCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCAACAGCGTCAGCATCTACTGACGTNACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCAtt                                                      >  2:504134/1‑300 (MQ=255)
                                                                                                                                                                                                                                            gcgcCCGTCTGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGTCGGGCGATAGCGTGCATTTTctc                                                 <  2:394166/300‑3 (MQ=255)
                                                                                                                                                                                                                                              gcCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTACTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCct                                               >  1:384812/1‑300 (MQ=255)
                                                                                                                                                                                                                                                  gTCGGCGACCACCACCGAGCGCGTGGCAATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACTTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAAGAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTCCTCctgc                                             <  2:470104/298‑1 (MQ=255)
                                                                                                                                                                                                                                                     ggCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGtgc                                        <  2:563704/300‑3 (MQ=255)
                                                                                                                                                                                                                                                        gcccacccacGAGCGCGTGGCGATCCGCGCGGCATCTTCCACCCTGTGTGACACCATCAATAACGTCCTTTTTTGCTGCTGGCCGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGACGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTCACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGAAAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCATGCTGTACcgc                                     >  2:427590/10‑300 (MQ=255)
                                                                                                                                                                                                                                                        gaccaccacCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACTTCATTTTTTGCTGCTTGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCcg                                      <  1:343824/292‑1 (MQ=255)
                                                                                                                                                                                                                                                        gaccaccacCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCcg                                      <  1:246201/292‑1 (MQ=255)
                                                                                                                                                                                                                                                        gaccaccacCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACCGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCctg                                              >  2:427584/8‑291 (MQ=255)
                                                                                                                                                                                                                                                                                         aTCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGATAACGCCAAGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCAATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCAGGATTTAGCCCCAGca     >  2:35968/1‑298 (MQ=255)
                                                                                                                                                                                                                                                                                           cTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGa  >  1:169332/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                           cTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTACTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGa  >  2:399109/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                           cTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGACAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGa  >  2:163959/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                           cTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTACTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGa  >  1:366539/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                           cTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTACTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCAATCTGACGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCCGccc    >  2:240445/1‑298 (MQ=255)
                                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                                    
ATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCACGCTGTGCGACACCATCAATAGCGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAACATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGTACCAGACAACGCGCTAACGCCACTCGCTGTCGCTGACCGCCGGAAAGCTCGCCCGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGA  >  NC_000913/72071‑72653

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