Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	74,069	C→T	L151L (CTG→CTA)	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,069	0	C	T	94.7%	109.5 / ‑1.5	40	L151L (CTG→CTA)	thiP	thiamine ABC transporter membrane subunit
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (14/22);  total (17/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.88e-01

GCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGG‑GCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGAC  >  NC_000913/73824‑74311                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    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)                                              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)                                                                                                                                                                                                        cGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCAATCGGAAGATTAAAAAATACGTGC‑GCCAGCAAAATACCTTGAAGGCCGTAGGGCGAAAAGGGCCACTACAGAACGAGCGATAGGAAGAGAGTTGCAAGCCAGCCCTGGCGACCATAGAAGCTAAGAATGCCGAAAAAAGCAACCAGCAACGGGAGGATTAAGGGCATTGCGCACAGACGCAACAGCGCGAGACGACACGGAAAgagac                                          >  1:52714/1‑245 (MQ=255)                                                                                                                                                                                                            aaTGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGC‑GCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGAt                                    >  1:45636/2‑251 (MQ=255)                                                                                                                                                                                                                   tCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGC‑GCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAgcgc                             >  1:160635/1‑251 (MQ=255)                                                                                                                                                                                                                                 taataaGCGGCTCGCCATCGGTAGATTAAAAAATACGTGC‑GCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAat                >  1:108403/1‑250 (MQ=255)                                                                                                                                                                                                                                   ataaGCGGCTCGCCATCGGTAGATTAAAAAATACGTGC‑GCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATc             <  1:66468/251‑1 (MQ=255)                                                                                                                                                                                                                                                 aTCGGTAGATTAAAAAATACGTGC‑GCAAGCAAAAAACCTTGCAGGCCGCAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGAGGTCCCAGCCAGCCCTGGCGACCAGAAACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGTAGGATTAGGGTCATTGCGCACAGACGCACCAGCCCCGGCCGACCCGGAAAGCGCCCGCGATCGAGCGCGCCGGCGCGGAATATCGCggggacgac  >  2:219409/1‑242 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    GCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGG‑GCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGAC  >  NC_000913/73824‑74311

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