Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	72,928	T→C	N532D (AAT→GAT)	thiP ←	thiamine ABC transporter membrane subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	72,928	0	T	C	97.1%	211.8 / ‑5.8	68	N532D (AAT→GAT)	thiP	thiamine ABC transporter membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base C (32/34);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGTGTGATCTACGCCATCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTT  >  NC_000913/72684‑73172                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       ggTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGAGTGATCACCGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCCTCGGCAAATGTTGGTAAAGCCCAGTGGTATCAGTCAGTTTTAACATTTCGccc                                                                                                                                                                                                                                                  >  2:425865/1‑251 (MQ=255) ggTTGTTGTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGccc                                                                                                                                                                                                                                                  <  1:291200/251‑1 (MQ=255)  gttgttTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGcccc                                                                                                                                                                                                                                                 <  1:44122/251‑1 (MQ=255)       tttcCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTa                                                                                                                                                                                                                                             <  1:207164/248‑1 (MQ=255)        ttcCTGAAACATCATCGACACCGGACGGAGTGACGGCGGGGTCGTGGTGTGAACAACGCCGTCGATAACCAGCGATCCGCTGGCTGGCGCCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGttg                                                                                                                                                                                                                                            <  2:524112/249‑3 (MQ=255)        ttcCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAg                                                                                                                                                                                                                                            <  1:83965/249‑1 (MQ=255)                ccAGCATCGACACAGGACGGCGTGCCGGCTGCGTCGTGGTGTGGTCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGAACAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTTAGCTAAAACGCATCGGTCAATGTTGGTAAAGCCAGGAGCTATCAGTCAGTTTTAACCTTTCGCCCCGGTAATTTTTccac                                                                                                                                                                                                                                   >  2:301675/2‑248 (MQ=255)                aCAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAAt                                                                                                                                                                                                                                   >  1:699840/1‑251 (MQ=255)                 cAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGAAATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATc                                                                                                                                                                                                                                  >  1:234547/1‑251 (MQ=255)                  aGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAATAAGCGAACCGCTGGCTGGCGTCAGAAAACAAGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGATCGACGCGTTCTACCGTTAAGCTAAAACGCATCGGCAACTGTTGGTAAAGACAGGTGATATCAGTCAGGTTTAACATCTCGCCCCGGTAGTTTTTCAAta                                                                                                                                                                                                                                  >  1:703972/1‑249 (MQ=255)                     atcgaacccgtccGGCGTGACGGCGCGATCGTGGTGTGATCAAGGCCGTCGATAGTCAGCGAGCCGCGGGCTGGCGTCATAAACCCGTCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCGTGCCGCGTGCTACCGTTAAGCGAAAACGCAGCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTGTTCAATCACgg                                                                                                                                                                                                                              <  2:491846/239‑1 (MQ=255)                           cccGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAAca                                                                                                                                                                                                                        <  1:530579/250‑1 (MQ=255)                            ccGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAAcag                                                                                                                                                                                                                       <  2:454654/251‑1 (MQ=255)                            ccGGACGGCGTGACGGCGGCGTCGTGGTGTGAGCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTGCAATCACGGTAAAcag                                                                                                                                                                                                                       <  1:219768/251‑1 (MQ=255)                                aCGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAAcagc                                                                                                                                                                                                                      >  1:17259/1‑248 (MQ=255)                                                  cgtggtgtgATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGCTCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCCCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGagtagcagc                                                                                                                                                                                                    >  2:366371/5‑248 (MQ=255)                                                      gtgtgATCAACGCCGTCGATAGTCAGCGCACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAAtt                                                                                                                                                                                              >  2:220108/1‑250 (MQ=255)                                                      gtgtgATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAACACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACCGAGTAGCAGCAGAAtt                                                                                                                                                                                              >  2:140136/1‑250 (MQ=255)                                                      gtgtgATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCCGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAAtt                                                                                                                                                                                              >  1:653924/1‑250 (MQ=255)                                                       tgtgATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATAGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTa                                                                                                                                                                                             >  1:469040/1‑250 (MQ=255)                                                        gtgATCAACGACGTCGATAGTCAGCGAACCGATGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCGCGGTAGTGTTTCAATCACGGTAAACAGTAGAAAACAGAGTAGCAGCAGAATTAAc                                                                                                                                                                                           <  1:237695/251‑1 (MQ=255)                                                             caaCGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGt                                                                                                                                                                                      <  1:513750/249‑1 (MQ=255)                                                             caaCGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTGACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGt                                                                                                                                                                                      <  2:83218/249‑1 (MQ=255)                                                              aaCGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTg                                                                                                                                                                                     >  1:702875/2‑251 (MQ=255)                                                               acgccggcGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGa                                                                                                                                                                                    <  1:563089/244‑1 (MQ=255)                                                                        aTAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGt                                                                                                                                                                           >  1:30549/1‑251 (MQ=255)                                                                           gTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCt                                                                                                                                                                        <  1:486554/251‑1 (MQ=255)                                                                               gCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTTCTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCt                                                                                                                                                                    <  2:653924/251‑1 (MQ=255)                                                                                     cGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGAt                                                                                                                                                               >  1:67130/1‑250 (MQ=255)                                                                                        gggctggcGGCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTGCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGTCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCTGTCCTGGCTGCGATAGGa                                                                                                                                                           <  1:156338/250‑1 (MQ=255)                                                                                         ggctggcGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAg                                                                                                                                                          <  1:73713/251‑1 (MQ=255)                                                                                                 tCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTTCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCGCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAAtt                                                                                                                                                    <  1:309548/249‑1 (MQ=255)                                                                                                    gAAAACCGGCTATCAAGTTCAGCAGGGTACTTG‑ACCCGGCGCCGCTTGGCCCGAGGATCGCTACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGTTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTg                                                                                                                                               <  2:26406/251‑1 (MQ=255)                                                                                                        aCCGGCGAGCAAGTTACGCAGGGTACTTTTCCCTG‑CGCCGCTTGGCCCGAGGATCTCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAAGTTGCTGGTag                                                                                                                                           <  2:93357/251‑1 (MQ=255)                                                                                                              gATCAAGTTCAGAAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCAACTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACAGCGCCGTACTGGCTGCGTTAGGAGCCAATTTGCTGGTAGAGATa                                                                                                                                      >  1:219478/1‑250 (MQ=255)                                                                                                                   agttCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGa                                                                                                                                                                >  1:365887/3‑219 (MQ=255)                                                                                                                   agttCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGa                                                                                                                                                                <  2:365887/217‑1 (MQ=255)                                                                                                                    gttCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACACGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGCGATAAAACGGc                                                                                                                               >  1:232299/2‑251 (MQ=255)                                                                                                                     ttCAGCGGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGGGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGc                                                                                                                               >  2:600197/1‑250 (MQ=255)                                                                                                                      tCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGc                                                                                                                               >  1:628535/1‑249 (MQ=255)                                                                                                                         gcagGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGt                                                                                                                          <  1:425865/251‑1 (MQ=255)                                                                                                                         gcagGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGt                                                                                                                          <  1:274904/251‑1 (MQ=255)                                                                                                                            gggTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCg                                                                                                                       <  2:30549/251‑1 (MQ=255)                                                                                                                                  ttttACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGt                                                                                                                                             <  1:344441/223‑1 (MQ=255)                                                                                                                                  ttttACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGGGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGt                                                                                                                                             >  2:344441/1‑223 (MQ=255)                                                                                                                                  ttttACCTG‑CGCCGCTTGGCACGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGAGCGGAAatc                                                                                                                 >  1:564100/1‑251 (MQ=255)                                                                                                                                            cgccgcTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTAcc                                                                                                        >  1:142364/1‑251 (MQ=255)                                                                                                                                              ccgcTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAACACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCg                                                                                                       >  1:429390/1‑250 (MQ=255)                                                                                                                                               cgcTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTACGCTAAAACGCATCGGAAAATGTTGGAAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGAGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATCAAACGGCAGGGTGCGGAAATCATCGTTACCGa                                                                                                      >  1:609750/1‑250 (MQ=255)                                                                                                                                                                          cgccgcGTTCTACCGGTAAGCTAAAACGCAACGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTGTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTGGCTGGTAGAGATAAAACGGCAGGGTGCGGAAAGCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑Ca                                                                          <  2:322304/249‑1 (MQ=255)                                                                                                                                                                           gccgcGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTCGCAGCAGAATTAACGCGGTGACCGCGACGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGAGGAAATCATCGTTACCGAAC‑AACGCCACCCCACCAAAATCAC‑AAATc                                                                       >  1:453709/1‑251 (MQ=255)                                                                                                                                                                               cgTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTAGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAAGCAC‑CAATCGAc                                                                    >  1:413406/1‑250 (MQ=255)                                                                                                                                                                                gTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCACCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGTAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCTGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAg                                                                  <  1:42736/251‑1 (MQ=255)                                                                                                                                                                                gTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCCCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAg                                                                  <  2:447914/251‑1 (MQ=255)                                                                                                                                                                                   ctaCCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTa                                                                                                                                                                                                                                             <  1:415417/77‑1 (MQ=255)                                                                                                                                                                                   ctaCCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTa                                                                                                                                                                                                                                             >  2:415417/3‑79 (MQ=255)                                                                                                                                                                                     aCCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCa                                                                <  1:621366/248‑1 (MQ=255)                                                                                                                                                                                         ttAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCCCCACGCCAAAATCAC‑CAATCCACAGCACGCaggc                                                         >  2:295334/1‑248 (MQ=255)                                                                                                                                                                                         ttAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCaggc                                                         >  1:301157/1‑248 (MQ=255)                                                                                                                                                                                              cTAAAACGCATCGGCCAATGTTGGTAAATCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAACAGAATTAACGCGGTGACCGCGCCGTCCTGGCCGCGAAAGGAGCCAATTTGCTGGTAGAGCTAAAACGGAAGGGTGCGGCAATCATCGTTACCTACCCAAC‑CCACCCCGCCAAAAT‑ACACAAACGACAACACGCA‑GCaaaaa                                                    >  2:356494/1‑248 (MQ=255)                                                                                                                                                                                                      cATCGGCAAATGTTGGTACAGCCAGGTGATATCAGTCCGTTTTAACCTCTCGACCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCCAGCAAAAACCAGCGcc                                            >  2:324342/1‑251 (MQ=255)                                                                                                                                                                                                         cGGCAAATGTTTGTAAATCCAGGTGAAATCAGGCAGTTGTAACATCGCGCCCCGGTAGTTTTTCAAGCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAAGTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACAACGCCAAAATCCC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTgc                                         <  2:653675/251‑1 (MQ=255)                                                                                                                                                                                                         cGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTg                                          <  1:602166/250‑1 (MQ=255)                                                                                                                                                                                                                       aaaGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGttt                           <  1:220108/251‑1 (MQ=255)                                                                                                                                                                                                                          gCCAGGTGATATCAGTCAGTTGTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAg                        <  1:65405/251‑1 (MQ=255)                                                                                                                                                                                                                            cAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAACCAGCGCCTGCGCCAGTGGACGTTTCCggg                      >  2:577812/1‑251 (MQ=255)                                                                                                                                                                                                                                         tcagtTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTcgac         >  1:560012/1‑248 (MQ=255)                                                                                                                                                                                                                                                tAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGCCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGCAATCATCGTTACCGAAC‑CACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTACAGGGCGCGCAGCTCGACCACtttt  >  2:251620/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       GGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGTGTGATCTACGCCATCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTT  >  NC_000913/72684‑73172

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