Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	13,029	C→T	D289D (GAC→GAT)	dnaK →	chaperone protein DnaK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	13,029	0	C	T	95.8%	139.1 / ‑2.5	50	D289D (GAC→GAT)	dnaK	chaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (23/23);  total (24/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.70e-01

TCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGGGGTGAAGACTTCGACAGCCGTCTGATCAACTATCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTC‑TTCCGCTCAGCAGACCGACGTTAACCTGCCATACATCACTGCAGACGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGT  >  NC_000913/12782‑13274                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                      ccgacgaAGTTGACGGCGAAAAAACCTTCGAAGTTCGGGCAACCAACGGTGATACCCACCTGGGTGGGGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTC‑TTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGAtgcg                                                                                                                                                                                                                                                    <  2:145208/250‑1 (MQ=255)   gacgaAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTC‑TTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGAc                                                                                                                                                                                                                                                  >  1:281268/1‑251 (MQ=255)    acgaAGTGGACGGCGAAAAAACCTTCGAAGTTCTGGCACCCAACGGTGATACCCACCTGGGTGGTGAAGACATCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTC‑TTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGAtgcgccc                                                                                                                                                                                                                                                 <  2:157497/251‑4 (MQ=255)         gTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATCCCCACCGGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTC‑TTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTcc                                                                                                                                                                                                                                            <  2:156207/251‑1 (MQ=255)                   aaaaaaCCGTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTC‑TTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATg                                                                                                                                                                                                                                  <  1:271454/251‑1 (MQ=255)                         ccTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCACTGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAACTCGAACTGTC‑TTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACACGAACAt                                                                                                                                                                                                                             >  2:295432/1‑250 (MQ=255)                          cTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTC‑TTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATc                                                                                                                                                                                                                            >  1:104077/1‑250 (MQ=255)                                      ggCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCGGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTC‑TTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTc                                                                                                                                                                                                                  <  2:442567/248‑1 (MQ=255)                                              acGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTC‑TTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTg                                                                                                                                                                                                       >  1:383068/1‑251 (MQ=255)                                                ggTGATACCCACCTGGGGGGTGAAGACTTCGACAGCCGTCTGATCAACTATCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTC‑TTCCGATCAGCAGACCGACGTTAACCTGCCATACATCACTGCAGACGCGACCGGTCCGAAACACATGAACCTCAAAGTGACTCGTGCGAAACTgg                                                                                                                                                                                                      >  2:215729/1‑250 (MQ=255)                                                   gATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTC‑TTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACAGGAACATCAAAGTGACTCGTGCGAAACTGGaa                                                                                                                                                                                                    >  2:505951/1‑249 (MQ=255)                                                                   gTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTC‑TTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGGGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTg                                                                                                                                                                                  >  2:85396/1‑251 (MQ=255)                                                                       aGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTC‑TTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTaa                                                                                                                                                                              <  2:104077/251‑1 (MQ=255)                                                                             cGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTC‑TTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGtt                                                                                                                                                                        <  1:26321/251‑1 (MQ=255)                                                                                   ccGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTC‑TTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTg                                                                                                                                                                  >  1:181632/1‑251 (MQ=255)                                                                                     gTCTGATCAACTACATGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTC‑TTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAg                                                                                                                                                                <  1:321609/251‑1 (MQ=255)                                                                                                        tGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTC‑TTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGc                                                                                                                                             >  1:327497/1‑251 (MQ=255)                                                                                                           agaaTTCAAGAAAGATCAGGGCATTGACCTGCGCAACGACCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTC‑TCCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAgg                                                                                                                                          <  2:251861/251‑1 (MQ=255)                                                                                                                 aaagaaagaccgggtaattgccatGCGCAACGATCCGCGGGCAATGCAGCGCCGTAAAGAAGCGGCAGAAAAAGCGAAAATAGAAC‑GTCCTTCCGCTCAGCAGACCGACGTTAACCAGCCGTACATCACTGCAGATGCGACCGGTCCGAATCACATGAACATCAAATTGACTCGTGCGAAACTGGAAATCCTGGTTGAAGATCGGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCTCTGCAGGACGcc                                                                                                                                     <  2:37806/227‑2 (MQ=255)                                                                                                                                cATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTC‑TTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAGGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTg                                                                                                                     >  1:67704/1‑251 (MQ=255)                                                                                                                                    gACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTC‑TTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGata                                                                                                                  <  1:207911/250‑1 (MQ=255)                                                                                                                                    gACCTGCGCAACGATCCGCTAGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTC‑TTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGGTCTGGTAAACCGTTCCATTGAGCCGCTGGAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGatat                                                                                                                 >  2:259451/1‑251 (MQ=255)                                                                                                                                                 aTCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTC‑TTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGtt                                                                                                                                                                                           >  1:465045/1‑164 (MQ=255)                                                                                                                                                 aTCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTC‑TTCCGCGCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACAGCAAAGTGACGCGTGCGAAACTGGAAAGCCTGGtt                                                                                                                                                                                           <  2:465045/164‑1 (MQ=255)                                                                                                                                                     gCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTC‑TTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCCCTGCAGGACGCTGGCCTGTCCGTATCTGCTATCGACGACGTTAGCCTCg                                                                                                >  2:358614/1‑251 (MQ=255)                                                                                                                                                       tGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTC‑TTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCAGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGtt                                                                                              >  1:367663/1‑251 (MQ=255)                                                                                                                                                               cAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAAACGAACTGGC‑TCCCGCTCAGCAGCCCGACGTTAACCTGCCGGACATCACTGCAGAGGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGATAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCa                                                                                      <  2:31468/251‑1 (MQ=255)                                                                                                                                                                        aaaGAAGCGGCAGAAAAAGCGAAAATCGAACTGTC‑TTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTAt                                                                             <  1:505454/251‑1 (MQ=255)                                                                                                                                                                          aaaagCGGCAGAAAAAGAGAAAATCGAACTGTC‑TTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGc                                                                           <  2:294850/249‑1 (MQ=255)                                                                                                                                                                               cGGCAGAAAAAGCGAAAATCGAACTGTC‑TTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATg                                                                      >  1:333524/1‑251 (MQ=255)                                                                                                                                                                               cGGCAGAAAAAGCGAAAATCGAACTGTC‑TTCCGCTCAGCAGACCGACGTTAACCTGCCATACATCACTGCAGACGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATg                                                                      <  2:67081/251‑1 (MQ=255)                                                                                                                                                                                ggCAGAAAAAGCGAAAATCGAACTGTC‑TTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATgg                                                                     >  2:250957/1‑251 (MQ=255)                                                                                                                                                                                         aGCGAAAATCGACCTGCC‑TTCCGCGCAGCAGACCGAGTTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCagaaga                                                            <  2:351602/251‑1 (MQ=255)                                                                                                                                                                                                        gTC‑TTACGCTCAGCAGACCGACGTTAACCTGCCGTACGTCACTGCAGATGCGACCGGTCCGCAACACATGAACATCAAAGTGACTCGTGCGAAACCGGAAAGCCTGGTTGAAGCTCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCCCTGCAGGACGCTGGCATGTCCGAATCTGATATCGACGAAGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCGAAAGAACGTTGATGAGTTCt                                             >  2:322362/1‑251 (MQ=255)                                                                                                                                                                                                             tCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCttt                                           <  1:505120/249‑1 (MQ=255)                                                                                                                                                                                                             gccgcgcggcagACGGAGGTTAACCTGCCGTTCATCACTGTAGAGGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAGCGGGAAAGCGGTGTGGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCTGTATCTGATATCGACGACGTTAGCCTCGGTGGTGGTCATACTCGTGTTGCAATGGGTCAGAAGAAAGTTGCTGAGTTCttt                                           <  2:250132/241‑1 (MQ=255)                                                                                                                                                                                                              ccgaccagaacaagaaCGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGTCTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCGCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGt                                        <  2:126536/236‑1 (MQ=255)                                                                                                                                                                                                                   cagcagACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAaga                                   <  1:97799/251‑1 (MQ=255)                                                                                                                                                                                                                    agcagACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGat                                                                                                                   >  1:251555/1‑170 (MQ=255)                                                                                                                                                                                                                    agcagACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCCGACTCGTATGCCAATGGATCAGAAGAAAGTTGATGAGTTATTTGGTAAag                                    >  2:147082/1‑249 (MQ=255)                                                                                                                                                                                                                         aCCGACCTAAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACGCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCcgcg                             <  2:162467/251‑1 (MQ=255)                                                                                                                                                                                                                           cGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAAACTCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTAGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGt                            >  1:47029/1‑250 (MQ=255)                                                                                                                                                                                                                                ttAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGAAGGACGCTGGCCTATCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGGGTTCTTTGGTAAAGAGCCGCGTAAAGAc                      >  2:289113/1‑251 (MQ=255)                                                                                                                                                                                                                                gtAACCTGCCGTACATCACGGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGAc                      <  1:56740/250‑1 (MQ=255)                                                                                                                                                                                                                                      tGCCGTACATCACTCCAGATGCGACCGGTCCGAAACACCTGAAAATCAACGTGACTCGTGCGACACTGGAAAGACTGGTTGAAGATCTGGACAACAGTTCCATTGAGCCGCTGAAAGTTGCACTACAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTGAAgcctttaa                 >  2:329335/1‑243 (MQ=255)                                                                                                                                                                                                                                       gccgtACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCCATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAAcc               >  2:425172/5‑251 (MQ=255)                                                                                                                                                                                                                                          gtACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAACGACGTTAACCCgg            >  1:264975/2‑251 (MQ=255)                                                                                                                                                                                                                                            acaacACTGCATATGCGACCGTTCCGAAAAACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTTGTTGAAGATCTGGTAAACCTTTCCATTGAGCCGTTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTTAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTACCCCGGAc          <  2:333524/247‑1 (MQ=255)                                                                                                                                                                                                                                             cATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACg         <  1:425172/251‑1 (MQ=255)                                                                                                                                                                                                                                                 aCTGCATATGCGACCGGGCCGAAACACATGAACAGCAAAGTGAATCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCGTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGc     <  2:355893/251‑1 (MQ=255)                                                                                                                                                                                                                                                   ggCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTg   <  2:383068/250‑1 (MQ=255)                                                                                                                                                                                                                                                    gCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGt  <  1:505951/251‑1 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                      TCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGGGGTGAAGACTTCGACAGCCGTCTGATCAACTATCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTC‑TTCCGCTCAGCAGACCGACGTTAACCTGCCATACATCACTGCAGACGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGT  >  NC_000913/12782‑13274

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