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→TTTynfEputative 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.

TAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTAATCATCCCGATCGCTT  >  NC_000913/1657926‑1658403
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 
tAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTTCCGTTTAcc                                                                                                                                                                                                                                         >  1:258345/1‑251 (MQ=255)
 aaTCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTAcc                                                                                                                                                                                                                                         <  1:441157/250‑1 (MQ=255)
   tCCCACATTATTACTCATAACCTTCTATTGATATGGCTTAATAATTCTTAACCCACAATGGGTCGACTCCCTCTATTGTTAGCGCGCTAAATATTAAATATATAAACTTTTCTATAACGATAAAGAACAGGGAGAGCGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGCTAGGTTCTCTGGCCCGGGCTGCAGGCGGTTCTTTTTTACCGCTTTAcccg                                                                                                                                                                                                                                       >  2:553905/1‑250 (MQ=255)
            attaCTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTTCCGTTTACCCTGCGCACTgcc                                                                                                                                                                                                                            >  2:102133/1‑250 (MQ=255)
                    tACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCATCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTTTCCTGCGCAATTcagcagcagg                                                                                                                                                                                                                   >  2:198058/1‑250 (MQ=255)
                      cccTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTTCCGTTTACCCCGCGCAAAGCAGCAGCAGCgg                                                                                                                                                                                                                 >  1:610293/1‑251 (MQ=255)
                      cccTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCCTCCCAATGCCGCAGCCGCgg                                                                                                                                                                                                                 >  1:637173/1‑251 (MQ=255)
                             aTTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGTCAGCAGGCGGTTTTTTTTTTCCGTTTACCCTGCGCACAGCagccgcagctggacaaa                                                                                                                                                                                                           >  1:246431/1‑242 (MQ=255)
                                gATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTTCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTacaaaaagg                                                                                                                                                                                                       >  1:433798/1‑246 (MQ=255)
                                                   tAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTg                                                                                                                                                                                    <  1:159545/250‑1 (MQ=255)
                                                    aaCCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGTCGCTGGCTGCAGGCGGTTTTTCTTTTCCGTTTTCCCCGCGCGAAGCCACCGCCCCGGTACAACAGGCCCGCGCAAAAgctgtttggg                                                                                                                                                                                   >  2:220060/1‑242 (MQ=255)
                                                     aCCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTgggg                                                                                                                                                                                 <  1:141274/251‑1 (MQ=255)
                                                        cAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTTTCGTTTACCCTGCGCAATGCAGCCGCAGCCGTACAACAGGCCCGCGCGAAAGAGGGCTggggggg                                                                                                                                                                              >  2:206072/1‑248 (MQ=255)
                                                               ggTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTcc                                                                                                                                                                       <  1:48331/251‑1 (MQ=255)
                                                                  agccTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGTTAAAGAACAGGGAGCGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTTCCGTTTTCCCTTCTCAAAGCAGCAGCCGCGGTTAAACAAGCCCACGAAAAAGTggatcagggggcccggtcccac                                                                                                                                                                     >  1:587273/4‑231 (MQ=255)
                                                                  aGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTc                                                                                                                                                                    <  1:570384/251‑1 (MQ=255)
                                                                     cTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTTCCGTTTACCCTGCGCAAAGCAGCAACAGCGGTACAACAGGCCCGCGAAAAAGAGGGCTGGGGTGCCTGTTCCGTCAAc                                                                                                                                                                 >  2:531542/1‑251 (MQ=255)
                                                                                           aaaTATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACtt                                                                                                                                            <  2:258345/250‑1 (MQ=255)
                                                                                           aaaTATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACt                                                                                                                                             <  1:398749/249‑1 (MQ=255)
                                                                                            aaTATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTTCCGTTTACCCTGCGCAATGGAGCCGCCACGGGAACAAACGCCCCCGAAAAAGTGGGCATGGGGGCCTGTTCCCTCAACTGTGGGAACCGCTGCGCCCTTCg                                                                                                                                          >  1:576124/1‑251 (MQ=255)
                                                                                               atTCAATATATAAACTTTTATATAACGATAAAGAACAGGGATTGAGTTATGTCCAAAAATGAACGAATGGTGGGCGTCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCGGGCGCTGGCTGCAGGCGGTTGTTCTTTGCCTTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGTGTGCCTGTTCCGCCAACTGTGGTAGCCGCTGTGCACTTCGTCt                                                                                                                                       <  2:7953/251‑1 (MQ=255)
                                                                                                    atatatAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCCGCGCAATGCAGCAGCAGCGGGACAACAAGCCCCCGGAAAAGGGGTCTGGGGGGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTAcctg                                                                                                                                  >  1:461104/1‑248 (MQ=255)
                                                                                                      atataAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGt                                                                                                                                 <  1:198058/250‑1 (MQ=255)
                                                                                                            aCTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGGCGCTGTGCACTTCGTCTACATGTTAAAGAt                                                                                                                          <  1:229797/251‑1 (MQ=255)
                                                                                                                   tataACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTTCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGGAAAAGGGGGCTGGGGGGCCTGTTCCGTCAACTGTGGGAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGaa                                                                                                                    >  1:506483/1‑250 (MQ=255)
                                                                                                                    ataACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGCACGAATGGTGGGCATCAGCCGCAGAACACTCGTTACATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTTCCCTTTTCCCTGCGCAATGCCGCAGCAACGGTACAACAGGCCCCCGCAAAAATGGTCTGGGGTGCCTTTTCCGTCAACTGTGGTAGCCCCTTTGCACTTCGTCTACATGTTTAAGATAATGAAgg                                                                                                                  >  1:204397/1‑250 (MQ=255)
                                                                                                                      aaCGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCTGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGGAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTg                                                                                                                 <  2:523455/250‑1 (MQ=255)
                                                                                                                       aCGATTAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTTGGGACAGCAGCAGCAGCGGTGCAACAGGCCCGCGAAAAAGGGGGGTGGGGTGGCTGTTCCGCCAAACGGGGTAGACCCTGTGCACTTCGTCTACATGTTTAAGATAATGAAGTGa                                                                                                                >  1:435777/1‑250 (MQ=255)
                                                                                                                          aTAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTTCCGTTTTCCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCCCGGAAAAATGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGGGCGATTCGTCTACCTGTTAAAGATAATGAAGGGACCTg                                                                                                            >  1:443563/1‑251 (MQ=255)
                                                                                                                                  cAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGaaa                                                                                                    <  1:274664/251‑1 (MQ=255)
                                                                                                                                    gggAGTGAGTTATGTCCAAAACTGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTATGGCGCTGGCTGCAGGCGGGTTTTCTTTGCCGTTTATCCCGCGCCATAGAGCAGCGACGGTGAAAAAGGGCCGCGCAAACGAAGGCTGGGGTGCCCGTTCCGTCAAATGTGGGAGCCGCCGTGCACTTCGTCCACATGTTTAAGATAAAGAAGTgtcccgggggggaac                                                                                                   >  2:75375/1‑236 (MQ=255)
                                                                                                                                        gtgagtTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACa                                                                                              <  1:423495/251‑1 (MQ=255)
                                                                                                                                          gagtTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGTAAACTGACAAt                                                                                            <  2:610965/251‑1 (MQ=255)
                                                                                                                                                      aaaaTGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTTCCCTGCAGCAAGCAGCCACCGCGGTACAACCAGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCCTCAACTGTGGGAGCCGCGGGGCACTTCGCTTACATGTTAAAGATAAAGAAGGGAACTGGGGGGGAACGGACACAACCGGccagg                                                                                  >  2:153030/1‑245 (MQ=255)
                                                                                                                                                           gAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTTCCGTTTACCCTGCGCAATGCAGCAGCAGAGGGACAAAAGGGCCCCGAAAAAGAGGTCTGGGGGGCCTGTTCCGTCCACTGTGGGAGCCGCCGTGCCCCTCGTCTCCATGTTAAAGATAATGAAGGGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAgg                                                                           >  1:85638/1‑249 (MQ=255)
                                                                                                                                                             aCGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTTCCGTTTACCCTGCGCAATGCAGCAGCAGCGGGACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGGGGAAACTGACAATACCGGCAGCGATGAGTACg                                                                        >  1:231859/1‑251 (MQ=255)
                                                                                                                                                                 aTGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCaa                                                                    <  2:506483/251‑1 (MQ=255)
                                                                                                                                                                    gtggGCATCAGCCGCAGAACACGCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCa                                                                 <  2:610293/251‑1 (MQ=255)
                                                                                                                                                                         cATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTTCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGGCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTTAAGATAATGAAGTGACCTGGGGGGAAACTGACAATACCGGCAGCGATGGGGACGGCAAACAACAgg                                                            >  1:31180/1‑251 (MQ=255)
                                                                                                                                                                          aTCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCAt                                                                <  1:448700/246‑1 (MQ=255)
                                                                                                                                                                          aTCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCCCCATGCAGCAGCAGCGGGACAACAGGCCCCCGAAAAAACGGGCTGGGGTGCCCTATCCGTCCACTTTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGGCCTGGGGGGAAACTGACAAAAACGGCAGCGATGAGGACGGGAACCAtagaat                                                           >  2:448700/1‑246 (MQ=255)
                                                                                                                                                                                  cAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCCGCGCAAAGCAGGAGCAACGGGACAACAGGCCCGCGGAAAAGGGGGCTGGGGGGCCTGTTCCGTCAACTGTGGGAGCCGCTGTGCACTTCGTCTACATGTTTAAGATAATGAAGTGACCTGGGTGGAAACTGACAATAACGGCAGCGGTGAGTACGGCCACCCTCAGGTAAGCGCCt                                                   >  1:403722/1‑251 (MQ=255)
                                                                                                                                                                                               aaaTCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTGCCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGGGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTGTGCGCGGTc                                       <  1:424156/250‑1 (MQ=255)
                                                                                                                                                                                                        gcgATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATccg                             <  1:438973/251‑1 (MQ=255)
                                                                                                                                                                                                          gATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTTCCGATTTCCCTGCGCCATTTAGCCGCAGCGGTACAACACGCCCGCGAAAAAGGGGAGTGGGGGGCCTGTTTCGGCAAATTTAGAAGCCGCTGTGCACCTCGTCTACATGTTTAAAATAAAGAAGTGACCTGGGGGGAAAATGACAATACCGGCAGCGAGTGTTACGGCAACCATCCGGTACCCGCCCTTTTTCGCCGTCGctccctccccc                           >  1:155080/1‑248 (MQ=255)
                                                                                                                                                                                                                    ctGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTAt                    <  1:102133/248‑1 (MQ=255)
                                                                                                                                                                                                                    ctGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTAt                    <  1:220060/248‑1 (MQ=255)
                                                                                                                                                                                                                        ggctggctgCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTAatcat             <  2:85638/250‑1 (MQ=255)
                                                                                                                                                                                                                        cgcggtctgCGGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTTTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCTCGGTCGCTCCATCCGCCGGCGTATTAatcat             <  2:548126/245‑1 (MQ=255)
                                                                                                                                                                                                                         gctggctgCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTTTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTAatcatc            <  1:505698/251‑1 (MQ=255)
                                                                                                                                                                                                                                  aGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTAATCATCCCGATCGCt   <  2:461104/251‑1 (MQ=255)
                                                                                                                                                                                                                                   ggCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTAATCATCCCGATCGCtt  >  1:13868/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 
TAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTAATCATCCCGATCGCTT  >  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: