Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 74,222 C→T L100L (TTG→TTA thiP ← thiamine ABC transporter membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091374,2220CT98.0% 168.1 / NA 50L100L (TTG→TTAthiPthiamine ABC transporter membrane subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (20/29);  minor base G (0/1);  total (20/30)
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

GCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCGAGAGAAACGCCTGCCAGAAGGAGAAGCGCACCACATGCCACAGATAGCTGTCCTGCCAGACTGCCACCCAGTCATCCTGCGGCGCGTTCCACCACAGGGCGAGAAACGCCGCCAGCGCAACCGCTACCACCAGCGTGGTGGC  >  NC_000913/73982‑74467
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     
gccgcGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTg                                                                                                                                                                                                                                             >  1:161782/5‑251 (MQ=255)
                ggggCGGCAAGCGGACGTTGTTCGCCGGGAATGTTGTCCAGTGCCTGGAGTAATAAGCGGCTAGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAg                                                                                                                                                                                                                             <  2:282270/250‑1 (MQ=255)
                 gggCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTACAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCAACGGGAGGCCTAAGGTCATTGCGCACAGACGCACCAgc                                                                                                                                                                                                                            >  1:296723/1‑251 (MQ=255)
                      gcaagtggacgtggtgagCCGGGAATGTTTTCAAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAAGACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACAGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAg                                                                                                                                                                                                                       <  2:95639/234‑1 (MQ=255)
                      gcAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACATGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCACAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAg                                                                                                                                                                                                                       <  2:153488/251‑1 (MQ=255)
                        aagctGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGATGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAGcc                                                                                                                                                                                                                     <  1:141042/247‑1 (MQ=255)
                        aagctGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAGcc                                                                                                                                                                                                                     >  2:20999/5‑251 (MQ=255)
                        aagcggacgtggTTCGCCGGGAAGGTTTTCCAGTGCCGGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAAGACGTGCGCCAGAAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACGCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAGcc                                                                                                                                                                                                                     <  2:104583/240‑1 (MQ=255)
                              acgttGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAGCCGACCCg                                                                                                                                                                                                               >  1:150944/5‑251 (MQ=255)
                                       gCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGcc                                                                                                                                                                                                      >  2:87074/1‑251 (MQ=255)
                                              aTGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGGAGAGTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACGCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGCCATTGCGCACAGACGCAACAGCGCCATCCGACCCGGAAAGCGCCTGCGATa                                                                                                                                                                                               <  2:183253/251‑1 (MQ=255)
                                                       aGTGCCTGGCGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATCCGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGATGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGcgcgcg                                                                                                                                                                                      >  1:258097/1‑251 (MQ=255)
                                                          gCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGAGTAAGGTCATTGCGCACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCggg                                                                                                                                                                                    <  2:179668/250‑1 (MQ=255)
                                                                           gCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGGGCGGGCGAGGAATATCGCGGGTa                                                                                                                                                                  <  1:251967/251‑1 (MQ=255)
                                                                                           ttaaaacatacgtgcgCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGGACGACAGAGAGCAgcgc                                                                                                                                                  >  2:71441/16‑248 (MQ=255)
                                                                                            tAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTACGACAGAGAGCAGCGcc                                                                                                                                                 >  1:173854/1‑251 (MQ=255)
                                                                                             aaaaaaTACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTACGACAGAGAGCAGCGCCg                                                                                                                                                <  1:20999/251‑1 (MQ=255)
                                                                                               aaaaatgttgcgCCAGCAAAAGTCCTAGCAGCCCGTAGGGCGAAAAGAGCCACACCAGCCCGAGCTATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCTCAGACGCGACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAAGATCGCGGGTACGACAGAGATCAGCGCCgag                                                                                                                                              <  2:191662/240‑1 (MQ=255)
                                                                                                aaatacgtgcgCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAGcc                                                                                                                                                                                                                     <  1:159191/169‑1 (MQ=255)
                                                                                                aaatacgtgcgCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTATTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAGcc                                                                                                                                                                                                                     >  2:159191/11‑179 (MQ=255)
                                                                                                  atacgcgcgCCAGCAAAATCCCTTGCACGACGTAGGGCGCAAAGGTCCACACCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCCTAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGCGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGGACGACAGATAGCAGCGCCgcg                                                                                                                                              >  2:135240/9‑246 (MQ=255)
                                                                                                    acgtgcgCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTACGACAGAGAGCAGCGCCGAGAGaaa                                                                                                                                         <  1:202371/245‑1 (MQ=255)
                                                                                                     cgtgcgCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTACGACAGAGAGCAGCGCCGAGAGAAAc                                                                                                                                        <  2:150944/246‑1 (MQ=255)
                                                                                                      gtgcgCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTACGACAGAGAGCAGCGCCGAGAGAAACg                                                                                                                                       <  1:2646/247‑1 (MQ=255)
                                                                                                           ccAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTACGACAGAGAGCAGCGCCGAGAGAAACGCCTGc                                                                                                                                  >  1:161175/1‑251 (MQ=255)
                                                                                                           ccAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTACGACAGAGAGCAGCGCCGAGAGAAACGCCTGc                                                                                                                                  >  1:207322/1‑251 (MQ=255)
                                                                                                                 aaaTACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTACGACAGAGAGCAGCGCCGAGAGAAACGCCTGCCAGa                                                                                                                              >  1:130772/1‑249 (MQ=255)
                                                                                                                  aagaCCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTACGACAGAGAGCAGCGCCGAGAGAAACGCCTGCCAGAAgg                                                                                                                           <  2:198416/248‑1 (MQ=255)
                                                                                                                     aCCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTACGACAGAGAGCAGCGCCGAGAGAAACGCCTGCCAGAAgg                                                                                                                           >  2:151470/1‑248 (MQ=255)
                                                                                                                      ccTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTACGACAGAGAGCAGCGCCGAGAGAAACGCCTGCCAGAAGgaga                                                                                                                        >  1:37861/1‑250 (MQ=255)
                                                                                                                                   aGGGCGAAAAGGGCCACGCCCGACCGGGCGATTGGCAGAGTGTGGCCAGCCAGCCCTGGCGAACATAGACGCTAAGGATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTACGACAGAGAGCAGCGCCGAGAGAAACGCCTGCCAGCAGTATAAGCGCACCACATGc                                                                                                          <  2:185914/251‑1 (MQ=255)
                                                                                                                                        gAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACTCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATGAAGGTCATTGCGCACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTACGACAGAGAGCAGCGCCGAGAGAAACGCCTGCCAGAAGGAGAAGCGCACCACATGCCACAg                                                                                                     <  2:289323/251‑1 (MQ=255)
                                                                                                                                              gTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTACGACAGAGAGCAGCGCCGAGAGAAACGCCTGCCAGAAGGAGAAGCGCACCACATGCCACAGATAGCt                                                                                               <  2:100498/251‑1 (MQ=255)
                                                                                                                                                  aCTCCAGACCGAGCGATTGGCAGAGGGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGACGAAAACAGCAAACAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGGACGACAGAGAGCAGCGCCGAGAGAAACGCCTGCCAGAAGGAGAAGCGCACCACATGCCACAGATAGCTGTcc                                                                                           >  1:283785/1‑251 (MQ=255)
                                                                                                                                                        gACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTACGACAGAGAGCAGCGCCGAGAGAAACGCCTGCCAGAAGGAGAAGCGCACCACATGCCACAGATAGCTGTCCTGCCAg                                                                                     <  2:215367/251‑1 (MQ=255)
                                                                                                                                                                 aTTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTACGACAGAGAGCAGCGCCGAGAGAAACGCCTGCCAGAAGGAGAAGCGCACCACATGCCACAGATAGCTGTCCTGCCAGACTGCCAcc                                                                            <  2:261145/251‑1 (MQ=255)
                                                                                                                                                                            gtTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTACGACAGAGAGCAGCGCCGAGAGAAACGCCTGCCAGAAGGAGAAGCGCACCACATGCCACAGATAGCTGTCCTGCCAGACTGCCACCCAGTCATCCt                                                                  <  1:320314/250‑1 (MQ=255)
                                                                                                                                                                              tGCCAGCCAGCCCTGGCGACCATAGCCGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCCCAGACGCAGCAGCGCCAGCCGACCCGTAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTCCGACAGAGAGCAGCGCCGATAGAAACGCCGGCCAGAAGGAGAAGCGCGCCACATGCCACAGATAGCTGGCCTGCCAGACTGCCCCCCAGTCATCCTgc                                                                <  2:141369/250‑1 (MQ=255)
                                                                                                                                                                               gccagccgCCCCCGGCGACCGGAGACGCTAAGAATGCCGAAAACAGCATCCGGCACCGGGAGGATTAAGGACGTTGCGCACAGCCGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGAGAGAGCGCGCGGGCGAGGAAGATCGCGGGTACGACAGAGAGCAGCGCGGAGAGAAACGCCTGCCAGAAGGAGAAGCGCACCACATGCCACAGATAGCTGTCCTGCCAGACAGCCACCCAGGCATCCTgcg                                                               <  2:47884/250‑1 (MQ=255)
                                                                                                                                                                                   gccagccCTGGGGAACATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTACGACAGAGAGCAGCGCCGAGAGAAACGCCTGCCAGAAGGAGAAGCGCACCACATGCCACAGATAGCTGTCCTGCCAGACTGCCACCCAGTCATCCTGCGGcgcg                                                          <  2:177129/251‑1 (MQ=255)
                                                                                                                                                                                         ccTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTACGACAGAGAGCAGCGCCGAGAGAAACGCCTGCCAGAAGGAGAAGCGCACCACATGCCACAGATAGCTGTCCTGCCAGACTGCCACCCAGTCATCCTGCGGCGCGTTccac                                                    <  1:100305/251‑1 (MQ=255)
                                                                                                                                                                                              cGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTACGACAGAGAGCAGCGCCGAGAGAAACGCCTGCCAGAAGGAGAAGCGCACCACATGCCACAGATAGCTGTCCTGCCAGACTGCCACCCAGTCATCCTGCGGCGCGTTCcaccac                                                 >  1:66413/1‑249 (MQ=255)
                                                                                                                                                                                                              gAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATGGCGCACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCGTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTACGACAGAGAGCAGCGCCGAGAGAAACGCCTGCCAGAAGGAGAAGCGCACCACATGCCACAGATAGCTGTCCTGCCAGACTGCCACCCAGTCATCCTGCGGCGCGTTCCACCACAGGGCGAGAAACgccgcc                               <  2:252097/251‑1 (MQ=255)
                                                                                                                                                                                                                 ggccgaaaccagcacacagcacCGGGGGGATTAAGGTCATTGCGCACAGCCGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTACGACAGAGAACAGCGCCGAGAGAAACGCCTGCCAGAAGGAGAAGCGCACCACATGCCACAGATAGCTGTCCTGCCAGACTGCCACCCAGTCATCCTGCGGCGCGTTCCACCCCAGGGCGAGAAACGCCGCCAgc                            <  2:219136/232‑1 (MQ=255)
                                                                                                                                                                                                                     gAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTACGCCAGAGAGCAGCGCCGAGAGAAACGCCTGCCAGAAGGAGAAGCGCACCACATGCCACAGATAGCTGTCCTGCCAGACTGCCACCCAGTCATCCTGCGGCGCGTTCCACCACAGGGCGAGAAACGCCGCCAGCGCa                         >  1:132333/1‑250 (MQ=255)
                                                                                                                                                                                                                                 cagcacCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTACGACAGAGAGCAGCGCCGAGAGAAACGCCTGCCAGAAGGAGAAGCGCACCACATGCCACAGATAGCTGTCCTGCCAGACTGCCACCCAGTCATCCTGCGGCGCGTTCCACCACAGGGCGAGAAACGCCGCCAGCGCAACCGCTaccaccc            >  2:25392/4‑250 (MQ=255)
                                                                                                                                                                                                                                  cagcacGGGAGGGATTAAGGTCATTGGTCACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCGGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTACGACAGAGAGCAGCGCCGAGAGAAACGCCTGCCAGAAGGAGAAGCGCACCACATGCCACAGATAGCTGGCCTGCCATACTGCCACCCAGTCAGCCTGCGGCGCGTTCCACCACAGGGCGAGAAACGCCGCAAGCGCAACCGCTAccacca            <  2:182573/246‑1 (MQ=255)
                                                                                                                                                                                                                                   ccacCGGGAGGATTAAGGTCAGTGCGCACAGACGCAACAGCGCCAGCCGACCCGAAAAGCGCCTGCGATAGGGCGCGCGGGCGAGGAATATCGCGGGTACGACAGAGAGCAGCTCCGAGAGAAACGCCTGCCAGAAGGAGAAGCGCACCACAGTCCACAGATAGCTGTCCTGCCAGACTGCCACCCAGTCATCCTGCGGCGCGTTCCACCACAGGGCGAGAAACGCCGCCAGCGCAACCGCTACCTCCAGc          <  2:45368/250‑1 (MQ=255)
                                                                                                                                                                                                                                          gaggaTTAAGGTCATTGCGCACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTACGACAGAGAGCAGCGCCGAGAGAAACGCCTGCCAGAAGGAGAAGCGCACCACATGCCACAGATAGCTGTCCTGCCAGACTGCCACCCAGTCATCCTGCGGCGCGTTCCACCACAGGGCGAGAAACGCCGCCAGCGCAACCGCTACCACCAGCGtggtgg   <  1:239759/251‑1 (MQ=255)
                                                                                                                                                                                                                                           aggaTTAAGGTCATTGCGCACAGACGCAACAGCTCCAGCCGACCCGGACAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTACGACAGAGAGCAGCGCCGAGAGAAACGCCTGCCAGAAGGAGAAGCGCACCACATGCCACAGATAGCTGACCTGCCAGACTGCCACCCCGTCATCCTGCGGCGCGTTCCACCACAGGGCGAGAAACGCCGCCAGCGCAACCGCTACCACCAGCGTGGTGgc  >  2:227045/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     
GCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCGAGAGAAACGCCTGCCAGAAGGAGAAGCGCACCACATGCCACAGATAGCTGTCCTGCCAGACTGCCACCCAGTCATCCTGCGGCGCGTTCCACCACAGGGCGAGAAACGCCGCCAGCGCAACCGCTACCACCAGCGTGGTGGC  >  NC_000913/73982‑74467

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