Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,287,710	0	G	A	26.2%	73.0 / 11.8	42	intergenic (+43/‑54)	ytcB/nrfA	protein YtcB/cytochrome c552 nitrite reductase
Reads supporting (aligned to +/- strand):  ref base G (30/1);  new base A (0/11);  total (30/12)
Fisher's exact test for biased strand distribution p-value = 2.80e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGAGGAACATGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAA‑CTCTAAAGTGGTATTTTACATG‑CACTTACAATTGA‑TTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATA‑GGGGCGAGCAATGT‑CATGACAGTG‑TAG‑GTGCGGTTACTGTC‑GTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAG  >  NC_000913/4287464‑4287955                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                      tCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTCATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGAGGAACATGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAA‑CTCTAAAGTGGTATTTTACATG‑CACTTACAATTGA‑TTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATA‑GGGGCGAGCAATGT‑CATGACAGTG‑TAG‑GAGCGGTTACTGTC‑GTgcaa                                                                                                                                                                                                                                                    >  2:506336/1‑247 (MQ=255)                             caggaaggtgaaTTTTTAAGGGTCATCAAAAAAGGTTAATTAATTTTTTAAGAAACATACAGTTATGCTTGTTGTTGAAAAAGAGGAAAATTCTGATTAA‑CTCTAAAGGGTTATTTTACATC‑CATTTCCAATTGA‑TTAAAGACAACATTCACAGTGTGGTTTTTTGTTCCACATA‑GGGGCGAGAAATTT‑CAGTCCGGTG‑TCG‑GGGCGTTTACTTTC‑GTAAAAAACAATAAAAACCGCCATTGCAACAATgg                                                                                                                                                                                                                       <  2:411713/241‑1 (MQ=255)                                               aaGGCCAGCAACAAAAGTTGTTTAATTCTTTGAGGAACATCCAGTTAGGCATGCTGTTGAAAAAGAGGAATATACTTACTAA‑CTCTAAAGTGGTATTTTACATG‑CACTTACAAT‑GACTTAAAGACAACATTCACAGTGTGTTTTTTTGTTACACATA‑GGGGCGAAAACTGC‑CATGACGGTG‑GAG‑GGGCGGTCACTGTC‑GTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGcc                                                                                                                                                                                                     <  1:434314/251‑1 (MQ=255)                                                                 tGATTAATTCTTTGAGGAACATGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAA‑CTCTAAAGTGGTATTTTACATG‑CACTTACAATTGA‑TTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATA‑GGGGCGAGCAATGT‑CATGACAGTG‑TAG‑GTGCGGTTACTGTC‑GTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATaaa                                                                                                                                                                                    >  1:117670/1‑250 (MQ=255)                                                                      aaTTCTTTGAGGAACATGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAA‑CTCTAAAGTGGTATTTTACATG‑CACTTACAATTGA‑TTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATA‑GGGGCGAGCAATGT‑CATGACAGTG‑TAG‑GTGCGGTTACTGTC‑GTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATaa                                                                                                                                                                               >  1:227854/1‑250 (MQ=255)                                                                            ttGAGGAACATGCAGTTATGAATGCTGTTGAAAAAGAGGAAGATACTGACTAA‑CTCTAAAGTGTTATTTTACATG‑CACTTACAATTGA‑TTAAAGACAACATTAACAGTGTGGTTATTTGTTCCACATA‑GGGGCGAGCAATTT‑CATGACAGTT‑TAG‑GTGCGGTTACTGCC‑GTAAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAacgcac                                                                                                                                                                         <  1:552664/250‑1 (MQ=255)                                                                                ggAACATGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAA‑CTCTAAAGTGGTATTTTACATG‑CACTTACAATTGA‑TTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATA‑GGGGCGAGCAATGT‑CATGACAGTG‑TAG‑GTGCGGTTACTGTC‑GTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGt                                                                                                                                                                    >  1:338567/1‑251 (MQ=255)                                                                                         aGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAA‑CTCTAAAGTGGTATTTTACATG‑CACTTACAATTGA‑TTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATA‑GGGGCGAGCAATGT‑CATGACAGGG‑TAG‑GTGCGGTTACTGTC‑GTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGc                                                                                                                                                           >  1:442895/1‑251 (MQ=255)                                                                                          gTTATGCATGCTGTTGAAAAAGAGGAAAATCCT‑ACTAAACTCTAAGGGGGTATTTTACAGG‑CCTTTAAATTTGA‑TTAAAGAAAACATCCAAAGTGGGGTTATTTGTTACAAATA‑GGGGCGAAAAAGGT‑CATAACGGTG‑TAG‑GGGCGGTTACTGGC‑GTAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCt                                                                                                                                                          <  1:45922/251‑1 (MQ=255)                                                                                             atgcatgcTGTTGAAAAAGAGGAAGATACTGACTAA‑CTCTAAAGTGGTATTTTACATG‑CACTTACAATTGA‑TTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATA‑GGGGCGAGCAATGT‑CATGACAGTG‑TAG‑GTGCGGTTACTGTC‑GTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTAt                                                                                                                                                       >  2:268060/1‑251 (MQ=255)                                                                                                aatgttctaaaaaaaaAGGAAGATATTTACTAA‑CTCAAAAGGGGTTTTTTACATG‑CATTTACAATAAT‑TTAAAGTCAACATCCATGTTGGTTTTTTTTGTAACAAATG‑GGGGCGAAAAATGT‑CATGAAAGTG‑TGT‑GTGCGTTTACGGTC‑GTAAAAAAAAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGa                                                                                                                                                    <  1:404340/241‑1 (MQ=255)                                                                                                                       tACTGACTAA‑CTCTAAAGTGGTATTTTACATG‑CACTTACAATTGA‑TTAAAGACAACATTCACAGTGTGTTTATTTGTTACACATA‑GGGGCGAGCAATGT‑CATGACAGTG‑TCG‑GTGCGGTTACTGTC‑GTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGCAGCCCCTATGACAAAGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCt                                                                                                                             >  1:332173/1‑251 (MQ=255)                                                                                                                          tGACTAA‑CTCTAAAGTGGTATTTTACATG‑CACTTACAATTGA‑TTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATA‑GGGGCGAGCAATGT‑CATGACAGTG‑TAG‑GTGCGGTTACTGTC‑GTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGtt                                                                                                                          >  2:8446/1‑251 (MQ=255)                                                                                                                           gtacaacccctAAAGTGGTATTTT‑CTTGGCATTTAAAATTGA‑TAAAAGAAAACATTCAAGGTGTGGTTATTTGT‑ACCCAAAAGGGGAGAGCAATTT‑CAT‑ACAGGGGTAG‑GTGCGGTTCTTGGC‑GTAAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGtt                                                                                                                          <  1:451296/242‑1 (MQ=255)                                                                                                                                      aaaGTGGTATTTTACATG‑CACTTACAATTGA‑TTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATA‑GGGGCGAGCAATGT‑CATGACAGTG‑TAG‑GTGCGGTTACTGTC‑GTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAAc                                                                                                                >  1:495425/1‑250 (MQ=255)                                                                                                                                         gTGGTATTTTACATG‑CACTTACAATTGA‑TTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATA‑GGGGCGCGCAATGT‑CATGACAGTG‑TAG‑GTGCGGTTACTGTC‑GTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACacac                                                                                                            >  1:556502/1‑248 (MQ=255)                                                                                                                                                          atttaaaaTTGA‑TAAAAGAAAAAATCCACGGTGTGTTTATTTGTAACAAAAG‑GGGGCGAGAAATGC‑CATGAAAGTG‑TAG‑GGGCGTTTACTGTC‑GTAAAAAACAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCaaaa                                                                                            <  1:221001/245‑1 (MQ=255)                                                                                                                                                            ttACAATTGA‑TTAAAGACAATATTCACAGTGTGGTTATTTGTTACACATA‑GGGGCGAGCAATGT‑CATGACAGTG‑TAG‑GTGCGGTTACTGTC‑GTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAAcc                                                                                          >  2:204804/1‑251 (MQ=255)                                                                                                                                                            ttACAATTGA‑TTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATA‑GGGGCGAGCAATGT‑CATGACAGTG‑AAG‑GTGCGGTTACTGTC‑GTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTATGTTCACGCTGAACAAACGGCTGCTCCCGCAAAAcc                                                                                          >  1:364769/1‑251 (MQ=255)                                                                                                                                                                 aattattaaatgtcaacaTTCACGGTTTGGTTATTTTTTACTCATA‑GGGGGGTGCAATGC‑TAGTACGTTG‑TGG‑GGGCGTTTATTGTG‑GTAAAAAACAAAAAAAACCGCCATTTCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGTTTCCTTTTTTCTTTTTCTCTTCTGTTCACGCTGAACAACCGGCTGCTCCCGCAAAACCTGTa                                                                                      <  2:501587/238‑1 (MQ=255)                                                                                                                                                                    gA‑TTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATA‑GGGGCGAGCAATGT‑CATGACAGTG‑TAG‑GTGCGGTTACTGTC‑GTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACtg                                                                                  >  1:31103/1‑251 (MQ=255)                                                                                                                                                                       ttAAAGACAACATTCACAGTGTGGTTATTTGTTACACATA‑GGGGCGAGCAATGT‑CATGACAGTG‑TAG‑GTGCGGTTACTGTC‑GTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACtgtg                                                                                >  1:527358/1‑251 (MQ=255)                                                                                                                                                                        tAAAGACAACATTCACAGTGTGGTTATTTGTTACACATA‑GGGGCGAGCAATGT‑CATGACAGTG‑TAG‑GTGCGGTTACTGTC‑GTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTgg                                                                               >  2:579197/1‑251 (MQ=255)                                                                                                                                                                        tAAAGACAACATTCACAGTGCGGGTATTTGTTACACATA‑GGGGAGAGCAATGT‑CATGACAGTG‑TAG‑GTGCGGTTACTGTC‑GTGAAAAACAATAAAAACAGCCATTGCAACATTGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGAACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGTAAAAACGGCTGCTCCCGCAAAACCTGTACCTGTgg                                                                               >  2:50564/1‑251 (MQ=255)                                                                                                                                                                            aacaacaTTCCAAGGGTGTTTTTTTGTTACACAAG‑GGGGGGAGCA‑TGTTAATGACAGTG‑TGG‑GTGGGTTTACTGCG‑GTAAAAAAAAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGc                                                                           <  1:440749/250‑1 (MQ=255)                                                                                                                                                                               aacaTTAACAGTGTGGTTATTTGTTACACATA‑GGGGCGAGCAATGT‑CATGACAGTG‑TAG‑GTGCGGTTACTGTC‑GTGAAAAACAATAAAAACCGCCATTGCAAAAATGGGGCAATTCGGATGCAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATATTAAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACACACGGCTGCTACCGCAAAACCTGTAAATGTGGAAGCga                                                                         >  1:334128/1‑250 (MQ=255)                                                                                                                                                                                      acaGTGTGGTTATTTGTTACACATA‑GGGGCGAGCAATGT‑CATGACAGTG‑TAG‑GTGCGGTTACTGTC‑GTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGa                                                                  >  1:497214/1‑250 (MQ=255)                                                                                                                                                                                      aaagtgggtttATTGGTTACAAAGA‑GGGGCAACAAATGT‑CATGACGGTG‑TGG‑GTGGGGTTACTGTC‑GTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGaa                                                                 <  1:268060/242‑1 (MQ=255)                                                                                                                                                                                                 aTTTGTTACACATA‑GGGGCGAGCAATGT‑CATGACAGTG‑TAG‑GTGCGGTTACTGTC‑GTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACTGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGcccc                                                      >  2:186555/1‑251 (MQ=255)                                                                                                                                                                                                    tGTTACACATA‑GGGGCGAGCAATGT‑CATGACAGTG‑TAG‑GTGCGGTTACTGTC‑GTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCACTTCGGATGAAGCCCATATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCAGCAAAACCTGTAACTGTGGAAGCAAAGAAAGACACCTTTGccctga                                                    >  2:8400/1‑247 (MQ=255)                                                                                                                                                                                                       tACACATA‑GGGGCGAGCAATGT‑CATGACAGTG‑TAG‑GTGCGGTTACTGTC‑GTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCgcagc                                                 >  1:158939/1‑250 (MQ=255)                                                                                                                                                                                                       tACACATA‑GGGGCGAGCAATGT‑CATGACAGTG‑TAG‑GTGCGGTTACTGTC‑GTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCgcagc                                                 >  1:187071/1‑250 (MQ=255)                                                                                                                                                                                                          acaTA‑GGGGCGAGCTACGT‑CATGACAGTG‑TAG‑AAGCGGTTACTGTC‑GTGAAAAAAAATAAAAACCGACATTGCAACAATGGGCCAAATCGCATGAACCACCTCTGACAAGGATAGATATAAACGCACGCCGAATCTTCAGATTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAACCGGCTGCTCCCGCAAAAACAGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCACGCAGCATcc                                             >  2:543561/1‑251 (MQ=255)                                                                                                                                                                                                            aTA‑GGGGCGAGCAATGT‑CATGACAGTG‑TAG‑GTGCGGTTACTGTC‑GTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCg                                           >  2:126575/1‑251 (MQ=255)                                                                                                                                                                                                              a‑GGGGCGAGCAATGT‑AATGACAGAG‑TAG‑GTGAGGTTACTGTC‑GAGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGACCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAACACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCCTCCCGAt                                         >  2:513952/1‑251 (MQ=255)                                                                                                                                                                                                                  ggCGAGCAATGT‑CATGCCAGTG‑TAG‑GTGCGGTTACTGTC‑GTGAAAAACAATAAAAACCGCCATTGAAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGACGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCa                                       >  2:128951/1‑250 (MQ=255)                                                                                                                                                                                                                    cGAGCAATGT‑CATGACAGTG‑TAG‑GTGCGGTTACTGTC‑GTGAACAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAat                                     >  2:271122/1‑250 (MQ=255)                                                                                                                                                                                                                     gAGCAATGT‑CATGACAGTG‑TAG‑GTGCGGTTACTGTC‑GTGAAAACCAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAatat                                   >  1:550858/1‑251 (MQ=255)                                                                                                                                                                                                                        cAATGT‑CATGACAGTG‑TAG‑GTGCGCTTACTCTC‑GTGAAAAACAATAAAAACCGACATAGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACCAGGATAAAAATAAACGCACACCGAATCTTCAGCTTATTGATTCCTTTTTTATTTTTCACTTCTGTTCCCGCTGAACAAACGACTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATAtctc                                >  2:268212/1‑251 (MQ=255)                                                                                                                                                                                                                               cATGACAGTG‑TTGTGTGCG‑TTACTCTCTGTAAAAAAAAATAAAACCCGCCATGACAACAATGGCGAAATTCGGATGAAGCCCCAATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAATCCTTTGCCCAGCAGCATCCCGATCAATATCGCTCCt                            <  2:18297/250‑1 (MQ=255)                                                                                                                                                                                                                                                     ttcTGTC‑GTAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGt       <  1:136692/248‑1 (MQ=255)                                                                                                                                                                                                                                                        tGTC‑GTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTc      >  1:280193/1‑248 (MQ=255)                                                                                                                                                                                                                                                         gTC‑GTGAAACACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCACTATGACAAGGAAAAAAATAAACGCACGCCGTATCTTCAGCTTATTGCTTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCCCAAAACCTGTAACTGTGGAAGCGCAGACTGAAACCTTTGCCCCCCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCagcg  >  2:46941/1‑249 (MQ=255)                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                      TCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGAGGAACATGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAA‑CTCTAAAGTGGTATTTTACATG‑CACTTACAATTGA‑TTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATA‑GGGGCGAGCAATGT‑CATGACAGTG‑TAG‑GTGCGGTTACTGTC‑GTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAG  >  NC_000913/4287464‑4287955

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