Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,556,205	0	A	C	26.9%	147.4 / 15.0	67	S215R (AGT→AGG) ‡	rtcA	RNA 3'‑terminal phosphate cyclase
Reads supporting (aligned to +/- strand):  ref base A (44/5);  new base C (0/18);  total (44/23)
Fisher's exact test for biased strand distribution p-value = 3.59e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGAC‑GGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCC  >  NC_000913/3556066‑3556343                                                                                                                                             |                                                                                                                                            ccaccTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGAC‑GGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGaaaaa                                                                                                                                             >  1:377347/1‑140 (MQ=255)   accTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGAC‑GGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGaaaaa                                                                                                                                             >  2:2507090/1‑138 (MQ=255)   accTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGAC‑GGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGaaaaa                                                                                                                                             <  1:2507090/138‑1 (MQ=255)   accTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGAC‑GGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑Ct                                                                                                                                          >  1:958390/1‑140 (MQ=255)    ccccGGCCCTGCCGCGTTTTTCCCCGACAAAAAAAAAGGGTTGGGTGTTTTTTTCCCTTTCGATTTAAAGCGAGAG‑GTTATTCCCCGGCCCCTGGTCGCGCGGCAGTTTATAAATTTTCTGTTCATG‑AGGGAAAAA‑cccc                                                                                                                                        <  2:958390/137‑1 (MQ=255)      tCGGCACTGACGCGTTTTCCCCCGACGACAAAAAAGCGTTCGGTGATATTTTCCCTTTCGACTTCAAGCGAGAC‑GGTATTCCCCGGCCCCTGGCCGCGCGGCAGTTTATGAATATTCTGTTCATGCA‑GGAAAAA‑cccccc                                                                                                                                      <  2:1763250/140‑1 (MQ=255)      ccGGCATTGCCGCGTTTTCCACCTGCGACAAAAAAGCTTTCGGAGATATTTTCATTTTCGCTTCCAAGCGAGAC‑GATATTACCCGCCCCCGGGCCGCGGGCCATGTTATGAATATTCTGTTCATGC‑GGGAAAAA‑GCcccc                                                                                                                                      <  2:2242304/139‑1 (MQ=255)             tGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGAC‑GGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAgt                                                                                                                                >  1:153684/1‑140 (MQ=255)               aCGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGAC‑GGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAgtgt                                                                                                                              >  1:2693573/1‑140 (MQ=255)                    cTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGAC‑GGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGa                                                                                                                         >  1:916262/1‑140 (MQ=255)                      tttCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGAC‑GGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGAtt                                                                                                                       >  1:2301508/1‑140 (MQ=255)                         cACCGCCGAAAAAAAAGCTTTCGGTGATTTTTTCCCTTTCGACTTCAAGCGAGAG‑GGTATTCCCCGGCCCCTGGTCGCGCGGCGGGTTATAAATTTTTTGTTCATGAAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCa                                                                                                                    <  2:2099553/140‑1 (MQ=255)                           ccacccacAAAAAAGCTTTGGGTGATTTTTTCACTTTCCATTTAAAGCGAGAC‑GGTTTTCCCCGGCCCCTGGCCGCGGGGCGGGTTATGAATATTCTGTTCAGGCAGGAAAAAA‑CCCCCCGCCAGTGTAGCGATTTCACg                                                                                                                  <  2:2271577/134‑1 (MQ=255)                           ccGACGACAAAAAGGCGTTCGGTGATTTTTTCCCTTTCGACTTCAAGCGAGAC‑GGTATTCCCCGGCCCTTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGAAAAAC‑CTCCCCGCCAGTGTAGCGATTTCACg                                                                                                                  <  2:2301508/140‑1 (MQ=255)                             gccgacAAAAAGGCGTGCGGTTATTTTTTCACTTTCGATTTCAAGCGAGAC‑GTTATTCCCCGGCCCCGGGTCGCGGGGCAGGTTATAAATATTTTTTTCAGGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCt                                                                                                                <  2:2693573/138‑1 (MQ=255)                             cccgaaAAAAACGGGTTGGGAAATTTTTTCCCTTCCGTCTTCAGCCGGGAC‑GGTATCCCCCGGCCCGGGGTCGCGTTGCGGTTTATGAATTTTCTGTTAATGCAGGAAAAAG‑CCCCCCGCCAGTGTAGCGATTTCACGCt                                                                                                                <  2:659222/138‑1 (MQ=255)                                 acAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGAC‑GGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGc                                                                                                            >  1:670424/1‑140 (MQ=255)                                        ggCGTTCGGTGATTTTTTCTCTTTCGACTTCAAGCGAGAG‑GGTATTCCCCGGCCCCTGGTCGCGCGGCAGGTTATAAATTTTCTGTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATgg                                                                                                     <  2:2377182/139‑1 (MQ=255)                                          ggTTCGGTATTTTTTTCCTTTTCGACTCCAAGCGAGCG‑GTTTTTCCCCGCCCCCGGGTCGCGGGGCAGTTAATAAATATTCTTTTCATG‑AGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgcg                                                                                                   <  2:670424/139‑1 (MQ=255)                                            ttCGGTGATATTTTCACTTTCGACTTCAAGCGAGAC‑GGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGcgcg                                                                                                 >  1:2327321/1‑140 (MQ=255)                                             tgggTGATTTTTTCTCTTCCGATTTAAGGCGAGGC‑GGTTTTCCCCGGCCCCGGGGCGGGCGGCGGTTTATAAATTTTTTGTTCAGGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCgg                                                                                                <  2:2749751/138‑1 (MQ=255)                                             tCGGTGATATTTTCACTTTCGACTTCAAGCGAGAC‑GGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCgg                                                                                                >  1:1830821/1‑140 (MQ=255)                                             tCGATGATATTTTCACTTTCGACTTCAAGCGAGAC‑GGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCgg                                                                                                >  2:291794/1‑140 (MQ=255)                                                 ttattttttCCTTTTCGATTTCAAGCGAGAC‑GGTATTCCCCGGCCCCTGGCCGCGCGGCAGTTTATAAATTTTCTGTTCAGGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGcaca                                                                                            <  2:2327321/135‑1 (MQ=255)                                                    ttttttCACTTTCGCCTTAAAGCGAGAC‑GTTATTCCCCGGCCCCTGGTCGCGCGGCAGTTTATGAATTTTCTGTTCAGGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCa                                                                                         <  2:849292/138‑1 (MQ=255)                                                    tatTTTCATTTGCGCCTTCAAGCGAGGC‑GGTTTTCCCGGCCCCCGGGTCGCGGGGCGGGTTATAATTTTTTTGTTAATGAAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCa                                                                                         <  2:2658583/140‑1 (MQ=255)                                                         tCACTTTCGACTTCAAGCGAGAC‑GGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTaa                                                                                    >  2:1653111/1‑140 (MQ=255)                                                           aCTTTCGACTTCAAGCGAGAC‑GGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATa                                                                                  >  2:2735991/1‑140 (MQ=255)                                                             tttCGACTTCAAGCGAGAC‑GGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGa                                                                                >  1:849292/1‑140 (MQ=255)                                                               tCGACTTCAAGCGAGAC‑GGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAAc                                                                              >  1:1378737/1‑140 (MQ=255)                                                                      cAAGGGAGAC‑GGTATTCCCCGGCCCCGGGCCGCGCGGCAGTTTATAAATTTTCTTTTCATGAAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCa                                                                       <  2:2318179/140‑1 (MQ=255)                                                                      aaggggaggC‑GTTATCACCCGGCCCCTGTCCGCGCGTCTGTTCATAAATATTCTCTTTATGCAGGAAAAGC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCa                                                                       <  2:1378737/135‑1 (MQ=255)                                                                      aaggCGAGAC‑GTTATTCCCCGGCCCCGGGCCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGAAAAAA‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCa                                                                       <  2:436691/137‑1 (MQ=255)                                                                           gagcC‑GGTATCCCCCGGCCCCGGGCCGCGCGGCAGGTTATAAATTTTCTGTTCATGCGGGGAAAAC‑CCCCCCGCCAGTGTAGCGATCTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCAt                                                                  <  2:141499/140‑1 (MQ=255)                                                                           gagaCGGGTATTCCCCGGCCCCTGGCCGCGCGGCAGTTTATAAATTTTCTGTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCa                                                                   <  2:2650048/140‑1 (MQ=255)                                                                           gagaC‑GTTATTCCCCGGCCCCTGGCCGGGCGGCAGGTTATGAATTTTCTGTTCATGAAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCAt                                                                  <  2:1830821/140‑1 (MQ=255)                                                                           gagaC‑GTTATTCCCCGCCCCCGGGCCGCGCGGCAGGTTATGAATTTTTTGTTCATGCAGGAAAAAA‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCAt                                                                  <  2:153684/140‑1 (MQ=255)                                                                            agacaggTATTCCCCGGCCCCTGGTCGCGCGGCAGGTTATGAATTTTCTGTTCAGGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCAt                                                                  <  2:1544423/135‑1 (MQ=255)                                                                              aC‑GGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTg                                                               <  1:355098/140‑1 (MQ=255)                                                                               aagttATTCCCCGCCCCCGGGTCGCGCGGCAGTTTATGAATTTTCTTTTCATGCATGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTg                                                               <  2:804270/136‑1 (MQ=255)                                                                                   tATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGcaca                                                           >  1:944806/1‑140 (MQ=255)                                                                                   tATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGcaca                                                           >  1:526725/1‑140 (MQ=255)                                                                                           ggCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTccc                                                   >  1:2104292/1‑140 (MQ=255)                                                                                             ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCATGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTccctg                                                 >  1:804270/1‑138 (MQ=255)                                                                                             ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATgctgtctc                                                 >  1:2650048/1‑133 (MQ=255)                                                                                             ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGtctgtct                                                 >  1:1544423/1‑134 (MQ=255)                                                                                             ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc                                                 >  1:1651664/1‑140 (MQ=255)                                                                                             ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc                                                 >  1:2588643/1‑140 (MQ=255)                                                                                             ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc                                                 >  2:150383/1‑140 (MQ=255)                                                                                             ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc                                                 >  1:2518035/1‑140 (MQ=255)                                                                                             ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc                                                 >  1:643739/1‑140 (MQ=255)                                                                                             ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc                                                 >  2:2472206/1‑140 (MQ=255)                                                                                                  ggccgcgggGCGGTTTATAAATTTTCTTTTCAGGCAGGGAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCg                                            <  2:2508924/137‑1 (MQ=255)                                                                                                    gccgggggCAGTTTATGAATTTTCTTTTCAGGAAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGcc                                          <  2:377347/134‑1 (MQ=255)                                                                                                       cgcgGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAg                                       >  2:751544/1‑140 (MQ=255)                                                                                                       cgcgGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAg                                       >  2:1366181/1‑140 (MQ=255)                                                                                                        gcgGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGt                                      >  1:2508924/1‑140 (MQ=255)                                                                                                            cAGTTTATAAATATTCTGTTCAGGAAGGGAAAAA‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCa                                  <  2:1971972/140‑1 (MQ=255)                                                                                                            cAGGTTATAAATATTCTGTTCATGCAGGAAAAAC‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCa                                  <  1:5776/140‑1 (MQ=255)                                                                                                             agtttATGAATATTCTGTTCATGAAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCaa                                 <  1:48701/137‑1 (MQ=255)                                                                                                                         ttCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCaa                                 <  1:652505/128‑1 (MQ=255)                                                                                                                         ttCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCaa                                 >  2:652505/1‑128 (MQ=255)                                                                                                                         ttCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGa                     >  1:1579745/1‑140 (MQ=255)                                                                                                                         ttCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGa                     >  2:2781650/1‑140 (MQ=255)                                                                                                                         ttCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGa                     >  2:1111221/1‑140 (MQ=255)                                                                                                                         ttCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGa                     >  1:1676190/1‑140 (MQ=255)                                                                                                                           ttGTTCAGGGGGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATg                   <  2:1676190/139‑1 (MQ=255)                                                                                                                              ttCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAATGATGCCa                >  1:1978713/1‑140 (MQ=255)                                                                                                                              ttCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCa                >  2:2137334/1‑140 (MQ=255)                                                                                                                              ttCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCa                >  2:1704082/1‑140 (MQ=255)                                                                                                                              ttCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCa                >  1:424127/1‑140 (MQ=255)                                                                                                                              ttCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCa                >  2:888258/1‑140 (MQ=255)                                                                                                                              ttCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCa                >  1:2220237/1‑140 (MQ=255)                                                                                                                               tCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCAc               >  1:418237/1‑140 (MQ=255)                                                                                                                                  tGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACggg            >  2:2032870/1‑140 (MQ=255)                                                                                                                                   gCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGc           >  2:1779709/1‑140 (MQ=255)                                                                                                                                   gCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGc           >  1:978632/1‑140 (MQ=255)                                                                                                                                      gggAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCgag        >  1:593281/1‑140 (MQ=255)                                                                                                                                            aa‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTcc  <  2:2104292/140‑1 (MQ=255)                                                                                                                                             |                                                                                                                                            CCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGAC‑GGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCC  >  NC_000913/3556066‑3556343

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