Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	13,014	A→G	P284P (CCA→CCG)	dnaK →	chaperone protein DnaK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	13,014	0	A	G	89.8%	121.8 / 5.3	49	P284P (CCA→CCG)	dnaK	chaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base A (3/2);  new base G (21/23);  total (24/25)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.46e-01

TTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGGGGTGAAGACTTCGACAGCCGTCTGATCAACTATCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTC‑TTCCGCTCAGCAGACCGACGTTAACCTGCCATACATCACTGCAGACGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAA  >  NC_000913/12768‑13259                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      ttttattatCGAAATCGACGAAGTTGACGGCGAAAAAACCCTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGGCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGCCCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTC‑TTCCGCTCAGCAGACCGACGTTAACCTGCCGTAc                                                                                                                                                                                                                                                    <  2:251452/247‑1 (MQ=255)        atCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGGGGTGAAGACTTCGACAGCCGTCTGATCAACTATCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTC‑TTCCGCTCAGCAGACCGACGTTAACCTGCCATACATCACTGc                                                                                                                                                                                                                                            >  1:67081/1‑251 (MQ=255)               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)                                                               gTGATACCCACCTGGGGGGTGAAGACTTCGACAGCCGTCTGATCAACTATTTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTC‑TTCCGCTCAGCAGACCGACGTTAACCTGCCATACAt                                                                                                                                                                                                                                                  >  2:383435/1‑190 (MQ=255)                                                               gTGATACCCACCTGGGGGGTGAAGACTTCGACAGCCGTCTGATCAACTATTTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTC‑TTCCGCTCAGCAGACCGACGTTAACCTGCCATACAt                                                                                                                                                                                                                                                  <  1:383435/190‑1 (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)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      TTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGGGGTGAAGACTTCGACAGCCGTCTGATCAACTATCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTC‑TTCCGCTCAGCAGACCGACGTTAACCTGCCATACATCACTGCAGACGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAA  >  NC_000913/12768‑13259

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