Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	74,015	C→T	Q169Q (CAG→CAA)	thiP ←	thiamine ABC transporter membrane subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	74,015	0	C	T	95.5%	129.0 / ‑1.5	44	Q169Q (CAG→CAA)	thiP	thiamine ABC transporter membrane subunit
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (16/26);  total (17/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.75e-01

CATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGG‑GCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGC  >  NC_000913/73769‑74255                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  cATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGAGCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCAGTCGACGAAGCGTAAAAAATGCTAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACgttgtt                                                                                                                                                                                                                                               <  1:15866/251‑1 (MQ=255)                    atcgCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGGTTTGTCGCCGTAACCACGGCCATTCGAAGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGGCGTTGGTCGCCGGGGATGTTTTcc                                                                                                                                                                                                                              >  2:47971/4‑249 (MQ=255)                        ccGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAATGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGcc                                                                                                                                                                                                                        <  2:186676/251‑1 (MQ=255)                         cgcgcgGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCt                                                                                                                                                                                                                       <  1:117618/251‑1 (MQ=255)                                gCAGGATCGTAGTCGTAACTCAGCGCATGATAGATTGCCAGCTCGATAGTGGTTGCCTGAGGACCGCCACCCAGCGATCCCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGCCGGGATTTGTCGCCGTAACCACGGCCATTAGCAGAAGCGGAAAAAATACCAGCCGCGCATACCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGGCTGTTTTCCACTGCCTGGAgac                                                                                                                                                                                                                 >  1:167935/1‑248 (MQ=255)                                             cGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGcc                                                                                                                                                                                                   <  1:36085/251‑1 (MQ=255)                                               tAACTCAGCGCCTGATAGCTTGCCAGCTCGATAGCGGTTGCCTGCGGACCGCCGACCAGCGATAGCACGGTGGCGAAGCTGGCGCAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCCCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCCGCCGCGCACCCCAAGCTGGGCCGCAACCTGAAGTTGTTCGCCCGGACTGTTTTCCAGCGGCTTGAGTAATAAGCGGCCCGcccc                                                                                                                                                                                                 >  2:22548/1‑249 (MQ=255)                                               tAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGACGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCCTCCCAAGCTGGGCGGCAAGCTGACGTTGTTCCCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCCCGCCac                                                                                                                                                                                                 >  2:134292/1‑250 (MQ=255)                                                aaCTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATc                                                                                                                                                                                                <  2:105725/251‑1 (MQ=255)                                                        gcCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGAtt                                                                                                                                                                                        >  1:101873/1‑251 (MQ=255)                                                           tGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCACCGGTAGATTaa                                                                                                                                                                                      >  2:140141/1‑250 (MQ=255)                                                            gATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTaaaa                                                                                                                                                                                    <  1:265232/251‑1 (MQ=255)                                                               aGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTaaaaaa                                                                                                                                                                                  <  1:211727/250‑1 (MQ=255)                                                                           cgagaGTGGTAGCCTGCGGACCGCCGCCCAGCGATAGCTCGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGGGCAGCAACCAGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCGGACGTTGATCGCCGGGAATGTTTTCCAGGGGCTGGAGTAATAAGCGGCTCGCGATTGGTTGATTaaaaaagtcgtgcgccag                                                                                                                                                                      <  2:106699/246‑13 (MQ=255)                                                                                 tGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGC‑GCCAGCAAAAt                                                                                                                                                               >  2:66468/1‑250 (MQ=255)                                                                                  ggttgCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGC‑GCCAGCAAAATa                                                                                                                                                              <  1:199553/246‑1 (MQ=255)                                                                                          gCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTTGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGC‑GCCAGCAAAATACCTTGCAg                                                                                                                                                      <  1:88452/250‑1 (MQ=255)                                                                                              aCCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCAt                                                                                                                                                                                                 <  2:143631/204‑1 (MQ=255)                                                                                              aCCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCAt                                                                                                                                                                                                 >  1:143631/1‑204 (MQ=255)                                                                                                     cccAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGC‑GCCAGCAAAATACCTTGCAGGCCGTAGGGCg                                                                                                                                           >  1:200403/1‑250 (MQ=255)                                                                                                        aGCGATAGAACGGTGGCACAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCAATCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCAGGGCGGCAAGCTGACGTTGTTCCCCGGGAATGGTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGGAGATTAAAAAATA‑AGGGCGCCAGCAAAATACCTTTCAGGACGTCGGGCaaaa                                                                                                                                        >  2:260341/1‑250 (MQ=255)                                                                                                               tcACGGTGGAGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGGTCGCCGGGAATGTTTTCCAGTGCCTGGCGGAATAAGGGGCCCGCCATCGGTGGATTAAAAAATACGTGC‑GCACGCAAAATACCTTGCCGGCCGTAGGGGGAAAAgggcccc                                                                                                                                >  2:91350/2‑246 (MQ=255)                                                                                                                cACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGC‑GCCATCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACt                                                                                                                               <  1:30140/251‑1 (MQ=255)                                                                                                                   ggTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGC‑GCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCa                                                                                                                            <  1:115670/251‑1 (MQ=255)                                                                                                                    gTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGC‑GCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAg                                                                                                                           <  1:69645/251‑1 (MQ=255)                                                                                                                            gCTGGGGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAAGAAGCGGCTCGCCAGCGGTAGATTAAAAAATACGTGC‑GCCAGCAAAATACCTTGCATGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCg                                                                                                                   <  2:239166/251‑1 (MQ=255)                                                                                                                                     aCAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGACGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGC‑GCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGc                                                                                                             >  1:113565/1‑248 (MQ=255)                                                                                                                                           cATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGC‑GCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTg                                                                                                    <  2:74792/251‑1 (MQ=255)                                                                                                                                             aaaaGAGAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTTGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGC‑GCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTgcc                                                                                                  <  2:192086/250‑1 (MQ=255)                                                                                                                                                 gATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGC‑GCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTgccagcc                                                                                              <  1:37524/251‑1 (MQ=255)                                                                                                                                                        gcagcaACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGC‑GCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTg                                                                                       <  1:209433/251‑1 (MQ=255)                                                                                                                                                         cagcaACCGGCGGGATTTGTCGCCGTAACCACGGCCAGTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGC‑GCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTgg                                                                                      <  2:218017/251‑1 (MQ=255)                                                                                                                                                                          tGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGC‑GCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTaa                                                                      <  2:104693/250‑1 (MQ=255)                                                                                                                                                                             cGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCCGCCGCGCATCCCAAGCTGGGCGGCAAGATGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGC‑GCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAAt                                                                  >  1:135052/1‑251 (MQ=255)                                                                                                                                                                                            cATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGC‑GCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGcaac                                                   <  2:38912/251‑4 (MQ=255)                                                                                                                                                                                              ttCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGC‑GCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCCGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAAcc                                                  >  2:233258/1‑250 (MQ=255)                                                                                                                                                                                                  acgaAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGC‑GCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAgcac                                             <  1:27196/251‑1 (MQ=255)                                                                                                                                                                                                   cgaAGCGGACAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGC‑GCCAGCAAAATACCTTGCAGGCCGTAGGGTGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCAcc                                            >  1:195042/1‑251 (MQ=255)                                                                                                                                                                                                                  gccagccgcGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGC‑GCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTc                             <  1:47971/243‑1 (MQ=255)                                                                                                                                                                                                                      gccgcGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGC‑GCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTg                         <  2:34322/247‑1 (MQ=255)                                                                                                                                                                                                                                ccAAGGTGGGCGGCAAGCTGGCGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGC‑GCCAGCAAAATACCTGGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGAGTAAGGGCCTTGCGCACAGACg               <  2:172917/251‑1 (MQ=255)                                                                                                                                                                                                                                ccAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGC‑GCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACg               >  2:59373/1‑251 (MQ=255)                                                                                                                                                                                                                                          cggcAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGG‑GCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGcc     >  2:235561/1‑251 (MQ=255)                                                                                                                                                                                                                                          cggcAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGC‑GCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGcc     <  1:140141/251‑1 (MQ=255)                                                                                                                                                                                                                                             cAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGG‑GCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGc  <  1:235561/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  CATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGG‑GCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGC  >  NC_000913/73769‑74255

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