Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,658,1670GT31.4% 79.2 / 15.2 52L33F (TTG→TTT) ynfEputative selenate reductase YnfE
Reads supporting (aligned to +/- strand):  ref base G (10/25);  new base T (16/0);  total (27/25)
Fisher's exact test for biased strand distribution p-value = 1.03e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCG‑TTTACCCTGC‑GCAATGCAGCAGC‑AGCGGTACAACAGGCCCGCGAAAAAGTGGTC‑TGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGA‑TGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTAATCATCCCGATCGCTT  >  NC_000913/1657926‑1658403
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 
tAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTTCCG‑TTTAcc                                                                                                                                                                                                                                         >  1:258345/1‑251 (MQ=255)
 aaTCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCG‑TTTAcc                                                                                                                                                                                                                                         <  1:441157/250‑1 (MQ=255)
   tCCCACATTATTACTCATAACCTTCTATTGATATGGCTTAATAATTCTTAACCCACAATGGGTCGACTCCCTCTATTGTTAGCGCGCTAAATATTAAATATATAAACTTTTCTATAACGATAAAGAACAGGGAGAGCGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGCTAGGTTCTCTGGCCCGGGCTGCAGGCGGTTCTTTTTTACCGCTTTAcccg                                                                                                                                                                                                                                       >  2:553905/1‑250 (MQ=255)
            attaCTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTTCCG‑TTTACCCTGC‑GCACTgcc                                                                                                                                                                                                                            >  2:102133/1‑250 (MQ=255)
                    tACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCATCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCG‑TTTTTCCTGC‑GCAATTcagcagc‑agg                                                                                                                                                                                                                   >  2:198058/1‑250 (MQ=255)
                      cccTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTTCCG‑TTTACCCCGC‑GCAAAGCAGCAGC‑AGCgg                                                                                                                                                                                                                 >  1:610293/1‑251 (MQ=255)
                      cccTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCG‑TTTACCCCTC‑CCAATGCCGCAGC‑CGCgg                                                                                                                                                                                                                 >  1:637173/1‑251 (MQ=255)
                             aTTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGTCAGCAGGCGGTTTTTTTTTTCCG‑TTTACCCTGC‑GCACAGCagccgc‑agctggacaaa                                                                                                                                                                                                           >  1:246431/1‑242 (MQ=255)
                                gATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTTCCG‑TTTACCCTGC‑GCAATGCAGCAGC‑AGCGGTacaaaaagg                                                                                                                                                                                                       >  1:433798/1‑246 (MQ=255)
                                                   tAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCG‑TTTACCCTGC‑GCAATGCAGCAGC‑AGCGGTACAACAGGCCCGCGAAAAAGTGGTC‑Tg                                                                                                                                                                                    <  1:159545/250‑1 (MQ=255)
                                                    aaCCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGTCGCTGGCTGCAGGCGGTTTTTCTTTTCCG‑TTTTCCCCGC‑GCGAAGCCACCGC‑CCCGGTACAACAGGCCCGCGCAAAAgctgtttggg                                                                                                                                                                                   >  2:220060/1‑242 (MQ=255)
                                                     aCCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCG‑TTTACCCTGC‑GCAATGCAGCAGC‑AGCGGTACAACAGGCCCGCGAAAAAGTGGTC‑Tgggg                                                                                                                                                                                 <  1:141274/251‑1 (MQ=255)
                                                        cAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTTTCG‑TTTACCCTGC‑GCAATGCAGCCGC‑AGCCGTACAACAGGCCCGCGCGAAAGAGGGC‑Tggggggg                                                                                                                                                                              >  2:206072/1‑248 (MQ=255)
                                                               ggTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCG‑TTTACCCTGC‑GCAATGCAGCAGC‑AGCGGTACAACAGGCCCGCGAAAAAGTGGTC‑TGGGGTGCCTGTTcc                                                                                                                                                                       <  1:48331/251‑1 (MQ=255)
                                                                  agccTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGTTAAAGAACAGGGAGCGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTTCCG‑TTTTCCCTTC‑TCAAAGCAGCAGC‑CGCGGTTAAACAAGCCCACGAAAAAGTggatcagggggcccggtcccac                                                                                                                                                                     >  1:587273/4‑231 (MQ=255)
                                                                  aGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCG‑TTTACCCTGC‑GCAATGCAGCAGC‑AGCGGTACAACAGGCCCGCGAAAAAGTGGTC‑TGGGGTGCCTGTTCCGTc                                                                                                                                                                    <  1:570384/251‑1 (MQ=255)
                                                                     cTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTTCCG‑TTTACCCTGC‑GCAAAGCAGCAAC‑AGCGGTACAACAGGCCCGCGAAAAAGAGGGC‑TGGGGTGCCTGTTCCGTCAAc                                                                                                                                                                 >  2:531542/1‑251 (MQ=255)
                                                                                           aaaTATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCG‑TTTACCCTGC‑GCAATGCAGCAGC‑AGCGGTACAACAGGCCCGCGAAAAAGTGGTC‑TGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACtt                                                                                                                                            <  2:258345/250‑1 (MQ=255)
                                                                                           aaaTATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCG‑TTTACCCTGC‑GCAATGCAGCAGC‑AGCGGTACAACAGGCCCGCGAAAAAGTGGTC‑TGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACt                                                                                                                                             <  1:398749/249‑1 (MQ=255)
                                                                                            aaTATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTTCCG‑TTTACCCTGC‑GCAATGGAGCCGCCACGGGAACAA‑ACGCCCCCGAAAAAGTGGGCATGGGG‑GCCTGTTCCCTCAACTGTGGGAACCGCTGCGCCCTTCg                                                                                                                                          >  1:576124/1‑251 (MQ=255)
                                                                                               atTCAATATATAAACTTTTATATAACGATAAAGAACAGGGATTGAGTTATGTCCAAAAATGAACGAATGGTGGGCGTCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCGGGCGCTGGCTGCAGGCGGTTGTTCTTTGCCT‑TTTACCCTGC‑GCAATGCAGCAGC‑AGCGGTACAACAGGCCCGCGAAAAAGTGGTC‑TGGTGTGCCTGTTCCGCCAACTGTGGTAGCCGCTGTGCACTTCGTCt                                                                                                                                       <  2:7953/251‑1 (MQ=255)
                                                                                                    atatatAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCG‑TTTACCCCGC‑GCAATGCAGCAGC‑AGCGGGACAACAAGCCCCCGGAAAAGGGGTC‑TGGGGGGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTAcctg                                                                                                                                  >  1:461104/1‑248 (MQ=255)
                                                                                                      atataAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCG‑TTTACCCTGC‑GCAATGCAGCAGC‑AGCGGTACAACAGGCCCGCGAAAAAGTGGTC‑TGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGt                                                                                                                                 <  1:198058/250‑1 (MQ=255)
                                                                                                            aCTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCG‑TTTACCCTGC‑GCAATGCAGCAGC‑AGCGGTACAACAGGCCCGCGAAAAAGTGGTC‑TGGGGTGCCTGTTCCGTCAACTGTGGTAGGCGCTGTGCACTTCGTCTACATGTTAAAGAt                                                                                                                          <  1:229797/251‑1 (MQ=255)
                                                                                                                   tataACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCG‑TTTACCCTTC‑GCAATGCAGCAGC‑AGCGGTACAACAGGCCCGCGGAAAAGGGGGC‑TGGGGGGCCTGTTCCGTCAACTGTGGGAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGaa                                                                                                                    >  1:506483/1‑250 (MQ=255)
                                                                                                                    ataACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGCACGAATGGTGGGCATCAGCCGCAGAACACTCGTTACATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTTCCC‑TTTTCCCTGC‑GCAATGCCGCAGC‑AACGGTACAACAGGCCCCCGCAAAAATGGTC‑TGGGGTGCCTTTTCCGTCAACTGTGGTAGCCCCTTTGCACTTCGTCTACATGTTTAAGATAATGAAgg                                                                                                                  >  1:204397/1‑250 (MQ=255)
                                                                                                                      aaCGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCTGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCG‑TTTACCCTGC‑GCAATGCAGCAGC‑AGCGGTACAACAGGCCCGCGAAAAAGTGGTC‑TGGGGTGCCTGTTCCGTCAACTGTGGGAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTg                                                                                                                 <  2:523455/250‑1 (MQ=255)
                                                                                                                       aCGATTAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCG‑TTTACCCTTG‑GGACAGCAGCAGC‑AGCGGTGCAACAGGCCCGCGAAAAAGGGGGG‑TGGGGTGGCTGTTCCGCCAAACGGGGTAGACCCTGTGCACTTCGTCTACATGTTTAAGATAATGAAGTGa                                                                                                                >  1:435777/1‑250 (MQ=255)
                                                                                                                          aTAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTTCCG‑TTTTCCCTGC‑GCAATGCAGCAGC‑AGCGGTACAACAGGCCCCCGGAAAAATGGTC‑TGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGGGCGATTCGTCTACCTGTTAAAGATAATGAAGGGACCTg                                                                                                            >  1:443563/1‑251 (MQ=255)
                                                                                                                                  cAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCG‑TTTACCCTGC‑GCAATGCAGCAGC‑AGCGGTACAACAGGCCCGCGAAAAAGTGGTC‑TGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGaaa                                                                                                    <  1:274664/251‑1 (MQ=255)
                                                                                                                                    gggAGTGAGTTATGTCCAAAACTGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTATGGCGCTGGCTGCAGGCGGGTTTTCTTTGCCG‑TTTATCCCGC‑GCCATAGAGCAGC‑GACGGTGAAAAAGGGCCGCGCAAACGAAGGC‑TGGGGTGCCCGTTCCGTCAAATGTGGGAGCCGCCGTGCACTTCGTCCACATGTTTAAGATAAAGAAGTgtcccgggggggaac                                                                                                   >  2:75375/1‑236 (MQ=255)
                                                                                                                                        gtgagtTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCG‑TTTACCCTGC‑GCAATGCAGCAGC‑AGCGGTACAACAGGCCCGCGAAAAAGTGGTC‑TGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACa                                                                                              <  1:423495/251‑1 (MQ=255)
                                                                                                                                          gagtTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCG‑TTTACCCTGC‑GCAATGCAGCAGC‑AGCGGTACAACAGGCCCGCGAAAAAGTGGTC‑TGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGTAAACTGACAAt                                                                                            <  2:610965/251‑1 (MQ=255)
                                                                                                                                                      aaaaTGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCG‑TTTTCCCTGCAGCAA‑GCAGCCAC‑CGCGGTACAACCAGCCCGCGAAAAAGTGGTC‑TGGGGTGCCTGTTCCCTCAACTGTGGGAGCCGCGGGGCACTTCGCTTACATGTTAAAGATAAAGAAGGGAACTGGGGGGGAACGGACACAACCGGccagg                                                                                  >  2:153030/1‑245 (MQ=255)
                                                                                                                                                           gAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTTCCG‑TTTACCCTGC‑GCAATGCAGCAGC‑AGAGGGACAAAAGGGCCCCGAAAAAGAGGTC‑TGGGGGGCCTGTTCCGTCCACTGTGGGAGCCGCCGTGCCCCTCGTCTCCATGTTAAAGATAATGAAGGGACCTGGGTGGAAACTGACAATACCGGCAGCGA‑TGAgg                                                                           >  1:85638/1‑249 (MQ=255)
                                                                                                                                                             aCGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTTCCG‑TTTACCCTGC‑GCAATGCAGCAGC‑AGCGGGACAACAGGCCCGCGAAAAAGTGGTC‑TGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGGGGAAACTGACAATACCGGCAGCGA‑TGAGTACg                                                                        >  1:231859/1‑251 (MQ=255)
                                                                                                                                                                 aTGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCG‑TTTACCCTGC‑GCAATGCAGCAGC‑AGCGGTACAACAGGCCCGCGAAAAAGTGGTC‑TGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGA‑TGAGTACGGCaa                                                                    <  2:506483/251‑1 (MQ=255)
                                                                                                                                                                    gtggGCATCAGCCGCAGAACACGCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCG‑TTTACCCTGC‑GCAATGCAGCAGC‑AGCGGTACAACAGGCCCGCGAAAAAGTGGTC‑TGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGA‑TGAGTACGGCAACCa                                                                 <  2:610293/251‑1 (MQ=255)
                                                                                                                                                                         cATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTTCCG‑TTTACCCTGC‑GCAATGCAGCAGC‑AGCGGTACAACAGGCCCGCGAAAAAGTGGGC‑TGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTTAAGATAATGAAGTGACCTGGGGGGAAACTGACAATACCGGCAGCGA‑TGGGGACGGCAAACAACAgg                                                            >  1:31180/1‑251 (MQ=255)
                                                                                                                                                                          aTCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCG‑TTTACCCTGC‑GCAATGCAGCAGC‑AGCGGTACAACAGGCCCGCGAAAAAGTGGTC‑TGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGA‑TGAGTACGGCAACCAt                                                                <  1:448700/246‑1 (MQ=255)
                                                                                                                                                                          aTCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCG‑TTTACCCTGC‑CCCATGCAGCAGC‑AGCGGGACAACAGGCCCCCGAAAAAACGGGC‑TGGGGTGCCCTATCCGTCCACTTTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGGCCTGGGGGGAAACTGACAAAAACGGCAGCGA‑TGAGGACGGGAACCAtagaat                                                           >  2:448700/1‑246 (MQ=255)
                                                                                                                                                                                  cAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCG‑TTTACCCCGC‑GCAAAGCAGGAGC‑AACGGGACAACAGGCCCGCGGAAAAGGGGGC‑TGGGGGGCCTGTTCCGTCAACTGTGGGAGCCGCTGTGCACTTCGTCTACATGTTTAAGATAATGAAGTGACCTGGGTGGAAACTGACAATAACGGCAGCGG‑TGAGTACGGCCACCCTCAGGTAAGCGCCt                                                   >  1:403722/1‑251 (MQ=255)
                                                                                                                                                                                               aaaTCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCG‑TTTGCCCTGC‑GCAATGCAGCAGC‑AGCGGTACAACAGGCCCGCGAAAAAGTGGTC‑TGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGGGACCTGGGTGGAAACTGACAATACCGGCAGCGA‑TGAGTACGGCAACCATCAGGTACGCGCCTGTGTGCGCGGTc                                       <  1:424156/250‑1 (MQ=255)
                                                                                                                                                                                                        gcgATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCG‑TTTACCCTGC‑GCAATGCAGCAGC‑AGCGGTACAACAGGCCCGCGAAAAAGTGGTC‑TGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGA‑TGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATccg                             <  1:438973/251‑1 (MQ=255)
                                                                                                                                                                                                          gATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTTCCG‑ATTTCCCTGC‑GCCATTTAGCCGC‑AGCGGTACAACACGCCCGCGAAAAAGGGGAG‑TGGGGGGCCTGTTTCGGCAAATTTAGAAGCCGCTGTGCACCTCGTCTACATGTTTAAAATAAAGAAGTGACCTGGGGGGAAAATGACAATACCGGCAGCGAGTG‑TTACGGCAACCATCCGGTACCCGCCCTTTTTCGCCGTCGctccctccccc                           >  1:155080/1‑248 (MQ=255)
                                                                                                                                                                                                                    ctGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCG‑TTTACCCTGC‑GCAATGCAGCAGC‑AGCGGTACAACAGGCCCGCGAAAAAGTGGTC‑TGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGA‑TGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTAt                    <  1:102133/248‑1 (MQ=255)
                                                                                                                                                                                                                    ctGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCG‑TTTACCCTGC‑GCAATGCAGCAGC‑AGCGGTACAACAGGCCCGCGAAAAAGTGGTC‑TGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGA‑TGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTAt                    <  1:220060/248‑1 (MQ=255)
                                                                                                                                                                                                                        ggctggctgCAGGCGGTTTTTCTTTGCCG‑TTTACCCTGC‑GCAATGCAGCAGC‑AGCGGTACAACAGGCCCGCGAAAAAGTGGTC‑TGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGA‑TGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTAatcat             <  2:85638/250‑1 (MQ=255)
                                                                                                                                                                                                                        cgcggtctgCGGGCGGTTTTTCTTTGCCG‑TTTACCCTGC‑GCAATGCAGCAGC‑AGCGGTACAACAGGCCCGCGAAAAAGTGGTT‑TGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGA‑TGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCTCGGTCGCTCCATCCGCCGGCGTATTAatcat             <  2:548126/245‑1 (MQ=255)
                                                                                                                                                                                                                         gctggctgCAGGCGGTTTTTCTTTGCCG‑TTTACCCTGC‑GCAATGCAGCAGC‑AGCGGTACAACAGGCCCGCGAAAAAGTGGTC‑TGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTTTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGA‑TGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTAatcatc            <  1:505698/251‑1 (MQ=255)
                                                                                                                                                                                                                                  aGGCGGTTTTTCTTTGCCG‑TTTACCCTGC‑GCAATGCAGCAGC‑AGCGGTACAACAGGCCCGCGAAAAAGTGGTC‑TGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGA‑TGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTAATCATCCCGATCGCt   <  2:461104/251‑1 (MQ=255)
                                                                                                                                                                                                                                   ggCGGTTTTTCTTTGCCG‑TTTACCCTGC‑GCAATGCAGCAGC‑AGCGGTACAACAGGCCCGCGAAAAAGTGGTC‑TGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGA‑TGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTAATCATCCCGATCGCtt  >  1:13868/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 
TAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCG‑TTTACCCTGC‑GCAATGCAGCAGC‑AGCGGTACAACAGGCCCGCGAAAAAGTGGTC‑TGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGA‑TGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTAATCATCCCGATCGCTT  >  NC_000913/1657926‑1658403

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