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	92.7%	109.6 / ‑0.3	41	P284P (CCA→CCG)	dnaK	chaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base G (20/18);  total (22/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.22e-01

CTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGGGGTGAAGACTTCGACAGCCGTCTGATCAACTATCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACC‑TGCCATACATCACTGCAGACGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCG  >  NC_000913/12770‑13245                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                        cTATTATCGAACTCGAAGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCGACGGTGATACCCACCTGGGGGGTGAAGACTTCGACAGCCGTCTGATCAACTATATGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGAAATGCAGCGCCTGAAAGACGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGT‑AACCATGCCAAACATc                                                                                                                                                                                                                                   >  2:200924/1‑251 (MQ=255)           aaTCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATAACCACCTGGGTGGTGAAGACTTCGACCACCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACC‑TGCCGTACATCACTGCAGatg                                                                                                                                                                                                                         >  1:241709/1‑249 (MQ=255)            aTCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACC‑TGCCGTACATCACTGCAGatgc                                                                                                                                                                                                                        >  1:68465/1‑248 (MQ=255)                   aatttGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGAGCATGAAAGAAGCGGCAGAAAAAGAGAACATCGAACTGTCTTCCGCTCAGCAGACCGACGTTACCC‑TGCCGTAcctaactgaaactgtgaccgtt                                                                                                                                                                                                                 >  2:297256/4‑230 (MQ=255)                         aCGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGGGAAGACTTCGAAAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACC‑TGCCGTACATCACTGCAGATGCGACCGGTCCGaaa                                                                                                                                                                                                           <  1:231965/251‑1 (MQ=255)                                  aaaCCTGCGAAGTTCCGGCAACCAACGTGGATACCCACCTGGGGGTTGAAGACTGCGACTGCCTGCTGAACAACTGTCTGGTTGAAAAATACAAGAAAGATCAGTGCATGGACCGGAGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGACCTGCCTTCCGCTCAGCAGACCGACGTTAACC‑TGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAAc                                                                                                                                                                                                  <  2:34224/251‑1 (MQ=255)                                           aaGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACC‑TGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTg                                                                                                                                                                                         >  1:277867/1‑251 (MQ=255)                                                cgggCAACCATCGGTGATACCCACCGGGGTGGTGAAGACTTCGACAGCCGTCTGTTCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACC‑TGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTc                                                                                                                                                                                     <  2:314493/248‑1 (MQ=255)                                                    caaccaacGGTGAGACCCGCCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGGCTTCCGCTCAGCAGACCGACGTTAACC‑TGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCg                                                                                                                                                                                <  2:59135/251‑1 (MQ=255)                                                                         tGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGTGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACAGACGTTAACC‑TGCCGTACATAACTGCAGATGCGACCGGTCACAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTgc                                                                                                                                                            >  2:296050/1‑249 (MQ=255)                                                                                  aaGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACC‑TGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGt                                                                                                                                                   >  1:143963/1‑250 (MQ=255)                                                                                           acagccgtattattaacgaccTGGTGGAAGAATGTAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCTCTGGCAATGCAGCGCCTGAAAGAAGCGTCAGAAAAAGCGAAAACTGACCTGTCTTCCGCTCAGCATACCGACGTTAACC‑TGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTcc                                                                                                                                         <  2:129176/231‑1 (MQ=255)                                                                                            cAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACC‑TGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCa                                                                                                                                        <  1:329133/251‑1 (MQ=255)                                                                                                               cTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACC‑TGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGc                                                                                                                     >  1:113222/1‑251 (MQ=255)                                                                                                                       agaaTTCAAGAAAGATCAGGGCATTGACCCGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGGCTTCCGCTCAGCAGACCGACGTTAACC‑TGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAgg                                                                                                             >  2:102808/1‑251 (MQ=255)                                                                                                                           ttCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACC‑TGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGc                                                                                                         >  1:33573/1‑251 (MQ=255)                                                                                                                                           gCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGACAATCGAAATGTCTTCCGCTCAGCAGACCGACGTTAACC‑TGCCGAACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGCAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCt                                                                                         >  2:84174/1‑251 (MQ=255)                                                                                                                                                gACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTGCCGCGCAGCAGACCGTCGTTAACC‑TGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCATGACGCTGGCCTGTCCGTATCTGatat                                                                                    <  1:264097/251‑1 (MQ=255)                                                                                                                                                            gagccGCTGGAAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGTACTTTCTGCCGCCAAGCAGACCGACGTTAACC‑TGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACAGCAAAGTGACTCGGGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCGGTCTGTATCTGATATCGACGACGTTAt                                                                        <  2:149389/248‑1 (MQ=255)                                                                                                                                                            gATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACC‑TGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTAt                                                                        <  1:84174/251‑1 (MQ=255)                                                                                                                                                                       aaTGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACC‑TGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGAAGACGTTATCCTCGTtggtg                                                             >  2:77624/1‑251 (MQ=255)                                                                                                                                                                        aTGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAAAC‑TGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTtggtgg                                                            >  1:16373/1‑251 (MQ=255)                                                                                                                                                                          gCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACC‑TGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTc                                                          <  1:267498/251‑1 (MQ=255)                                                                                                                                                                                   gAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACC‑TGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGt                                                  >  1:127346/1‑250 (MQ=255)                                                                                                                                                                                   gAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACC‑TGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGt                                                  >  1:182568/1‑250 (MQ=255)                                                                                                                                                                                      agaagCGGCAGAAAACGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACC‑TGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGt                                                                                                 >  1:302871/1‑200 (MQ=255)                                                                                                                                                                                      agaagCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCTCTCAGCAGTCCGACGTTAACC‑TGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGt                                                                                                 <  2:302871/200‑1 (MQ=255)                                                                                                                                                                                          gCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACC‑TGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAAt                                          >  1:292992/1‑251 (MQ=255)                                                                                                                                                                                          gCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACC‑TGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAAt                                          <  1:58674/251‑1 (MQ=255)                                                                                                                                                                                          gCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACC‑TGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAAt                                          <  1:346991/251‑1 (MQ=255)                                                                                                                                                                                                 aaaaaGCGAAAAGCGAGCTGTCTTCCGCTCAGCAGACCGACGTTAACC‑TGCCGTACATCACTGCAGATGCGACCGCTGCGAAACACATGAACATCAAAGTTACGCGTGCGAAACTGGAAAGCCTGGTTGAATATCTGGTAATCCGTTCTATTGACGCATTGAAAGTTGCACTGCAGTGCGCTGGCCTGTCCAGATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAAAGGTTCag                                   <  1:249383/251‑1 (MQ=255)                                                                                                                                                                                                   aaaGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACC‑TGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCagaa                                 <  1:375824/251‑1 (MQ=255)                                                                                                                                                                                                      gCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACC‑TGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAgaagaa                              <  1:9093/251‑1 (MQ=255)                                                                                                                                                                                                         aaaaTCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACC‑TGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGt                           <  1:131090/251‑1 (MQ=255)                                                                                                                                                                                                          aaaTCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACC‑TGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTCATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGtt                          >  1:673/1‑251 (MQ=255)                                                                                                                                                                                                                       ttCCGCTCAGCAGACCGACGTTAACC‑TGCCACACATCACTGCAGACGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTg             >  1:238487/1‑251 (MQ=255)                                                                                                                                                                                                                        tCCGCTCAGCAGACCGACGTTAACC‑TGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCtt               <  1:106077/248‑1 (MQ=255)                                                                                                                                                                                                                        tCCGCTCAGCAGACCGACGTTAACC‑TGCCATACATCACTGCAGACGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCCGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCttt              <  1:200924/249‑1 (MQ=255)                                                                                                                                                                                                                          cGCTCAGCAGACCGACGTTAACC‑TGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATcga                                                                                 >  1:222527/1‑180 (MQ=255)                                                                                                                                                                                                                                  agACCGACGTTAACC‑TGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCcg  >  1:275851/1‑251 (MQ=255)                                                                                                                                                                                                                                  agACCGACGTTAACC‑TGCCGTACATCAATGCAGATGCGACCGGTCAGAAACACACAAACATAAAAGTGACTCGTGCGAAACTGGAAAGCCTGGGTGAAGAGCTGGTAAACCGTTCTATTGAGCCGCTGAAAGTTGCAATGCAGGAAGCTGGCCTGTCCGTATGTGCTAAAGACGACGTTATCATCGTAGGTGGGCAGCCTCGCAGTCCAATGGATCATAAGAATGTTGTTGTGTTCTTTGGTAAAGTGCcg  >  2:253235/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                        CTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGGGGTGAAGACTTCGACAGCCGTCTGATCAACTATCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACC‑TGCCATACATCACTGCAGACGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCG  >  NC_000913/12770‑13245

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