Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	224,989	A→T	noncoding (1219/1542 nt)	rrsH →	16S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	224,989	0	A	T	90.2%	10.6 / 20.0	92	noncoding (1219/1542 nt)	rrsH	16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (4/5);  new base T (21/62);  total (25/67)
Fisher's exact test for biased strand distribution p-value = 2.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

GTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGC‑GGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGGCGCATA  >  NC_000913/224854‑225018                                                                                                                                         |                              gTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGc                           <  2:873082/140‑1 (MQ=2) gTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGc                           <  2:386029/140‑1 (MQ=2) gTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGc                           <  2:228743/140‑1 (MQ=2)  ttGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAAACGGCGGGAGGTGGGGAAGACGTCAAATCATCATGCCCCTTATGACCTGGGCt                          >  2:592137/1‑140 (MQ=2)  ttGGGTTAAGTCCCCCAACGAGCGCACCCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGc                           <  2:1165549/140‑1 (MQ=2)   tGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTa                         <  2:1193238/140‑1 (MQ=2)   tGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTa                         <  2:4573/140‑1 (MQ=2)   tGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTa                         <  1:831289/140‑1 (MQ=2)   cgggTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTa                         <  2:607671/139‑1 (MQ=2)    gggTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGTGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTac                        <  1:321697/140‑1 (MQ=2)    gggTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTac                        <  2:1231169/140‑1 (MQ=2)    gggTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTac                        <  2:233155/140‑1 (MQ=2)    gggTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTac                        <  1:514722/140‑1 (MQ=2)    gggTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTac                        <  1:196907/140‑1 (MQ=2)    gggTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTac                        <  1:914195/140‑1 (MQ=2)    gggTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTac                        <  2:721512/140‑1 (MQ=2)     ggTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGC‑GGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTaca                       >  2:1602805/1‑140 (MQ=32)     ggTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTac                        <  1:1133627/140‑1 (MQ=2)     ggTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTac                        <  2:832369/140‑1 (MQ=2)      gTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacac                      <  2:1097507/140‑1 (MQ=2)      gTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCTGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacac                      <  1:502462/140‑1 (MQ=2)      gTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacac                      <  1:568283/140‑1 (MQ=2)      gTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacac                      <  2:14476/140‑1 (MQ=2)      gTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacac                      <  1:788965/140‑1 (MQ=2)      gTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacac                      <  1:330673/140‑1 (MQ=2)      gTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacac                      <  1:404948/140‑1 (MQ=2)      gTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacac                      <  2:519490/140‑1 (MQ=2)      gTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacac                      <  2:177081/140‑1 (MQ=2)       ttAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGC‑GGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTacaca                     >  1:1248413/1‑140 (MQ=33)       ttAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacaca                     <  2:644917/140‑1 (MQ=2)       ttAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacaca                     <  2:778466/140‑1 (MQ=2)       ttAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCGAATCATCATGCCCCTTATGACCTGGGCTacaca                     >  1:430534/1‑140 (MQ=2)       ttAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacaca                     >  1:515424/1‑140 (MQ=2)       ttAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacaca                     <  2:515424/140‑1 (MQ=2)        tAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacaca                     <  2:1004608/140‑1 (MQ=2)         aaGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCATATCATCATGCCCCTTATGACCTGGGCTACACACg                   >  2:53356/1‑140 (MQ=2)         aaGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACg                   <  2:665025/140‑1 (MQ=2)         aaGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACg                   <  2:468821/140‑1 (MQ=2)         aaGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACg                   <  1:861461/140‑1 (MQ=2)         aaGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACg                   >  1:196622/1‑140 (MQ=2)          aGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGC‑GGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGt                  >  2:1771625/1‑140 (MQ=32)          aGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTac                        >  1:356619/1‑135 (MQ=2)          aGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTac                        <  2:356619/135‑1 (MQ=2)          aGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGc                           >  2:626431/1‑131 (MQ=2)          aGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGc                           <  1:626431/131‑1 (MQ=2)          aGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTac                        >  1:1048760/1‑134 (MQ=2)          aGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTac                        <  2:1048760/134‑1 (MQ=2)          aGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGt                  >  2:1174757/1‑140 (MQ=2)          aGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGt                  >  1:358268/1‑140 (MQ=2)           gTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGta                 <  1:156688/140‑2 (MQ=2)           gTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGta                 <  2:1042529/140‑2 (MQ=2)           gTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGta                 <  2:910073/140‑2 (MQ=2)           gTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGta                 <  1:729707/140‑2 (MQ=2)           gTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGta                 <  1:828308/140‑2 (MQ=2)           gTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGta                 >  2:277854/1‑139 (MQ=2)           gTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGta                 <  1:1465920/140‑2 (MQ=2)           gTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGta                 <  1:1088887/140‑2 (MQ=2)           gTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGta                 >  1:969799/1‑139 (MQ=2)            tCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTa                         >  2:1116951/1‑132 (MQ=2)            tCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTa                         <  1:1116951/132‑1 (MQ=2)            tCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGta                 >  2:374426/1‑139 (MQ=2)            tCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGtac                >  2:770546/1‑138 (MQ=2)              ccGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTa              >  1:1150884/1‑140 (MQ=2)               cGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacac                      >  2:1001282/1‑131 (MQ=2)               cGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacac                      <  1:1001282/131‑1 (MQ=2)               cGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGtactac             >  2:372530/1‑140 (MQ=2)               cGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGtactac             >  2:106839/1‑140 (MQ=2)                gCAACGAGCGCAACCCTTATCTTTTGTTGCCAGC‑GGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACa            <  1:1540853/140‑1 (MQ=32)                gCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACa            <  1:526550/140‑1 (MQ=2)                gCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACa            <  2:1211099/140‑1 (MQ=2)                gCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACa            <  1:372530/140‑1 (MQ=2)                gCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACa            <  2:1078955/140‑1 (MQ=2)                gCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACa            <  2:342908/140‑1 (MQ=2)                gCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACa            <  1:71439/140‑1 (MQ=2)                 cAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACa            <  1:1067416/140‑1 (MQ=2)                 cAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACaa           >  2:394804/1‑140 (MQ=2)                    cGAGCGCAACCCTTATCTTTTGTTGCCAGC‑GGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATgg        <  1:1497198/140‑1 (MQ=32)                    cGAGCGCAACCCTTATCTTTTGTTGCCAGC‑GGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATgg        <  1:1721801/140‑1 (MQ=32)                     gAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACAATGgc       <  2:977036/140‑1 (MQ=2)                     gAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACAATGgc       <  2:948868/140‑1 (MQ=2)                     gAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACAATGgc       <  2:870557/140‑1 (MQ=2)                     gAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACAATGgc       <  2:869305/140‑1 (MQ=2)                     gAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACAATGgc       <  2:359469/140‑1 (MQ=2)                     gAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACAATGgc       <  2:466811/140‑1 (MQ=2)                     gAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACAATGgc       <  1:1131151/140‑1 (MQ=2)                      ggcgaAACCCTTATCTTTTGTTGCCAGC‑GGTCCGGCCGGGAACTCGAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGgcg      <  2:1243256/139‑1 (MQ=12)                      aGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACAATGgc       <  2:294938/140‑1 (MQ=2)                      aGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTac                        <  2:401745/122‑1 (MQ=2)                      aGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTac                        >  1:401745/1‑122 (MQ=2)                      aGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACAATGgcc      >  1:1193388/1‑139 (MQ=2)                          cAACCCTTATCTTTTGTTGCCAGC‑GGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGGCGCATa  <  2:1410256/140‑1 (MQ=32)                          cAACCCTTATCTTTTGTTGCCAGC‑GGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGGCGCATa  >  1:1739724/1‑140 (MQ=33)                              ccctATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACAATGgc       <  2:427183/128‑1 (MQ=2)                                                                                                                                         |                              GTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGC‑GGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGGCGCATA  >  NC_000913/224854‑225018

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