Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,442,824	C→T	S148F (TCC→TTC)	tamA →	translocation and assembly module subunit TamA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,442,824	0	C	T	98.2%	176.0 / ‑6.0	54	S148F (TCC→TTC)	tamA	translocation and assembly module subunit TamA
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (18/35);  total (18/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GCGAAGGTCTGAAAGCGCTGGGTTATTACCAGCCGACCATTGAATTTGATCTCCGTCCACCGCCAAAGAAAGGGCGGCAGGTATTGATCGCCAAAGTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGACCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAAGAGGGCGATGAGTACGAATCGAAAGATCTGGCAGAACTGA  >  NC_000913/4442578‑4443060                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                             gcgAAGGCCTGAAAGCGCTGGGTTATTACCAGCCGACCATTGAATTTGATCTCCGTCCACCGCCAAAGAAAGGGCGGCAGGTATTGATCGCCAAAGTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGACCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTTCtt                                                                                                                                                                                                                                           >  1:453596/1‑250 (MQ=255)    aaGGTCTGAAAGCGCTGGGTTATTACCAGCCGACCATTGAATTTGATCTCCGTCCACCGCCAAAGAAAGGGCGGCAGGTATTGATCGCCAAAGTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGACCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTTCTTAAcc                                                                                                                                                                                                                                       >  2:46907/1‑251 (MQ=255)    aaGGTCTGAAAGCGCTGGGTTATTACCAGCCGACCATTGAATTTGATCTCCGTCCACCGCCAAAGAAAGGGCGGCAGGTATTGATCGCCAAAGTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGACCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTTCTTAAcc                                                                                                                                                                                                                                       <  2:462038/251‑1 (MQ=255)       gTCTGAAAGCGCTGGGTTATTACCAGCCGACCATTGAATTTGATCTCCGTCCACCGCCAAAGAAAGGGCGGCAGGTATTGATCGCCAAAGTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGACCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTTCTTAACCAGc                                                                                                                                                                                                                                    >  1:236842/1‑251 (MQ=255)        tCTGAAAGCGCTGGGTTATTACCAGCCGACCATTGAATTTGATCTCCGTCCACCGCCAAAGAAAGGGCGGCAGGTATTGATCGCCAAAGTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGACCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTTCTTAACCAGc                                                                                                                                                                                                                                    >  2:132130/1‑250 (MQ=255)             aaGCTCTGGGTTATTACCAGCCGACCATTGAATTTGATCTCCGTCCACCGCCAAAGAAAGGGCGGCAGGTATTGATCGCCAAAGTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGACCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTTCTTAACCAGCAttgcg                                                                                                                                                                                                                              <  1:286702/251‑1 (MQ=255)                  cTGGGTTATTACCAGCCGACCATTGAATTTGATCTCCGTCCACCGCCAAAGAAAGGGCGGCAGGTATTGATCGCCAAAGTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGACCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTTCTTAACCAGCAttgcgttgcg                                                                                                                                                                                                                         >  1:73824/1‑251 (MQ=255)                         attaCCAGCCGACCATTGAATTTGATCTCCGTCCACCGCCAAAGAAAGGGCGGCAGGTATTGATCGCCAAAGTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGACCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTTCTTAACCAGCATTGCGTTGCGTAAAGGt                                                                                                                                                                                                                  >  1:7128/1‑251 (MQ=255)                         attaCCAGCCGACCATTGAATTTGATCGCCGTCCACCGCCAAAGAAAGGGCGGCAGGTATTGATCGCCAAAGTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGACCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTTCTTAACCAGCATTGCGTTGCGTAAAgg                                                                                                                                                                                                                   <  1:637113/250‑1 (MQ=255)                               aGCCGACCATTGACTTTGATCTACGTCCACCGCCACAGAAAGTGCGCCAGGGATTGACGGCCATAGTCACGCCAGGCGCGCCGGTGTTACTTGGCGGCACCGATGTGGTATCGCGCGGCGGCGCGCGGCACGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGGACGGTACTGAACCAAGGCGATTATGAAAATTTCAAAAAGTTCTTAACCAGCATTGCGGTGCGTAAAGGTTATTTc                                                                                                                                                                                                            >  2:419517/1‑251 (MQ=255)                                 ccGACCATTGAATTTGATCTCCGTCCACCGCCAAAGAAAGGGCGGCAGGTATTGATCGCCAAAGTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGACCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTTCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCg                                                                                                                                                                                                           >  1:642155/1‑250 (MQ=255)                                   gACCATTGAATTTGATCTCCGTCCACCGCCAAAGAAAGGGCGGCCGGAATTGATCGCCACAGTCACGCCAGGCGTGCCGGCGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGAGCGCGGACCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGAACTGAACCAGGACGATTATGAAAATTTCAAAAAGTTCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGAt                                                                                                                                                                                                         >  2:643497/1‑250 (MQ=255)                                           aaTTTGATCTCCGTCCACCGCCAAAGAAAGGGCGGCAGGTATTGATCGCCAAAGTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGACCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTTCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAAttt                                                                                                                                                                                                <  2:426220/251‑1 (MQ=255)                                                     ccGTCCACCGCCAAAGAAAGGGCGGCAGGTATTGGTCGCCAAAGTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGACCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTTCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAgc                                                                                                                                                                                        <  1:560094/249‑1 (MQ=255)                                                     ccGTCCACCGCCAAAGAAAGGGCGGCAGGTATTGATCGCCAAAGTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGACCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTTCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCaaa                                                                                                                                                                                          <  2:337170/247‑1 (MQ=255)                                                                cAANGAAAGGGCGGCAGGTATTGATCGCCAAAGTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGACCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTTGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTTCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCTGCTGGGCAtt                                                                                                                                                                           <  2:603365/251‑1 (MQ=255)                                                                                ggTATTGATCGCCAAAGTCACGCCAGGCGTGCCGGTGCTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGACCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTTCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCAt                                                                                                                                                            >  1:284874/1‑250 (MQ=255)                                                                                      gATCGCCAAAGTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGACCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTTCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCt                                                                                                                                                     >  1:562766/1‑251 (MQ=255)                                                                                       aTCGCCAAAGTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGACCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTTCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGGCTGCATAAAGCCtt                                                                                                                                                    <  1:445870/251‑1 (MQ=255)                                                                                                       cAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGACCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTTCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATACAGCCTTCTGGGATATTGATtat                                                                                                                                    >  2:239242/1‑251 (MQ=255)                                                                                                              gCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGACCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTTCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTg                                                                                                                             <  1:225623/251‑1 (MQ=255)                                                                                                                 ggTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGACCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTTCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCg                                                                                                                          <  1:90866/251‑1 (MQ=255)                                                                                                                       aagtGGCGGCACCTATGTGGTATTGCGCGGCGGCGCGCGGACCGATAAAGACCATTTGAAATTGCTCGATACTCTCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTTCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCTGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATGGATTAGAACAGTGGCTAActtt                                                                                                                    <  2:252711/248‑4 (MQ=255)                                                                                                                       aaTTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGACCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTTCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGtt                                                                                                                    <  1:74683/251‑1 (MQ=255)                                                                                                                        aTTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGACCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTTCTTAACCAGCATTGCGTTGCGTAAAGGTTATGTCGATAGCGAACTTACCAAATCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTTTGGGATATTGATTATAACATTGGCGAACGTTa                                                                                                                   <  1:232312/251‑1 (MQ=255)                                                                                                                               gcACCGATGTGGTATTGCGCGGCGGCGCGCGGACCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTGCAAAAAGTTCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCttt                                                                                                            <  1:588987/251‑1 (MQ=255)                                                                                                                                  ccGATGTGGTATTGCGCGGCGGCGCGCGGACCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTTCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTggg                                                                                                         <  1:152874/251‑1 (MQ=255)                                                                                                                                          gTATTGCGCGGCGGAGCGCGGAACGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTTCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGAc                                                                                                 <  2:359761/251‑1 (MQ=255)                                                                                                                                              tGCGCGGCGGCGCGCGGACCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTTCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCtt                                                                                              <  1:56728/250‑1 (MQ=255)                                                                                                                                                     cggccCGCGTACCGATAAAGACTATTTGAAATTGCTCGATACGCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTTCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTGCTGGGATATTGATTATAACAGTGGCGACCGTTACCTCTTTGGGGATGTGACCTTTCAAGGAt                                                                                      <  2:81797/251‑1 (MQ=255)                                                                                                                                                         gcgcgGACCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTTCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATcaca                                                                                  <  2:87914/251‑1 (MQ=255)                                                                                                                                                         gcgcgGACCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGAACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTTCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATcac                                                                                   >  1:409938/1‑250 (MQ=255)                                                                                                                                                           gcgGACCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTTCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACaaa                                                                                <  1:186021/251‑1 (MQ=255)                                                                                                                                                             ggACCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTTCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTa                                                                                                                                                                                               >  2:187473/1‑138 (MQ=255)                                                                                                                                                             ggACCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTTCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTa                                                                                                                                                                                               <  1:187473/138‑1 (MQ=255)                                                                                                                                                                  gATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTTCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGa                                                                         <  1:132130/251‑1 (MQ=255)                                                                                                                                                                     aaaGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTTCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATg                                                                       >  2:160594/1‑250 (MQ=255)                                                                                                                                                                      aaGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTTCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGaa                                                                     <  1:367302/251‑1 (MQ=255)                                                                                                                                                                                   aTTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTTCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATc                                                        <  2:7128/251‑1 (MQ=255)                                                                                                                                                                                             aCTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTTCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGtt                                              <  1:448873/251‑1 (MQ=255)                                                                                                                                                                                             aCTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTTCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAAGCTGGTGCCGt                                               <  2:73824/250‑1 (MQ=255)                                                                                                                                                                                             aCTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGt                                               <  1:259118/250‑1 (MQ=255)                                                                                                                                                                                                cGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTTCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGttt                                             <  1:59383/249‑1 (MQ=255)                                                                                                                                                                                                cGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTTCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTaa                                           <  1:30324/251‑1 (MQ=255)                                                                                                                                                                                                cGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTTCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTa                                            <  2:236804/250‑1 (MQ=255)                                                                                                                                                                                                     ggCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTTCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAAGAgg                                      >  1:212557/1‑251 (MQ=255)                                                                                                                                                                                                           tGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTTCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAAGAGGGCGAt                                 <  1:219440/250‑1 (MQ=255)                                                                                                                                                                                                                 ggTACTGAACCAGGTCGATTATGAAAATTTCAAAAAGTTCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTGAAAGATCGCGATGAGTACg                          <  2:453596/251‑1 (MQ=255)                                                                                                                                                                                                                 ggTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTTCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAAGAGGGCGATGAGTACg                          <  2:562766/251‑1 (MQ=255)                                                                                                                                                                                                                       gAACCAGGGCGATTATGAAAATTTCAAAAAGTTCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAAGAGGGCGATGAGTACGAATc                      <  1:28357/249‑1 (MQ=255)                                                                                                                                                                                                                        aaCCAGGGCGATTATGAAAATTTCAAAAAGTTCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAAGAGGGCGATGAGTACGAATCGaa                   <  2:212557/251‑1 (MQ=255)                                                                                                                                                                                                                                 gATTATGAAAATTTCAAAAAGTTCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAAGAGGGCGATGAGTACGAATCGAAAGATCTGGc          >  2:282082/1‑251 (MQ=255)                                                                                                                                                                                                                                 gATTATGAAAATTTCAAAAAGTTCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAAGAGGGCGATGAGTACGAATCGAAAGATCTGGc          >  1:516002/1‑251 (MQ=255)                                                                                                                                                                                                                                         aaaTTTCAAAAAGTTCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTGTAAAGAGGGCGATGAGTACGAATCGAAAGATCTGGCAGAACTGa  <  1:46907/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                             GCGAAGGTCTGAAAGCGCTGGGTTATTACCAGCCGACCATTGAATTTGATCTCCGTCCACCGCCAAAGAAAGGGCGGCAGGTATTGATCGCCAAAGTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGACCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCCTTAACCAGCATTGCGTTGCGTAAAGGTTATTTCGATAGCGAATTTACCAAAGCGCAGCTGGGCATTGCGCTCGGCCTGCATAAAGCCTTCTGGGATATTGATTATAACAGTGGCGAACGTTACCGCTTTGGGCATGTGACCTTTGAAGGATCACAAATCCGCGATGAATACCTGCAAAATCTGGTGCCGTTTAAAGAGGGCGATGAGTACGAATCGAAAGATCTGGCAGAACTGA  >  NC_000913/4442578‑4443060

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