Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	729,707	G→A	G42E (GGG→GAG)	rhsC →	rhs element protein RhsC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	729,707	0	G	A	94.9%	95.0 / ‑6.6	40	G42E (GGG→GAG)	rhsC	rhs element protein RhsC
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (22/15);  total (25/15)
Fisher's exact test for biased strand distribution p-value = 5.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCG‑GGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCAGAAGCCTGTATTTTG  >  NC_000913/729457‑729952                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                      tGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCcggcg‑g                                                                                                                                                                                                                                                       >  1:274619/1‑251 (MQ=255)           ataAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCG‑AGGTGACGTc                                                                                                                                                                                                                                              <  1:379187/250‑1 (MQ=255)                     tCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCG‑AGGTGACGTCCGGCCCTCCgg                                                                                                                                                                                                                                   >  1:55874/1‑251 (MQ=255)                      ccTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGGACGCATCGGTGCCCCCCCCCGCGTGGCCTGTTCGGGGTGGCCCCGCG‑CGGGGACGCCCGGccccccggg                                                                                                                                                                                                                                  >  2:371396/1‑244 (MQ=17)                         aaTTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGAAGCATTGTCCAGGGGTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCCGTTCGGTGTGGCCCGGCG‑GGGGGACGGCCGGCCCTCCGGTCa                                                                                                                                                                                                                                >  2:253516/1‑250 (MQ=17)                             gTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCG‑AGGTGACGTCCGGCCATCCGGTCAATCCg                                                                                                                                                                                                                           <  1:351195/251‑1 (MQ=255)                                       gACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCG‑AGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCg                                                                                                                                                                                                                    <  2:231576/248‑1 (MQ=255)                                          aGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCG‑AGGTGACGTCCGGCCCTCCGGTCAATCCGCCGCTCGGGGCa                                                                                                                                                                                                               >  1:157615/1‑250 (MQ=255)                                           gATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGGTCAGCCGGGGTACGCATCGGTGCCCCCCCCGGCGTGGCCTGTTCGGTGTGCCCCGGCG‑AGGGGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCaa                                                                                                                                                                                                              >  2:46371/1‑250 (MQ=255)                                                                          gTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCCGGGTTCAGCCGGGGTACGCCTCGGTGCCCCCACCGGCGTGGCCTGTTCGGGGGGCCCCGGCG‑AGGGGACGTCCCGCCATCCGGTCAATCCGCTGCTCGGGGCAAAGGTCCTTCCCCGTGGAACCGACATCGCCCTg                                                                                                                                                                              >  1:346490/1‑251 (MQ=14)                                                                             gggTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCCTTGTCCCGGGGTCCGCCCGGGGACGCACCGGTGCCCCCACCGGCGTGGCCTGTTCGGGGGGGCCCGGCGAGGGGGA‑GTCCGGCCATCCGGTCCCTCCGCCGCTCGGTGCCAAAGTCCCTCCCGGGGAAACCGACCACGCCCTgccc                                                                                                                                                                           >  2:305286/1‑251 (MQ=14)                                                                                        gTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCG‑AGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTg                                                                                                                                                                  <  1:307911/249‑1 (MQ=255)                                                                                                  aaaTAACCCGATAAAACACGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGGTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCCGCG‑AGGTGACGTCCGGCCCTCCGGTCAATCCCCCGCCCGGTGCAAAGGGCCTTCCCGGTGAAACCCCCATCGCCCCGCCCCGCCCCCCGCCGTTCCTtctc                                                                                                                                                      >  2:40024/1‑251 (MQ=14)                                                                                                    aTAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGAGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGAATTGTCCAGGGTTCCACCGGGGTACGAATCGGTGGCCCCACCGGCGTGGCCTGTTCGGAGTGCCCCCGCG‑AGGTGCAGTCCGGGCCTCCGGGCAATCCGCTGCTCCGTGCAAAGGTCCTTCCCGGTGAAACCCCCATCGCCCTGCCCGGCCCGCTCCCGTTCCTtctccc                                                                                                                                                    >  1:45209/1‑249 (MQ=14)                                                                                                                   aaGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCCGTTCGGTGTGCCCCGGCG‑AGGGGACGTCCGGCCATCCGGTCCATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAAACGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCCCCTAcagc                                                                                                                                       >  1:219279/1‑249 (MQ=255)                                                                                                                           gCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCG‑AGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGa                                                                                                                             <  1:355047/251‑1 (MQ=255)                                                                                                                                  ggCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGCGCCCCGGCG‑AGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAAc                                                                                                                      <  1:62607/250‑1 (MQ=255)                                                                                                                                   gCGGAAAACCGGCGGCGCGTCAGTGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCG‑AGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACg                                                                                                                     <  1:35660/251‑1 (MQ=255)                                                                                                                                    cGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCG‑AGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGc                                                                                                                    <  1:46371/251‑1 (MQ=255)                                                                                                                                           ccGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCG‑AGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGGTCATTCTCTCCCGCCCCTACAGCAGTTACCGGACAAAAACGCCCGCGcc                                                                                                             >  1:277023/1‑251 (MQ=255)                                                                                                                                                           gCGACATGACGCAGTATGGCGGTAGCATTTTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCCGTTCGGTGTGCCCCGGCG‑AGGGGACGTCCGGCCATCCGGTCCATCCGCTTCCCGGTGCAAAAGGCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCCCCTCCAGCAGTTACCGGACAAAAACGCCCGCGCCGGGGGGGGGCCTCGGc                                                                                             >  2:4827/1‑251 (MQ=14)                                                                                                                                                               cATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCG‑AGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACa                                                                                                                                                                                      <  1:307055/158‑1 (MQ=255)                                                                                                                                                                             gCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCG‑AGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCt                                                                           <  2:274619/251‑1 (MQ=255)                                                                                                                                                                                gTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGGGTGCCCCGGCG‑AGGGGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCGGGAAAATGCCTGCg                                                                        >  2:376654/1‑251 (MQ=255)                                                                                                                                                                                  aGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCG‑AGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCgg                                                                       >  1:276875/1‑250 (MQ=255)                                                                                                                                                                                       tGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCG‑AGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCCGGCAAAAGCCTGCGGACATcc                                                                 >  2:244647/1‑251 (MQ=255)                                                                                                                                                                                        gTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCG‑AGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATccc                                                                >  2:248788/1‑250 (MQ=255)                                                                                                                                                                                                ttCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCG‑AGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGGGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTAccgc                                                        >  2:167881/1‑248 (MQ=255)                                                                                                                                                                                                      cGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCG‑AGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGGGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGGGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTTcgcg                                                  >  2:63849/1‑251 (MQ=255)                                                                                                                                                                                                            ccGCACCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCG‑AGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAaca                                            <  2:28601/250‑1 (MQ=25)                                                                                                                                                                                                            aCGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCG‑AGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCCTTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCATCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAac                                             >  1:202885/1‑250 (MQ=255)                                                                                                                                                                                                              gCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCG‑AGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGCAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGAGATCCGCTTACCGCTGCGCGATAacacc                                          >  2:326459/1‑250 (MQ=255)                                                                                                                                                                                                                  cGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCG‑AGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTg                                       <  1:202301/250‑1 (MQ=25)                                                                                                                                                                                                                        ccccACCGGCGTGGCCTGTTCGGTGTGCCCCGGCG‑AGGGGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGACCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTCCCGGACAAAAACGCCCGCGCCGGGGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGCTATCCGCTTACAGCTGCGCGATAACACACTGATACTcc                                >  2:1212/1‑250 (MQ=255)                                                                                                                                                                                                                                 ggTGGCCTGTTCGGTGTGCCCCGGCG‑AGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAAcg                       <  2:55874/250‑1 (MQ=21)                                                                                                                                                                                                                                    ggCCTGTTCGGTGTGCCCCGGCG‑AGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAAcggc                     <  2:247931/250‑1 (MQ=25)                                                                                                                                                                                                                                            cGGTGTGCCCCGGCG‑AGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCAggatc            >  2:294288/1‑247 (MQ=21)                                                                                                                                                                                                                                             ggTGTGCCCCGGCG‑AGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGAAACTCAGTGATAACGGCGGCAGAAGcc           >  1:146108/1‑251 (MQ=21)                                                                                                                                                                                                                                                gtgCCCCGGCG‑AGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCCGAAGCCTGt        >  1:114685/1‑251 (MQ=21)                                                                                                                                                                                                                                                   ccccGGCG‑AGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCAGAAGCCTGTAtt     >  1:78880/1‑251 (MQ=21)                                                                                                                                                                                                                                                   ccccGGCG‑AGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCAGAAGCCTGTAtt     <  1:142773/251‑1 (MQ=21)                                                                                                                                                                                                                                                      cggcg‑aGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCAGAAGCCTGTAttttt  >  2:289095/1‑250 (MQ=21)                                                                                                                                                                                                                                                      cggcg‑aGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGATGCGAGCTAACACACTGATACTCAGGGATAACGGCGGCAGCAGCCTGTActttt  >  2:289080/1‑250 (MQ=11)                                                                                                                                                                                                                                                      cggcg‑aGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGACCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCatgacgcctgagtttg  >  2:339992/1‑236 (MQ=255)                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                      TGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCG‑GGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCAGAAGCCTGTATTTTG  >  NC_000913/729457‑729952

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