Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,556,207	0	T	C	47.5%	63.9 / 33.5	62	S215G (AGT→GGT)	rtcA	RNA 3'‑terminal phosphate cyclase
Reads supporting (aligned to +/- strand):  ref base T (25/7);  new base C (0/29);  total (25/37)
Fisher's exact test for biased strand distribution p-value = 3.85e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGC‑GTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGG‑CCCCTGGTCG‑CGCGGCAGGTTATGAATATTCTGTTCA‑TGCAGGGAAAA‑ACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACA  >  NC_000913/3556058‑3556355                                                                                                                                                           |                                                                                                                                                     tgCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGC‑GTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGG‑CCCCTGGTCG‑CGCGGCAGGTTATGAATATTCTGTTCA‑TGCAGGGAAAA‑ACTc                                                                                                                                                     <  2:417325/151‑1 (MQ=255)       aaCCCCCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGC‑GTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGG‑CCCCTGGTCG‑CGCGGCAGGTTATGAATATTCTGTTCA‑TGCAGGGAAAA‑ACTCCCCGcc                                                                                                                                               <  1:997954/151‑1 (MQ=255)       aaCCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGC‑GTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGG‑CCCCTGGTCG‑CGCGGCAGGTTATGAATATTCTGTTCA‑TGCAGGGAAAA‑ACTCCCCGcc                                                                                                                                               <  1:7259/151‑1 (MQ=255)            ccTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGC‑GTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGG‑CCCCTGGTCG‑CGCGGCAGGTTATGAATATTCTGTTCA‑TGCAGGGAAAA‑ACTCCCCGcc                                                                                                                                               <  2:1192624/146‑1 (MQ=255)            ccTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGC‑GTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGG‑CCCCTGGTCG‑CGCGGCAGGTTATGAATATTCTGTTCA‑TGCAGGGAAAA‑ACTCCCCGcc                                                                                                                                               >  1:1193105/1‑146 (MQ=255)                 gCACTGACGCGCTTTTCACCGACGACAAAAAAGC‑GTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGG‑CCCCTGGTCG‑CGCGGCAGGTTATGAATATTCTGTTCA‑TGCAGGGAAAA‑ACTCCCCGCCAGTGTAGCGa                                                                                                                                     >  2:1233612/1‑151 (MQ=255)                 gCACTGACGCGCTTTTCACCGACGACAAAAAAGC‑GTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGG‑CCCCTGGTCG‑CGCGGCAGGTTATGAATATTCTGTTCA‑TGCAGGGAAAA‑ACTCCCCGCCAGTGTAGCGa                                                                                                                                     >  1:1070338/1‑151 (MQ=255)                  cACTGACGCGTTTTTCACCGACGACAAAAAAGC‑GTTCGGTAATATTTTCTCTTCCGTCTTCAAGCGAGACGGTATCCCCCGG‑CCCCTGGTCG‑CGCGGCAGTTTATAAATATTTTGTTCA‑GGCAGGAAAAA‑CCCCCCCGCCAGTGTAGCGa                                                                                                                                     <  2:1119164/150‑1 (MQ=255)                          ctttttcccccGACAAAAAAAAGGC‑GTCCGGTTATTTTTTCCTTTTCGCCTTAAAGCGAGAGGGTTTTCCCCGC‑CCCCGGGCCG‑CGCGGCAGTTTATAAATTTTCTTTTCA‑TGCGGGAAAAA‑CCCCCCCGCCAGTGTAGCGATTTCACGc                                                                                                                             <  1:127629/141‑1 (MQ=255)                              tttCACCGACGACAAAAAAGC‑GTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGG‑CCCCTGGTCG‑CGCGGCAGGTTATGAATATTCTGTTCA‑TGCAGGGAAAA‑ACTCCCCGCCAGTGTAGCGATTTCACGCTCAg                                                                                                                         >  2:1045981/1‑150 (MQ=255)                              tttCACCGACGACAAAAAAGC‑GTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGG‑CCCCTGGTCG‑CGCGGCAGGTTATGAATATTCTGTTCA‑TGCAGGGAAAA‑ACTCCCCGCCAGTGTAGCGATTTCACGCTCAg                                                                                                                         >  1:609108/1‑150 (MQ=255)                              ttcccccGACGAAAAAAAAGG‑GTTGGGTGATATTTTCCCTTTCGACTTAAAGGGAGACGGTATTCCCCGC‑CCCCGGGTCG‑CGCGGAGGGTTATAAATATTTTTTTCA‑TGCAGGAAAAA‑ACCCCCCGCCAGTGTAGCGATTTCACGCTCAGc                                                                                                                        <  1:1107434/146‑1 (MQ=255)                                  cccGACGA‑AAAAAAGT‑TGTCGGTGATTTTTCCCCTTTCGTCTTCAAGCGAGCGGGTATTCCCCGC‑CCCCTGGCCG‑CGCGGCGGTTTATGAATTTTCTGTACA‑TGCAGGGAAAA‑ACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACAt                                                                                                                   <  2:390797/150‑1 (MQ=255)                                   ccGACGAAAAAAAGGC‑GTTGGGTGATTTTTTCCCTTTCGACTTCAAGGGAGCCGGTTTTCCCCGG‑CCCCTGGTCG‑CGCGGCAGTTTATAAATTTTCTGTTCA‑TGCAGGAAAAA‑ACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACAt                                                                                                                   <  1:752376/151‑1 (MQ=255)                                   ccGACAAAAAAAAAGC‑GTTGGGTTATTTTTTCTCTTTCGCCTTCAAGCGAGGCGTTTTTCCCCGC‑CCCCTGGTCG‑CGCGGCAGGTTATAAATTTTCTGTCCT‑TGCAGGAAAAA‑CCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACAt                                                                                                                   <  1:992094/151‑1 (MQ=255)                                    ccgacaaaaaaaaGCTGTGGGGTGATTTTTCCCCTTCCGACTCAAAGCGAGAGGTTATTCCCCGC‑CCCCGGGGCG‑CGCGGCGGTTTATAAATGTTCTGTTCA‑TGCAGGAAAAA‑CCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACAt                                                                                                                   <  1:932937/144‑1 (MQ=255)                                         aaaaaaaaGC‑GTCCGGTGATTTTTCCCCTTCCGACTTAAAGCGAGCCGTTATCCCCCGC‑CCCCTGGGCG‑CGCGGCGGTTTATAAATATTCTGTTCA‑GGCAGGAAAAA‑CCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgcgc                                                                                                              <  2:106197/148‑1 (MQ=255)                                                     tgcGGTTATTTTTTCCTTTTCGACTTAAAGCGAAACGGTATTCCCCGC‑CCCCGGGCCG‑CGGGGCAGGTTATAAATTTTTTGTTCA‑TGCAGGGAAAA‑CCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCt                                                                                                  <  1:220181/149‑1 (MQ=255)                                                          tcatttttccccTTCCGATTCCAAGGGAGACGGTATTACCC‑G‑CCCCCGGTCGCCGCGCAAGTTTAGAAATATTCTATCCA‑GGCAGGAAAAA‑CCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAg                                                                                             <  2:28150/140‑1 (MQ=255)                                                                tttCCCTTGCGACTTCAAGGGAGAGGGTATCCCCCGG‑CCCCTGGCCG‑CGCGGCGGTTTAAAAATTTTCTTTTCA‑GGCAGGAAAAA‑CCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTtc                                                                                       <  1:525368/151‑1 (MQ=255)                                                                      ttccgccTTCAAGCAAGACGGTATTCCCCGC‑CCCCTGGCCG‑CGCGGCAGGTTATGAATATTCTGTTCA‑TGCAGGAAAAA‑CCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACg                                                                                 <  1:358663/145‑1 (MQ=255)                                                                               cAAGCGAGCCGGTTTTCCCCGC‑CCCCGGGGCG‑CGGGGCGGGTTAGAAATATTCTGTTCA‑GGCGGGAAAAA‑ACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGcac                                                                        <  1:614469/151‑1 (MQ=255)                                                                               cAAGCGAGACGTTATTCCCCGC‑CCCGTGGTCG‑CGCGGCAGTTTATAATTATTCTGTCCA‑GGCAGGAAAAA‑CCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGcac                                                                        <  1:756075/151‑1 (MQ=255)                                                                               aaggggagaGGGTTTTCCCCGC‑CCCCTGGTCG‑CGCGGGAGTTTATGATTTTTCTGTTCA‑GGCAGGAAAAA‑CCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGcac                                                                        <  1:1040235/146‑1 (MQ=255)                                                                               aaagggggaCGTTATTCCCCGG‑CCCCGGGTCG‑CGCGGCGGTTTATAATTTTTCTTTCCA‑TGAGGGGAAAA‑CCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGcac                                                                        <  1:687066/144‑1 (MQ=255)                                                                                aaGCGGAACGGTTTTCCCCGGCCCCCTGGCCG‑CCCGGGGGTTTAGGAAT‑TTCTTTTTAGTGGAGGAAACA‑CCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGcac                                                                        <  2:505447/151‑1 (MQ=255)                                                                                   cgagatgtaatcccccGC‑CCCGTGGTCG‑CGCGGCGGGTTTTAATTTTTCTTTTCA‑TGAAGTAAAAA‑CCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCCCACCAGCTAATAGAACTTCTCCACGCATCTGCACAATg                                                                    <  2:947605/138‑1 (MQ=255)                                                                                     ggcgggtttcccccGG‑CCCCGGGGCG‑CGGGGCGGTTTATAAATATTCTTTTCA‑GGAAGGAAAAA‑CCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGt                                                                   <  1:782979/139‑1 (MQ=255)                                                                                     aggcGGTTTTCCCCGC‑CCCCTGGCCG‑CGCGGCAGTTGATAAATTTTCTGTTCA‑AGCAGGAAAAA‑CCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGt                                                                   <  1:1138016/147‑1 (MQ=255)                                                                                       ccGGTATTCCCCGG‑CCCCGGGTCG‑CGCGGCAGTTTATGAATATTCTGTTCA‑TGCAGGAAAAA‑ACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGtt                                                                  <  1:1044102/148‑1 (MQ=255)                                                                                       acgtttttccccGC‑CCCCTGGTCG‑CGCGGCAGGTTTTAAATTTTCTGTTCA‑TGCAGGAAAAA‑CCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTcc                                                                <  1:1164862/142‑1 (MQ=255)                                                                                         ggTATTACCCGG‑CCCCTGGTCG‑CGCGGCAGGTTATGAATATTCTGTTCA‑TGCAGGGAAAA‑ACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcg                                                              >  1:944553/1‑151 (MQ=255)                                                                                          gTATTACCCGG‑CCCCTGGTCG‑CGCGGCAGGTTATGAATATTCTGTTCA‑TGCAGGGAAAA‑ACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcg                                                              >  2:796804/1‑150 (MQ=255)                                                                                          gTATTACCCGG‑CCCCTGGTCG‑CGCGGCAGGTTATGAATATTCTGTTCA‑TGCAGGGAAAA‑ACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcg                                                              >  2:107869/1‑150 (MQ=255)                                                                                          gTATTACCCGG‑CCCCTGGTCG‑CGCGGCAGGTTATGAATATTCTGTTCA‑TGCAGGGAAAA‑ACTCCCCGCCAGCGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcg                                                              >  2:926819/1‑150 (MQ=255)                                                                                          gTATTACCCGG‑CCCCTGGTCG‑CGCGGCAGGTTATGAATATTCTGTTCA‑TCCAGGGAAAA‑ACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATCTTCCcg                                                              >  2:1200991/1‑150 (MQ=255)                                                                                                           ggTCG‑CGCGGCAGGTTATGAATATTCTGTTCA‑TGCAGGGAAAA‑ACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGACACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCGAGTTGCaa                                             >  2:320527/1‑151 (MQ=255)                                                                                                                      cgttttagaaaTATTCTTTTCA‑GCCAGAAAATG‑ACACCCCGCCAGTGTAGCTATTTCACGCTCAGCAACATGGCGCTGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAAc                                   <  2:654057/142‑1 (MQ=255)                                                                                                                      cggtttATGATTTTTTTTTTCA‑TGGAGGAAAAA‑CCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAAc                                   <  1:886941/147‑1 (MQ=255)                                                                                                                      cgggTTTTAAATATTCTGTTCG‑TGCAGGGAAAA‑CCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCGCGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAAc                                   <  1:1073900/149‑1 (MQ=255)                                                                                                                       ggtttATGAATATTTTGTTCA‑GGCGGGAAAAA‑CCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAAc                                   <  1:846630/147‑1 (MQ=255)                                                                                                                       gggTTATAAATATTCTTTTAA‑TGGGGGAAAAA‑CCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAAc                                   <  1:369387/149‑1 (MQ=255)                                                                                                                        ggTTATGCATATTCTGTTCA‑TGCAGGGAAAA‑ACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGa                                 >  2:187818/1‑151 (MQ=255)                                                                                                                        ggTTATGAATATTCTGTTCA‑TGCAGGGAAAA‑ACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGa                                 >  2:1248446/1‑151 (MQ=255)                                                                                                                        agaaagaataatTCTTTTAA‑TG‑AGGGAAAACCCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAAc                                   <  2:565629/139‑1 (MQ=255)                                                                                                                         gTTATGAATTTTCTGTCGA‑GGAGGAGAAAA‑CCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGa                                 <  2:80853/150‑1 (MQ=255)                                                                                                                         gTTATGAATATTCTGTTCA‑TGCGGGAAAAA‑CCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGTTGTTAAAc                                   <  1:598242/148‑1 (MQ=255)                                                                                                                         gTTATGAATATTCTGTTCA‑TGCAGGGAAAA‑ACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGCGTTAAACGAt                                >  2:58854/1‑151 (MQ=255)                                                                                                                           tataAATTTTCTGTTCA‑TGCGGGGAAAA‑CCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGc                              <  1:630807/151‑1 (MQ=255)                                                                                                                           tATGAATATTCTGTTCA‑TGCAGGGAAAA‑ACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATg                               >  2:626341/1‑150 (MQ=255)                                                                                                                            aTGAATATTCTGTTCA‑TGCAGGGAAAA‑ACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGc                              >  1:476079/1‑150 (MQ=255)                                                                                                                            aTGAATATTCTGTTCA‑TGCAGGGAAAA‑ACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGc                              >  1:963186/1‑150 (MQ=255)                                                                                                                            aTGAATATTCTGTTCA‑TGCAGGGAAAA‑ACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGc                              >  2:648483/1‑150 (MQ=255)                                                                                                                            aTGAATATTCTGTTCA‑TGCAGGGAAAA‑ACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGc                              >  1:738819/1‑150 (MQ=255)                                                                                                                              gAATATTCTGTTCA‑TGCAGGGAAAA‑ACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCCCGCCAAGTTGCAAGGTGTTAAAc                                   >  1:565862/1‑143 (MQ=255)                                                                                                                                 tatTCTGTTCA‑TGCAGGGAAAA‑ACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACggg                        >  2:636578/1‑151 (MQ=255)                                                                                                                                  atTCTGTTCA‑TGCAGGGAAAA‑ACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACggg                        >  1:229439/1‑150 (MQ=255)                                                                                                                                       tttaa‑TGCAGGAAAAA‑CCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGAc                  <  2:979375/147‑1 (MQ=255)                                                                                                                                                   gAAAA‑ACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGAGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGcca       >  1:1234626/1‑151 (MQ=255)                                                                                                                                                    aaaa‑cCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGccac      <  1:1234133/151‑1 (MQ=255)                                                                                                                                                      aa‑aCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGccacca    >  1:431709/1‑151 (MQ=255)                                                                                                                                                         aCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACCCCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACcaca  >  2:442635/1‑151 (MQ=255)                                                                                                                                                           |                                                                                                                                                     TGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGC‑GTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGG‑CCCCTGGTCG‑CGCGGCAGGTTATGAATATTCTGTTCA‑TGCAGGGAAAA‑ACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACA  >  NC_000913/3556058‑3556355

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