Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,401,459	0	T	C	27.6%	121.0 / 14.4	58	D625G (GAT→GGT)	csrD	targeting factor for csrBC sRNA degradation
Reads supporting (aligned to +/- strand):  ref base T (33/9);  new base C (0/16);  total (33/25)
Fisher's exact test for biased strand distribution p-value = 2.55e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GTCCAGAAAGGCTTAATGCATCAGCAGCATTCCTCGCCTGTTTGCAACCGCGTAAGGTCATTCATCTGCATATGTGCTGCAACAATCTGGCGAGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTA‑AACCGAGTATCTTTGTGAATATTTTTTCAC‑GTTAGTATCAAGTGGC‑TGTGAGGACGCGAAAAAA‑TCCCCTTG‑CCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAG  >  NC_000913/3401232‑3401679                                                                                                                                                                                                                                        |                                                                                                                                                                                                                              gTCCAGAAAGGCTTAATGCATCAGCAGCATTCCTCGCCTGTTTGCAACCGCGTAAGGTCATTCATCTGCATATGTGCTGCAACAATCTGGCGAGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTA‑AACCGAGTATCTTTGTGAATATTTTTTCAC‑GTTAGTATCAAGTGGC‑TGTGAGGACGCGAAAAAA‑TCCCCTTG‑ccc                                                                                                                                                                                                                    <  1:392401/238‑1 (MQ=255) gTCCAGAAAGGCTTAATGCATCAGCAGCATTCCTCGCCTGTTTGCAACCGCGTAAGGTCATTCATCTGCATATGTGCTGCAACAATCTGGCGAGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTA‑AACCGAGTATCTTTGTGAATATTTTTTCAC‑GTTAGTATAAAGTGGC‑TGTGAGGACGCGAAAAAA‑CCCCCTTG‑ccc                                                                                                                                                                                                                    <  2:245931/238‑1 (MQ=255)       aaaGGCTTAATGCATCAGCAGCATTCCTCGCCTGTTTGCAACCGCGTAAGGTCATTCATCTGCATATGTGCTGCAACAATCTGGCGAGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGCAAACGGATTAAACGGCAGGTTA‑AACCGAGTATCTTTGTGAATATTTTTTCAC‑GTTAGTATCAAGTGGC‑TGTGAGGACGCGAAAAAA‑TCCCCTTG‑CCCGCCGGTa                                                                                                                                                                                                             >  1:359273/1‑239 (MQ=255)            cTTAATGCATCAGCAGCATTCCTCGCCTGTTTGCAACCGCGTAAGGTCATTCATCTGCATATGTGCTGCAACAATCTGGCGAGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTA‑AACCGAGTATCTTTGTGAATATTTTTTCAC‑GTTAGTATCAAGTGGC‑TGTGAGGACGCGAAAAAA‑TCCCCTTG‑CCCGCCTGTAACAcc                                                                                                                                                                                                        >  2:74468/1‑239 (MQ=255)                     tCAGCAGCATTCCTCGCCTGTTTGCAACCGCGTAAGGTCATTCATCTGCATATGTGCTGCAACAATCTGGCGAGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTA‑AACCGAGTATCTTTGTGAATATTTTTTCAC‑GTTAGTATCAAGTGGC‑TGTGAGGACGCGAAAAAA‑CCCCCTTG‑CCCGCCTGTAACACCGCGCTGAAt                                                                                                                                                                                               <  1:74468/239‑1 (MQ=255)                                         ttgccacccGCGTAAGTTCATTCATCGGAATATGTGCTGCAAAAATCGGGCGAGCAGACAAGCCCACTCCCATGAGACCCAGCGCGCTTTATTCTACGTAAAAACGGATTAAACGGCAGGTAA‑AACCGAGTATCTTTGTAAATATTTTTTCAC‑GTTAGTATCAAGTGGC‑TGTGAGGACGCGAAAAAA‑CCCCCTTG‑ccc                                                                                                                                                                                                                    <  2:275918/191‑1 (MQ=255)                                                   cgTAAGGTCATTCATCTGCATATGTGCTGCAACAATCTGGCGAGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTA‑AACCGAGTATCTTTGTGAATATTTTTTCAC‑TTTAGTATCAAGTGGC‑TGTGAGGACGCGAAAAAA‑TCCCCTTg                                                                                                                                                                                                                        <  1:483231/185‑1 (MQ=255)                                                   cgTAAGGTCATTCATCTGCATATGTGCTGCAACAATCTGGCGAGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTA‑AACCGAGTATCTTTGTGAATATTTTTTCAC‑GTTAGTATCAAGTGGC‑TGTGAGGACGCGAAAAAA‑TCCCCTTg                                                                                                                                                                                                                        >  2:483231/1‑185 (MQ=255)                                                   cgTAAGGTCATTCATCTGCATATGTGCTGCAACAATCTGGCGAGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTA‑AACCGAGTATCTTTGTGAATATTTTTTCAC‑GTTAGTATCAAGTGGC‑TGTGAGGACGCGAAAAAA‑TCCCCTTG‑CCCGCCTGTAaca                                                                                                                                                                                                          >  1:275918/1‑198 (MQ=255)                                                                            gctgcAACAATCTGGCGAGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTA‑AACCGAGTATCTTTGTGAATATTTTTTCAC‑GTTAGTATCAAGTGGC‑TGTGAGGACGCGAAAAAA‑TCCCCTTG‑CCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAAACTGGGt                                                                                                                                            >  1:883532/1‑235 (MQ=255)                                                                                     aTCTGGCGAGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTA‑AACCGAGTATCTTTGTGAATATTTTTTCAC‑GTTAGTATCAAGTGGC‑TGTGAGGACGCGAAAAAA‑TCCCCTTG‑CCCGCCTGTAaca                                                                                                                                                                                                          <  1:148998/164‑1 (MQ=255)                                                                                     aTCTGGCGAGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTA‑AACCGAGTATCTTTGTGAATATTTTTTCAC‑GTTAGTATCAAGTGGC‑TGTGAGGACGCGAAAAAA‑TCCCCTTG‑CCCGCCTGTAaca                                                                                                                                                                                                          >  2:148998/1‑164 (MQ=255)                                                                                     aTCTGGCGAGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTA‑AACCGAGTATCTTTGTGAATATTTTTTCAC‑GTTAGTATCAAGTGGC‑TGTGAGGACGCGAAAAAA‑TCCCCTTG‑CCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAacgcacg                                                                                                                                                                >  1:316291/1‑206 (MQ=255)                                                                                     aTCTGGCGAGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTA‑AACCGAGTATCTTTGTGAATATTTTTTCAC‑GTTAGTATCAAGTGGC‑TGTGAGGACGCGAAAAAA‑TCCCCTTG‑CCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGa                                                                                                                                >  1:596202/1‑238 (MQ=255)                                                                                               cAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTA‑AACCGAGTATCTTTGTGAATATTTTTTCAC‑GTTAGTATCAAGTGGC‑TGTGAGGACGCGAAAAAA‑TCCCCTTG‑CCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTcc                                                                                                                      >  1:182323/1‑238 (MQ=255)                                                                                                      ccacaCCCCCATGAGCCGCGGCCAGAATTATTCTACGTAAAAACGGATTAACCGGCAGGTTA‑ACCCGAGTATCTTTGAGAATTTTTTTTCCC‑GTTATTATCAAGTGGC‑TGGGAGGACGCAAAAAAA‑CCCCCTGC‑CCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCg                                                                                                                                                       <  2:637271/197‑1 (MQ=255)                                                                                                                aTGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTA‑AACCGAGTATCTTTGTGAATATTTTTTCAC‑GTTAGTATCAAGTGGC‑TGTGAGGACGCGAAAAAA‑TCCCCTTG‑CCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTggg                                                                                                                                             >  1:637271/1‑198 (MQ=255)                                                                                                                 tGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTA‑AACCGAGTATCTTTGTGAATATTTTTTCAC‑GTTAGTATCAAGTGGC‑TGTGAGGACGCGAAAAAA‑TCCCCTTG‑ccc                                                                                                                                                                                                                    <  1:824303/126‑1 (MQ=255)                                                                                                                 tGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTA‑AACCGAGTATCTTTGTGAATATTTTTTCAC‑GTTAGTATCAAGTGGC‑TGTGAGGACGCGAAAAAA‑TCCCCTTG‑ccc                                                                                                                                                                                                                    >  2:824303/1‑126 (MQ=255)                                                                                                                      cGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTA‑AACCGAGTATCTTTGTGAATATTTTTTCAC‑GTTAGTATCAAGTGGC‑TGTGAGGACGCGAAAAAA‑TCCCCTTG‑CCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGCGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTgg                                                                                              >  1:645130/1‑239 (MQ=255)                                                                                                                      cGCACCTCGTATTATTCTACGGAAAATCGGATAAAGCGGAATGTTA‑AACCGATTATTTTTGGAAATTTTTTTCCCC‑GTTAGTATCAAGGGGC‑TGGGAGGACGCAAAAAAA‑CCCCCTTG‑CCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTg                                                                                               <  2:182323/238‑1 (MQ=255)                                                                                                                                attattCTGAGTGAAAACGGATGAAACCGCAAGTTAGAAGGTAGTTTCTTTAGAAATTTTTTTCCAC‑GTTATTA‑CAAGGGGCGTGTGAGGCCGGGAAAAAC‑CCCCCTCC‑CCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTtcaaagcc                                                                                     <  2:133202/239‑7 (MQ=255)                                                                                                                                attattCTACGTGAAAACGGATTAAACGGCAGTTTA‑ACCCGAGTATCTTTGTAAATATTTTTTCAC‑GTTAGTATCAAGGGGC‑TGTGAGGACGCAAAAAAA‑TCCCCTTG‑CCCGCCTGTAACACcgcg                                                                                                                                                                                                     <  2:758278/126‑1 (MQ=255)                                                                                                                                         cGTGAAAACGGATTAAACGGCAGGTTA‑AACCGAGTATCTTTGTGAATATTTTTTCAC‑GTTAGTATCAAGTGGC‑TGTGAGGACGCGAAAAAA‑TCCCCTTG‑CCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCaa                                                                           >  1:486522/1‑239 (MQ=255)                                                                                                                                         cGTGAAAACGGATTAAACGGCAGGTTA‑AACCGAGTATCTTTGTGAATATTTTTTCAC‑GTTAGTATCAAGTGGC‑TGTGAGGACGCGAAAAAA‑TCCCCTTG‑CCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCCTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCaa                                                                           >  1:342476/1‑239 (MQ=255)                                                                                                                                          gTGAAAACGGATTAAACGGCAGGTTA‑AACCGAGTATCTTTGTGAATATTTTTTCAC‑GTTAGTATCAAGTGGC‑TGTGAGGACGCGAAAAAA‑TCCCCTTG‑CCCGCCTGTAACACCGCGCTGAATCAgg                                                                                                                                                                                           >  1:758278/1‑126 (MQ=255)                                                                                                                                                     ttAAACCGTAGGGTA‑AACCGTGTATTTTGGTTAATTTTTTTTCCCGGTGGGTAGCAAGTGGC‑TGT‑AGGCAGCCAAAAAACCCCCCTTG‑CCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGAc                                                                 <  2:486522/238‑1 (MQ=255)                                                                                                                                                                          aGTATCTTTGTGAATATTTTTTCAC‑GTTAGTATCAAGTGGC‑TGTGAGGACGCGAAAAAA‑TCCCCTTG‑CCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTccacca                                                                                                                  >  2:450808/1‑168 (MQ=255)                                                                                                                                                                          aGTATCTTTGTGAATATTTTTTCAC‑GTTAGTATCAAGTGGC‑TGTGAGGACGCGAAAAAA‑TCCCCTTG‑CCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTccacca                                                                                                                  <  1:450808/168‑1 (MQ=255)                                                                                                                                                                          aGTATCTTTGTGAATATTTTTTCAC‑GTTAGTATCAAGTGGC‑TGTGAGGACGCGAAAAAA‑TCCCCTTG‑CCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGt                                           >  1:744903/1‑239 (MQ=255)                                                                                                                                                                           gTATCTTTGTGAATATTTTTTCAC‑GTTAGTATCAAGTGGC‑TGTGAGGACGCGAAAAAA‑TCCCCTTG‑CCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTGTCAATGTTTCTGACCAGCCCCGGATCGAGCTTGCGtt                                          >  2:931160/1‑238 (MQ=255)                                                                                                                                                                            tATTTTGGTAAATTTTTTTCCAC‑GTAAGTACCAAGGGGC‑TGTGAGACCGCAAAAAAC‑CCCACTGC‑CCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGt                                           <  2:584967/237‑1 (MQ=255)                                                                                                                                                                                    taaaTTTTTTTCCAC‑GTTAGTATCAGGGGGC‑TGGGGGGACGCAAAAAAA‑CCCCCTTG‑CCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACAt                                 <  1:931160/237‑1 (MQ=255)                                                                                                                                                                                      aaTATTTTTTCAC‑GTTAGTATCAAGTGGC‑TGTGAGGACGCGAAAAAA‑TCCCCTTG‑CCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACAt                                 >  1:584967/1‑237 (MQ=255)                                                                                                                                                                                       ttatTTTTGCAT‑GTTAGTATCAGGTGGC‑TGTGAGCACGCG‑AAAAC‑CCCCCTTGCCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGGGCTTGAGTAACTCAACATTaa                              <  2:645130/238‑1 (MQ=255)                                                                                                                                                                                                    tttattaCCAGGTGGC‑TGGGAGG‑CCCGAAAAAACCCCCCTTG‑CCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGc                                                                                                               <  2:750888/145‑1 (MQ=255)                                                                                                                                                                                                      tAGTATCAAGTGGC‑TGTGAGGACGCGAAAAAA‑TCCCCTTG‑CCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGCCCAGCCCCGGATGGAGCTTGAGTAACTCACAATTACGTTCTTTGATCCa                 >  2:626426/1‑238 (MQ=255)                                                                                                                                                                                                        ttataAAGTGGC‑TGTGAGGACGCAAAAAAC‑CCCCCTTC‑CCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGAc                                                                 <  2:63264/187‑1 (MQ=255)                                                                                                                                                                                                        ttaacatgTGGT‑TGGGAGGCCGGAAAAAAC‑CCCCCTTG‑CCCGCCTGTAACCCCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTtct                                                                              <  2:905587/168‑1 (MQ=255)                                                                                                                                                                                                         tagcaggTGGC‑GGTGAGGCCGCGAAAAAC‑CCCCCTTG‑CCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGa                                                       <  2:589927/191‑1 (MQ=255)                                                                                                                                                                                                           tCAGGTGGC‑TGTGAGAACGCAAAAAAC‑CCCCCTTC‑CCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAg                                                                                                                               <  2:537298/123‑1 (MQ=255)                                                                                                                                                                                                              ggggcc‑TGTGAGGCCGCAAAAAAC‑CCCCCTGG‑CCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCg                                                                                   <  2:826780/159‑1 (MQ=255)                                                                                                                                                                                                              aGTGGC‑TGTGAGGACGCGAAAAAA‑TCCCCTTG‑ccc                                                                                                                                                                                                                    >  1:188355/1‑35 (MQ=255)                                                                                                                                                                                                              aGTGGC‑TGTGAGGACGCGAAAAAA‑TCCCCTTG‑ccc                                                                                                                                                                                                                    <  2:188355/35‑1 (MQ=255)                                                                                                                                                                                                              aGTGGC‑TGTGAGGACGCGAAAAAA‑TCCCCTTG‑CCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTcc                                                                                                                                    >  2:592175/1‑115 (MQ=255)                                                                                                                                                                                                              aGTGGC‑TGTGAGGACGCGAAAAAA‑TCCCCTTG‑CCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTcc                                                                                                                                    <  1:592175/115‑1 (MQ=255)                                                                                                                                                                                                              aGTGGC‑TGTGAGGACGCGAAAAAA‑TCCCCTTG‑CCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACcag                                                                                                     >  1:750888/1‑146 (MQ=255)                                                                                                                                                                                                              aGTGGC‑TGTGAGGACGCGAAAAAA‑TCCCCTTG‑CCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACt        >  2:310029/1‑239 (MQ=255)                                                                                                                                                                                                              aGTGGC‑TGTGAGGACGCGAAAAAA‑TCCCCTTG‑CCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACt        >  1:352259/1‑239 (MQ=255)                                                                                                                                                                                                                  gC‑TGTGAGGACGCGAAAAAA‑TCCCCTTG‑CCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCgg                                                                                                                                                            >  2:668090/1‑87 (MQ=255)                                                                                                                                                                                                                  gC‑TGTGAGGACGCGAAAAAA‑TCCCCTTG‑CCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCgg                                                                                                                                                            <  1:668090/87‑1 (MQ=255)                                                                                                                                                                                                                  gC‑TGTGAGGACGCGAAAAAA‑TCCCCTTG‑CCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCt                                                                    >  1:905587/1‑175 (MQ=255)                                                                                                                                                                                                                  gC‑TGTGAGGACGCGAAAAAA‑TCCCCTTG‑CCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGa                                                       >  1:63264/1‑188 (MQ=255)                                                                                                                                                                                                                   c‑TGTGAGGACGCGAAAAAA‑TCCCCTTG‑CCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGa                                             >  1:589927/1‑197 (MQ=255)                                                                                                                                                                                                                   c‑TGTGAGGACGCGAAAAAA‑TCCCCTTG‑CCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGATGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCa   >  1:529497/1‑239 (MQ=255)                                                                                                                                                                                                                     tgcgAGGACACGAAAAAA‑TCCCCTTG‑CCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGCTCCCGGAGCAAGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCCCAATGTTTCTGACCAGCCCCGGATCGAGCATGAGTACATCTGCCATAGGGTCCTTGATACAACAGGTACtaccaag  >  2:804248/4‑233 (MQ=255)                                                                                                                                                                                                                      gtgAGGACGCGAAAAAA‑TCCCCTTG‑CCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTcca                                                                                                                     >  1:537298/1‑123 (MQ=255)                                                                                                                                                                                                                         aGGACGCGAAAAAA‑TCCCCTTG‑CCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATg                                                                         >  1:826780/1‑164 (MQ=255)                                                                                                                                                                                                                                        |                                                                                                                                                                                                                              GTCCAGAAAGGCTTAATGCATCAGCAGCATTCCTCGCCTGTTTGCAACCGCGTAAGGTCATTCATCTGCATATGTGCTGCAACAATCTGGCGAGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTA‑AACCGAGTATCTTTGTGAATATTTTTTCAC‑GTTAGTATCAAGTGGC‑TGTGAGGACGCGAAAAAA‑TCCCCTTG‑CCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAG  >  NC_000913/3401232‑3401679

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