Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	12,879	T→C	Y239Y (TAT→TAC)	dnaK →	chaperone protein DnaK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	12,879	0	T	C	91.1%	144.6 / 4.9	56	Y239Y (TAT→TAC)	dnaK	chaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base T (3/2);  new base C (21/30);  total (24/32)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.13e-01

GGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGCACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGGGGTGAAGACTTCGACAGCCGTCTGATCAACTATCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCATACATCACTGCAGACGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAG  >  NC_000913/12634‑13116                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              ggCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTGCTGGCAGCCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGt                                                                                                                                                                                                                                          <  1:235995/251‑1 (MQ=255)          gCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAAtt                                                                                                                                                                                                                                 <  1:436727/251‑1 (MQ=255)           cTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGATTTc                                                                                                                                                                                                                                >  2:444713/1‑251 (MQ=255)             ggTCTGGAATTAAAACGGTTCATCAACGAACCGACCACAGCTGCGCTGGCTTACGTTCTGGACAAAGGTACTGGCAACCGTACTAGCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTGGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTATGATCAACGACCTGGTTGAAGAATTCaa                                                                                                                                                                                                                              <  2:466947/251‑1 (MQ=255)              gTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGAGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTGGAAGAATTCacg                                                                                                                                                                                                                             >  2:410941/1‑249 (MQ=255)                ttGGAAGTAAAACGTATCATCAACACACCGACCGCAGCTGCGCTGGCTTACGTTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACATTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATTCCCACCTGGGTGGTGAAGACTGCGACAGCCGTCTGATCAACTGCCTGTTTGAAGAATTCaaga                                                                                                                                                                                                                            <  2:37624/249‑1 (MQ=255)                     aGTAAAACGTATCATCAACGAACCGACCGCAGATGCGATGGCGTACGGGCGGGAAAAAGGGACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCTCCTGGGTGGTGAAGACTGCGACAGACGTCTGATCAACTACCTGGTTGAAGAATTCaagaaagagc                                                                                                                                                                                                                      <  2:118273/251‑3 (MQ=255)                      gTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCa                                                                                                                                                                                                                     <  1:79492/251‑1 (MQ=255)                          aacttATCATCAACGAACCGACCGCAGCTGCGCTTGCTTACTGTCTGGACAAATGTTCTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCCGGGc                                                                                                                                                                                                                 <  2:359977/247‑1 (MQ=255)                           aCGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGc                                                                                                                                                                                                                 >  1:156207/1‑250 (MQ=255)                           aCGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGc                                                                                                                                                                                                                 >  1:452357/1‑250 (MQ=255)                            cGTATCAACAACGAACCGTCCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGTCCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCAt                                                                                                                                                                                                               <  1:183988/251‑1 (MQ=255)                                      aCGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTgcgc                                                                                                                                                                                                     >  1:294850/1‑251 (MQ=255)                                        gAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGAGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACAAACGGTGATACACACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCa                                                                                                                                                                                                    >  1:50896/1‑250 (MQ=255)                                           ccgaccgCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGa                                                                                                                                                                                                <  1:464188/251‑1 (MQ=255)                                                                tACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCt                                                                                                                                                                           >  1:499822/1‑251 (MQ=255)                                                                tACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCA‑GAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGcc                                                                                                                                                                            <  1:48052/249‑1 (MQ=255)                                                                        ggACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAaag                                                                                                                                                                      <  1:117177/248‑1 (MQ=255)                                                                           cAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCgg                                                                                                                                                                >  2:97799/1‑251 (MQ=255)                                                                           cAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATAAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCgg                                                                                                                                                                >  2:322642/1‑251 (MQ=255)                                                                                        aaCCGTATTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGTGCCTGAAAGAAGCGGCAGAAAAAGCGaa                                                                                                                                                   <  1:159682/251‑1 (MQ=255)                                                                                          ccgtgttATCGCGGTTTATGACCTGGGTGGTGGTTCTTTCGATATTGCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGGTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGaaaa                                                                                                                                                 <  2:506924/245‑1 (MQ=255)                                                                                          ccGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGaaaa                                                                                                                                                 <  2:499822/251‑1 (MQ=255)                                                                                                 aTCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAAGCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACt                                                                                                                                          <  1:443461/251‑1 (MQ=255)                                                                                                 aTCGCGGTTTATGACCGTGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCATCCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTCCCTGGTTGAAGAAGTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCGGAAAGAAGCGGCAGAAAAAGCGAAAATCGaa                                                                                                                                            <  1:427324/249‑1 (MQ=255)                                                                                                   cgcgGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTg                                                                                                                                         <  1:444713/250‑1 (MQ=255)                                                                                                   cgcgGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGt                                                                                                                                        <  1:322642/251‑1 (MQ=255)                                                                                                       gTTTAGGACCTGGGTGGTGGTACTTTCGATATTTCGAGTATCGAAATCGACGAAGTTGACGGCGAAAAAATCGTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGGGGTGAAGACTTCGACAGCCGTCTGATCAACTATCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTTAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTc                                                                                                                                    <  2:185339/251‑1 (MQ=255)                                                                                                         ttATGACCTGGGTGGTGGTACTTTCGATCTTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGGGGTGAAGACTTAGACAGCCGTCTGATCAACTCCCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACATTCCGCTGGCAATCAAGAGCCTGAAAGACGCGGCAGAAAACCCGCAAATCGAACAGTCttacg                                                                                                                                  >  2:485647/1‑248 (MQ=255)                                                                                                             gACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGACAAAGCGAAAATCGAACTGTCTTCCGCTc                                                                                                                               >  1:26650/1‑250 (MQ=255)                                                                                                               ccTGGGTGGTTGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCGTCGAAGTTCTGGCAACCATCGGTGATACCCACCTGTGTGGTGAAGACGTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTcatc                                                                                                                            <  1:306927/251‑3 (MQ=255)                                                                                                               ccTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTcagc                                                                                                                            <  1:306924/251‑1 (MQ=255)                                                                                                                 tGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTcagc                                                                                                                            >  2:207911/1‑249 (MQ=255)                                                                                                                   ggtggtggtACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGCAAATCGAACTGTCTTCCGCTCAGCAGAc                                                                                                                        >  1:31468/1‑251 (MQ=255)                                                                                                                    gtggtggtACTTTCGAGATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTGCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACGTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGAcc                                                                                                                       <  2:452357/251‑1 (MQ=255)                                                                                                                          gtACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGtt                                                                                                                 >  1:145208/1‑251 (MQ=255)                                                                                                                                       ttttattatCGAAATCGACGAAGTTGACGGCGAAAAAACCCTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGGCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGCCCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTAc                                                                                                     <  2:251452/247‑1 (MQ=255)                                                                                                                                              atCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGGGGTGAAGACTTCGACAGCCGTCTGATCAACTATCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCATACATCACTGc                                                                                             >  1:67081/1‑251 (MQ=255)                                                                                                                                                     ccgacgaAGTTGACGGCGAAAAAACCTTCGAAGTTCGGGCAACCAACGGTGATACCCACCTGGGTGGGGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGAtgcg                                                                                      <  2:145208/250‑1 (MQ=255)                                                                                                                                                       gacgaAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGAc                                                                                    >  1:281268/1‑251 (MQ=255)                                                                                                                                                        acgaAGTGGACGGCGAAAAAACCTTCGAAGTTCTGGCACCCAACGGTGATACCCACCTGGGTGGTGAAGACATCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGAtgcgccc                                                                                   <  2:157497/251‑4 (MQ=255)                                                                                                                                                             gTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATCCCCACCGGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTcc                                                                              <  2:156207/251‑1 (MQ=255)                                                                                                                                                                       aaaaaaCCGTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATg                                                                    <  1:271454/251‑1 (MQ=255)                                                                                                                                                                             ccTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCACTGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAACTCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACACGAACAt                                                               >  2:295432/1‑250 (MQ=255)                                                                                                                                                                              cTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATc                                                              >  1:104077/1‑250 (MQ=255)                                                                                                                                                                                          ggCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCGGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTc                                                    <  2:442567/248‑1 (MQ=255)                                                                                                                                                                                                  acGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTg                                         >  1:383068/1‑251 (MQ=255)                                                                                                                                                                                                    ggTGATACCCACCTGGGGGGTGAAGACTTCGACAGCCGTCTGATCAACTATCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGATCAGCAGACCGACGTTAACCTGCCATACATCACTGCAGACGCGACCGGTCCGAAACACATGAACCTCAAAGTGACTCGTGCGAAACTgg                                        >  2:215729/1‑250 (MQ=255)                                                                                                                                                                                                     gTGATACCCACCTGGGGGGTGAAGACTTCGACAGCCGTCTGATCAACTATTTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCATACAt                                                                                                   <  1:383435/190‑1 (MQ=255)                                                                                                                                                                                                     gTGATACCCACCTGGGGGGTGAAGACTTCGACAGCCGTCTGATCAACTATTTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCATACAt                                                                                                   >  2:383435/1‑190 (MQ=255)                                                                                                                                                                                                       gATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACAGGAACATCAAAGTGACTCGTGCGAAACTGGaa                                      >  2:505951/1‑249 (MQ=255)                                                                                                                                                                                                                       gTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGGGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTg                    >  2:85396/1‑251 (MQ=255)                                                                                                                                                                                                                           aGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTaa                <  2:104077/251‑1 (MQ=255)                                                                                                                                                                                                                                 cGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGtt          <  1:26321/251‑1 (MQ=255)                                                                                                                                                                                                                                       ccGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTg    >  1:181632/1‑251 (MQ=255)                                                                                                                                                                                                                                         gTCTGATCAACTACATGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAg  <  1:321609/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              GGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGCACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGGGGTGAAGACTTCGACAGCCGTCTGATCAACTATCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCATACATCACTGCAGACGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAG  >  NC_000913/12634‑13116

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