Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,556,207	0	T	C	27.9%	132.9 / 14.6	61	S215G (AGT→GGT)	rtcA	RNA 3'‑terminal phosphate cyclase
Reads supporting (aligned to +/- strand):  ref base T (39/5);  new base C (1/16);  total (40/21)
Fisher's exact test for biased strand distribution p-value = 1.53e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGG‑TGATATTTTCACTTTCG‑ACTTCAAG‑CGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑C‑TCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGC  >  NC_000913/3555961‑3556456                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                           cgcccgccAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGG‑TGATATTTTCACTTTCG‑ACTTCAAG‑CGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑C‑Tcccc                                                                                                                                                                                                                                                        <  1:436641/251‑1 (MQ=255)     cgccAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAGCTCTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGG‑GGATAGTTGCACTTTCG‑ACTTCAAG‑CGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGGTCATGCAGGTTAAAA‑C‑CCCCC‑Gccc                                                                                                                                                                                                                                                   >  2:240517/1‑250 (MQ=255)                 cGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGG‑TGATATTTTCACTTTCG‑ACTTCAAG‑CGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑C‑TCCCC‑GCCAGTGTAGCGAttt                                                                                                                                                                                                                                       >  2:605206/1‑251 (MQ=255)                  ggTAGCACCCGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGG‑TGATATTTTCACTTTCG‑ACTTCAAG‑CGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑C‑TCCCC‑GCCAGTGTAGCGATTTc                                                                                                                                                                                                                                      >  1:125862/1‑251 (MQ=255)                     aGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGG‑TGATATTTTCACTTTCG‑ACTTCAAG‑CGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑C‑TCCCC‑GCCAGTGTAGCGATTTCACg                                                                                                                                                                                                                                   >  1:248436/1‑251 (MQ=255)                             ttGGTCAGCAAAATTTCCCCCCCCCCCCGCCGTGTTTCCCAGGTGGCTTTCCCCCTCTTCCCCAAACGGGGCCCCAACCACCCCGGCATGGCCGCGTTTTTCACCGAC‑AAAAAAAAGGGGTCGGGTAATTTTTTCCCTTCCG‑CCTCAAAG‑CGAGAGGGTATTCCCCGGCCCCTGGCCGCGCGGCAGTTTAGAAATTTTCTGTCCAGGAAGGAAAAAC‑C‑CCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCa                                                                                                                                                                                                                            <  1:308332/250‑1 (MQ=255)                              tggtaagcaaaatttccccccACCCCCCCCGGGCTTGCCAGGTAGCGTTCCCCCTCTTCCACCACCTGGCCCCCAACCACCTCGGCATTGACGCGCTTTCCACCGACAAAAAAAAAGCTTTGGT‑TGATTTTTTCCCTTTCG‑ACTCCAAG‑CGAGACTTTATTCCCCGGCCCCGGGTCGCGCGGCAGTTTATAAATATTCTGTCCATGCAGGGAAAAC‑C‑CCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCaa                                                                                                                                                                                                                           <  1:138590/234‑1 (MQ=255)                              tGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGG‑TGATATTTTCACTTTCG‑ACTTCAAG‑CGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑C‑TCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAAc                                                                                                                                                                                                                          >  2:34886/1‑251 (MQ=255)                                gTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGG‑TGATATTTTCACTTTCG‑ACTTCAAG‑CGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑C‑TCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACAt                                                                                                                                                                                                                        >  2:508649/1‑251 (MQ=255)                                gTCAGAAAGATTTCCCCCCACTCCCGCCTTGCTTGCCAGGTAGCTTTCCCCCTTTTTCCCCAACTGGCCCGAAACCCCCGCGGCAGTGACGCGCTTTCCCCCGCCGAAAAAAAAGCTTTCGG‑TTTTATTTTC‑CTTTCGGACTTCAAG‑CGAGACGTTATTCCCGGGCCCCGGGCCGCGGGGCAGTTTATGAATATTTTGTTCATGCAGGAAAAAC‑C‑CCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACAt                                                                                                                                                                                                                        <  1:642949/251‑1 (MQ=255)                                                       ccgccgTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGG‑TGATATTTTCACTTTCG‑ACTTCAAG‑CGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑C‑TCCCC‑GCCAGTGTAGCGAATTCACGCTCAGCAACATGGCGCGGCACACCAGCTAAt                                                                                                                                                                                                    >  2:634829/1‑248 (MQ=255)                                                         gccgTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGG‑TGATATTTTCACTTTCG‑ACTTCAAG‑CGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑C‑TCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAAt                                                                                                                                                                                                    >  2:340815/1‑246 (MQ=255)                                                         gccgTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGG‑TGATATTTTCACTTTCG‑ACTTCAAG‑CGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑C‑TCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAAc                                                                                                                                                                                               >  1:630538/1‑251 (MQ=255)                                                           cgTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGG‑TGATATTTTCACTTTCG‑ACTTCAAG‑CGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑C‑TCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACtt                                                                                                                                                                                             >  1:213341/1‑251 (MQ=255)                                                                  gCCAGGGAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGG‑TGATATTTTCACTTTCG‑ACTTCAAG‑CGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑C‑TCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACg                                                                                                                                                                                      >  2:466788/1‑251 (MQ=255)                                                                  cccAGGTGGCTTTTCCCCTCTTCCACCACCTGTGCCGCACCCCCCTGGCCACTGCCCCGCTTTCCACCAAGGAAAAAAAAGCGTGCGG‑TGATATTT‑CCCTTTCC‑GGTTCAAGGCGAGAGGTTATTCCCCGGCCCCTGGCCGCGCGGCGGGTTAGAAATTTTCTGTTCATGCAGGAAAAAC‑C‑TCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACg                                                                                                                                                                                      <  2:60126/250‑1 (MQ=255)                                                                     aGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGG‑TGATATTTTCACTTTCG‑ACTTCAAG‑CGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑C‑TCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCAt                                                                                                                                                                                   >  1:129401/1‑251 (MQ=255)                                                                       gTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGG‑TGATATTTTCACTTTCG‑ACTTCAAG‑CGAGCCGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑C‑TCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGCCACACCAGCTAAAAGAACTTCTCCACGCATCt                                                                                                                                                                                 >  2:267308/1‑251 (MQ=255)                                                                               ccccctcttcccccaacgttgCCGAAACCCCCTCGGCACTGCCGCGTTTTCCCCCGACGAAAAAAAAGCGTTGGG‑TGATATTTTCCCTTTCG‑ACTTAAGG‑GGAGACGGTATTCCCCGGCCCCGGGTCGCGCGGCAGGTTATAAATTTTCTGTTAATGCAGGAAAAAC‑C‑CCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAAt                                                                                                                                                                          <  1:376488/231‑1 (MQ=255)                                                                                         cACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGG‑TGATATTTTCACTTTCG‑ACTTCAAG‑CGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑C‑TCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCgcgcc                                                                                                                                                               >  2:481961/1‑250 (MQ=255)                                                                                         cACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGG‑TGATATTTTCACTTTCG‑ACTTCAAG‑CGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑C‑TCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCt                                                                                                                                                               >  1:76195/1‑251 (MQ=255)                                                                                            cAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGG‑TGATATTTTCACTTTCG‑ACTTCAAG‑CGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑C‑TCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCCCGc                                                                                                                                                            >  2:603741/1‑251 (MQ=255)                                                                                                       caccccccccGGCACTGACGCGCTTTCCCCCGACGAAAAAAAGGCGTGCGG‑TGATTTTTTCCCTTCCG‑ACTTCAAG‑CGAGACGTTATTCCCCGGCCCCTGGCCGCGCGGCAGTTTATGAATATTCTGTTCATGCAGGAAAAAC‑C‑CCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAg                                                                                                                                                 <  1:252799/242‑1 (MQ=255)                                                                                                       caacccctcggctctgggcccgcTTTCCACCGACGAAAAAAAAGCGTTCGG‑TGATTTTTTCCTTTTCG‑CCTCCAGG‑CGGGACGGTATCCCCCGGCCCGGGGCCGCGCGGCAGGTTATAAATATTCTGTT‑AAGAGGGGAAAAAAC‑CCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAg                                                                                                                                                 <  1:582186/231‑1 (MQ=255)                                                                                                            accTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGG‑TGATATTTTCACTTTCG‑ACTTCAAG‑CGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑C‑TCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGtt                                                                                                                                            >  1:153782/1‑251 (MQ=255)                                                                                                                       gACGCGCTTTTCACCGACGACAAAAAAGCGTTCGG‑TGATATTTTCACTTTCG‑ACTTCAAG‑CGAGACGGTATTACCCAGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑C‑TCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGcc                                                                                                                                  >  2:493908/1‑250 (MQ=255)                                                                                                                              tttccACCGACGAAAAAAAAGCGTTGGG‑TGATTTTTTCACTTCCG‑ACATCAAG‑CGAGGCGTTATTCCCCGGCCCCTTGCCGCGCGGAAGGTTATGAATATTCTGTTAATGGAGGAAAAAC‑C‑CCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCg                                                                                                                           <  2:213341/246‑1 (MQ=255)                                                                                                                                   aCCGACGACAAAAAAGCGTTCGG‑TGATATTTTCACTTTCG‑ACTTCAAG‑CGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑C‑TCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACtt                                                                                                                     >  1:288622/1‑251 (MQ=255)                                                                                                                                      gacgacAAAAAAGCGTTCGG‑TGATATTTTCACTTTCG‑ACTTCAAG‑CGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑C‑TCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCg                                                                                                                  >  2:541778/1‑251 (MQ=255)                                                                                                                                           aaaaaagGCGTTCGG‑TGATATTTCCCCTTCCG‑ACTACAGG‑CGAGACGTTATTCCCCGGCCCTTGGTCGCGCGGCAGTTTATGAATTTTCTTTTCATGCAGGAAAAAT‑C‑CCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGcc                                                                                                             <  1:419251/250‑1 (MQ=255)                                                                                                                                             aaaaaGCGTTCGG‑TGATATTTTCACTTTCG‑ACTTCAAG‑CGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑C‑Tcccc                                                                                                                                                                                                                                                        >  1:184123/1‑111 (MQ=255)                                                                                                                                             aaaaaGCGTTCGG‑TGATAGTTTCACTTTCG‑ACTTCAAG‑CGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCATGGAAAAA‑C‑Tcccc                                                                                                                                                                                                                                                        <  2:184123/111‑1 (MQ=255)                                                                                                                                              aaaaGCGTTCGG‑TGATATTTTCACTTTCG‑ACTTCAAG‑CGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑C‑TCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGccacc                                                                                                          >  2:83682/1‑251 (MQ=255)                                                                                                                                               aaaagcttcgggttattttttCACTTTCC‑ACCCAAAG‑CGAAACGTTATGCCCCGGCCCCTGGTCGCGCGGCAGTTTATAATTATGCTTTCCATGC‑GGGAAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGccacc                                                                                                          <  1:337424/234‑1 (MQ=255)                                                                                                                                               aaaGCGTTCGG‑TGATATTTTCACTTTCG‑ACTTCAAG‑CGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑C‑TCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGccacc                                                                                                          >  1:461415/1‑250 (MQ=255)                                                                                                                                                 aGCGTTCGG‑TGATATTTTCACTTTCG‑ACTTCAAG‑CGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑C‑TCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCcaccac                                                                                                        >  2:276191/1‑250 (MQ=255)                                                                                                                                                 aGCGTTCGG‑TGATATTTTCACTTTCG‑ACTTCAAG‑CGAGACGGTATTACCCGGACCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑C‑TCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGATGCcaccac                                                                                                        >  1:546985/1‑250 (MQ=255)                                                                                                                                                     ttCGG‑TGATATTTTCACTTTCG‑ACTTCAAG‑CGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑C‑TCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGc                                                                                                   >  2:441821/1‑251 (MQ=255)                                                                                                                                                     ttCGG‑TGATATTTTCACTTTCG‑ACTTCAAG‑CGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑C‑TCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGc                                                                                                   >  1:380792/1‑251 (MQ=255)                                                                                                                                                       cGG‑TGATATTTTCACTTTCG‑ACTTCAAG‑CGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑C‑TCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGATGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCt                                                                                                 >  2:383589/1‑251 (MQ=255)                                                                                                                                                           taatttttcctttttCG‑ACTCCAAC‑AGAGACGTTATTCCCCGGCCCCTGCTCGCGCGGCAGTTTATGAATTTTCTGGACATGAAGGAAAAAC‑C‑CCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTcc                                                                                               <  1:383589/238‑1 (MQ=255)                                                                                                                                                             atatTTTCACTTTCG‑ACTTCAAG‑CGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑C‑TCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGc                                                                                            >  2:300384/1‑251 (MQ=255)                                                                                                                                                                   tCACTTTCG‑ACTTCAAG‑CGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑C‑TCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCCAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGt                                                                                       >  2:589994/1‑250 (MQ=255)                                                                                                                                                                    cACTTTCG‑ACTTCAAG‑CGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑C‑TCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTa                                                                                      >  1:504369/1‑250 (MQ=255)                                                                                                                                                                           g‑ACTTCAAG‑CGAGACGGGATTCCCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑C‑TCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGt                                                                              >  1:127595/1‑251 (MQ=255)                                                                                                                                                                           caccttaaaG‑CAAGAGGGTATTCCCCGCCCCCGGGTCGCGCGGCAGGTTATGAATATTCTTTTCAGGCAGGAAAAAC‑C‑CCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAAcc                                                                                <  1:481961/243‑1 (MQ=255)                                                                                                                                                                                acAAG‑CGAGCCGTTATTCCCCGGCCCCGGGCCGCGCGGAAGGTAATGAATTTTCTTTCCATGGAGGAAAAAC‑C‑CCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTgg                                                                            <  1:276191/248‑1 (MQ=255)                                                                                                                                                                                               aTTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGCAAAA‑C‑TCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCCTCTGCACAATGTTCCCGCGCTCGCCAAGTTGCACGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACAGCgctgttt                                                            >  2:511551/1‑245 (MQ=255)                                                                                                                                                                                                      ggCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑C‑TCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCtgatg                                                     >  2:458443/1‑251 (MQ=255)                                                                                                                                                                                                             ggcggcggggcattttATGAATTTTCTGTTAAGGCAGGAAAAAA‑C‑CCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTcc                                               <  1:605206/236‑1 (MQ=255)                                                                                                                                                                                                                cgcgcgGCAGTTTATAAATATTCTTTTCATGAAGGAAAAAC‑C‑TCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTAtt                                           <  1:34886/251‑1 (MQ=255)                                                                                                                                                                                                                   gcgGCAGGTTATGAATATTCTGTCAATGCAGGAAAAAC‑C‑CCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTc                                        <  1:300384/251‑1 (MQ=255)                                                                                                                                                                                                                     ggCAGGTTATGAATATTCTGTTAATGCAGGGAAAAA‑C‑TCCCC‑GCCAGTGTAGCGTTTTCACGCTCAGCAACAAGACGCGGCACACCAGCTAATAGCACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTACACGATGCCACGGGCGAGACTTCCGTTGCCACCAAACCGCCTCCGGCACGGTAAAAACCGCGGCGAAACAGCGTGGTTTGCTGATGACTTCCTATTTTCGc                                      >  2:555086/1‑251 (MQ=255)                                                                                                                                                                                                                        aGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑C‑TCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGcc                                     >  1:287719/1‑249 (MQ=255)                                                                                                                                                                                                                         ggTTATGAATATTCTGTTCATGCAGGGAAAAA‑C‑TCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGCCTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCcagc                                  >  2:583053/1‑251 (MQ=255)                                                                                                                                                                                                                             atgatttttCTGTTCAGGCAGGAAAAAC‑C‑CCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCagcagc                               <  1:267308/243‑1 (MQ=255)                                                                                                                                                                                                                                 atatTCTGTTCATGAAGGAAAAAC‑C‑CCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTc                          <  1:458443/251‑1 (MQ=255)                                                                                                                                                                                                                                    ttCTGTTCATGCAGGGAAAAA‑C‑TCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTcc                         >  1:521512/1‑249 (MQ=255)                                                                                                                                                                                                                                                 gggAAAAA‑CCTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGAt           <  1:441821/251‑1 (MQ=255)                                                                                                                                                                                                                                                          c‑TCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGc  >  2:11732/1‑251 (MQ=255)                                                                                                                                                                                                                                                          c‑TCCCC‑GCCAGGGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGc  >  2:586376/1‑251 (MQ=255)                                                                                                                                                                                                                                                            tCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCCCCAGCACCCGGCGGCTAAAATCCGc  >  2:165485/1‑250 (MQ=255)                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                           CGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGG‑TGATATTTTCACTTTCG‑ACTTCAAG‑CGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑C‑TCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGC  >  NC_000913/3555961‑3556456

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