Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,328,4650TG33.9% 116.2 / 11.5 65T216T (ACA→ACC) yciKputative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (12/31);  new base G (22/0);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 4.68e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CAAATGCCGCGGTGGTTTTGCCTTTTCCATTGCCGGTAAAGACGATGATAATACCGCGTTCATCCTGGGCCTGGGCCACACGAGCATCTACTTTTTCTTTCACTCGCTGCTGACGCTGTTGGTAGCGTTCATCACTCATTGGGAAATTCCTGGTTTACGGCCCGGTTGGGCGTCAAAGGTCATGCCGGTTTTACGGCGGCTGTCATCGCCCATCAGCCAGAGGTAGAGCGGCATGATATCAGCGGGTGTTTT‑AAGTTTCTGTGGATCTTCGGTCGGGAAGGCGCTGGCACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCAT‑CGGCCAGTACCTGCATCATCCCTTCGGTGG‑CAAATTTCGACGCTGC‑ATATGCACCCCAGTTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGT  >  NC_000913/1328219‑1328711
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                          
cAAATGCCGCGGTGGTTTTGCCTTTTCCATTGCCGGTAAAGACGATGATAATACCGCGTTCATCCTGGGCCTGGGCCACACGAGCATCTACTTTTTCTTTCACTCGCTGCTGACGCTGTTGGTAGCGTTCATCACTCATTGGGAAATTCCTGGTTTACGGCCCGGTTGGGCGTCAAAGGTCATGCCGGTTTTACGGCGGCTGTCATCGCCCATCAGCCAGAGGTAGAGCGGCATGATATCAGCGGGTGttt                                                                                                                                                                                                                                                        <  1:293751/251‑1 (MQ=255)
   aTGCCGCGGTGGTTTTGCCTTTTCCATTGCCGGTAAAGACGATGATAATACCGCGTTCATCCTGGGCCTGGGCCACACGAGCATCTACTTTTTCTTTCACTCGCTGCTGACGCTGTTGGTAGCGTTCATCACTCATTGGGAAATTCCTGGTTTACGGCCCGGTTGGGCGTCAAAGGTCATGCCGGTTTTACGGCGGCTGTCATAGCCCATCAGCCAGAGGTAGAGCGGCATGATATCAGAGGGGGTTTT‑a                                                                                                                                                                                                                                                     >  1:293606/1‑250 (MQ=255)
      ccGCGGTGGTTTTGCCTTTTCCATTGCCGGTAACGACGATGATAATACCGCGTTCATCCTGGGCCTGGGCCACACGAGCATCTACTTTTTCTTTCACTCGCTGCTGACGCTGTTGGTAGCGTTCATCACTCATTGGGAAATTCCTGGTTTACGGCCCGGTTGGGCGTCAAAGGTCATGCCGGTTTTACGGCGGCTTTCATCGCCCCTCAGCCAGAGGTAGAGCGGCCTGATATCAGCGGGGGTTTT‑AAGtt                                                                                                                                                                                                                                                 >  2:255737/1‑251 (MQ=255)
         cgGTGGTTTTGCCTTTTCCATTGCCGGTAAAGACGATGATAATACCGCGTTCATCCTGGGCCTGGGCCACACGAGCATCTACTTTTTCTTTCACTCGCTGCTGACGCTGTTGGTAGCGTTCATCACTCATTGGGAAATTCCTGGTTTACGGCCCGGTTGGGCGTCAAAGGTCATGCCGGTTTTACGGCGGCTGTCATCGCCCATCAGCCAGAGGTAGAGCGGCATGATATCAGCGGGTGTTTT‑AAGTTTCt                                                                                                                                                                                                                                              <  1:480290/251‑1 (MQ=255)
            tggtTTTGCCTTTTCCATTGCCGGTAAAGACGATGATAATACCGCGTTCATCCTGGGCCTGGGCCACACGAGCATCTACTTTTTCTTTCACTCGCTGCTGACGCTGTTGGTAGCGTTCATCACTCATTGGGAAATTCCTGGTTTACGGCCCGGTTGGGCGTCAAAGGTCATGCCGGTTTTACGGCGGCTGTCATCGCCCATCAGCCAGAGGGAGAGCGGCATGATATCAGCGGGGGTTTT‑AAGTTTTtgtg                                                                                                                                                                                                                                           >  1:94412/1‑251 (MQ=255)
            tggtTTTGCCTTTTCCATTGCCGGTAAAGACGATGATAATACCGCGTTCATCCTGGGCCTGGGCCACACGAGCATCTACTTTTTCTTTCACTCGCTGCTGACGCTGTTGGTAGCGTTCATCACTCATTGGGAAATTCCTGGTTTACGGCCCGGTTGGGCGTCAAAGGTCATGCCGGTTTTACGGCGGCCGTCATCTCCCCTCCGCCCGAGGGAGAGCGGCATGATATCAGCGGGGGtttttaatttttttg                                                                                                                                                                                                                                            >  1:219892/1‑240 (MQ=255)
             ggtTTTGCCTTTTCCATTGCCGGTAAAGACGATGATAATACCGCGTTCATCCTGGGCCTGGGCCACACGAGCATCTACTTTTTCTTTCACTCGCTGCTGACGCTGTTGGTAGCGTTCATCACTCATTGGGAAATTCCTGGTTTACGGCCCGGTTGGGCGTCAAAGGTCATGCCGGTTTTACGGCGGCTGTCATCGCCCATCAGCCAGAGGTAGAGCGGCATGATATCAGCGGGTGTTTT‑AAGTTTCtgtg                                                                                                                                                                                                                                           <  1:304741/250‑1 (MQ=255)
               ttttGCCTTTTCCATTGCCGGTAAAGACGATGATAATACCGCGTTCATCCTGGGCCTGGGCCACACGAGCATCTACTTTTTCTTTCACTCGCTGCTGACGCTGTTGGTAGCGTTCATCACTCATTGGGAAATTCCTGGTTTACGGCCCGGTTGGGCGTCAAAGGTCATGCCGGTTTTACGGCGGCTGTCATCGCCCCTCAGCCAGAGGGAGATCGGCATGCTATCCACGGGGGGTTT‑TAGTTTCtgggggt                                                                                                                                                                                                                                        >  1:110007/1‑246 (MQ=255)
                    ccTTTTCCATTGCCGGTAAAGGCGATGATAATACCGCGTTCATCCTGGGCCTGGGCCACACGAGCATCTACTTTTTCTTTCACTCGCTGCTGACGCTGTTGGTAGCGTTCATCACTCATTGGGAAATTCCTGGTTTACGGCCCGGTTGGGCGTCAAAGGTCATGCCGGTTTTACGGCGGGTGTCATCGCCCCTCCGCCCGAGGGAGAGGGGGCTGTTATTAGCGGGGGTGTT‑TAGTTTCttggggactttt                                                                                                                                                                                                                                   >  2:452354/1‑240 (MQ=255)
                              tGCCGGTAAAGACGATGATAATACCGCGTTCATCCTGGGCCTGGGCCACACGAGCATCTACTTTTTCTTTCACTCGCTGCTGACGCTGTTGGTAGCGTTCATCACTCATTGGGAAATTCCTGGTTTACGGCCCGGTTGGGCGTCAAAGGTCATGCCGGTTTTACGGCGGCTGTCATCGCCCATCAGCCACGGGTAGAGCGGCCTGATATCAGCGGGGGGTTT‑AAGTTTTTGTGGGTCTTCGGTCGGGAAgg                                                                                                                                                                                                                         >  1:392994/1‑251 (MQ=255)
                               gCCGGTAAAGACGATGATAATACCGCGTTCATCCTGGGCCTGGGCCACACGAGCATCTACTTTTTCTTTCACTCGCTGCTGACGCTGTTGGTAGCGTTCATCACTCATTGGGAAATTCCTGGTTTACGGCCCGGTTGGGCGTCAAAGGTCATGCCGGTTTTACGGCGGCTGTCATCGCCCATCAGCCAGAGGTAGAGCGGCATGATATCAGCGGGTGTTTT‑AAGTTTCTGTGGATCTTCGGTCGGGAAgg                                                                                                                                                                                                                         <  2:127170/250‑1 (MQ=255)
                                  ggTAAAGACGATGATAATACCGCGTTCATCCTGGGCCTGGGCCACACGAGCATCTACTTTTTCTTTCACTCGCTGCTGACGCTGTTGGTAGCGTTCATCACTCATTGGGAAATTCCTGGTTTACGGCCCGGTTGGGCGTCAAAGGTCATGCCGGTTTTACGGCGGCTGTCATCGCCCATCAGCCAGAGGTAGAGCGGCATGATATCAGCGGGTGTTTT‑AAGTTTCTGTGGATCTTCGGTCGGGAAGGCGCt                                                                                                                                                                                                                     <  1:511104/251‑1 (MQ=255)
                                           gatgatAATACCGCGTTCATCCTGGGCCTGGGCCACACGAGCATCTACTTTTTCTTTCACTCGCTGCTGACGCTGTTGGTAGCGTTCATCACTCATTGGGAAATTCCTGGTTTACGGCCCGGTTGGGCGTCAAAGGTCATGCCGGTTTTACGGCGGCTGTCATCGCCCATCAGCCAGAGGTAGAGCGGCATGATATCAGCGGGTGTTTT‑AAGTTTCTGTGGATCTTCGGTCGGGAAGGCGCTGGCACGCAt                                                                                                                                                                                                            <  1:569699/251‑1 (MQ=255)
                                              gatAATACCGCGTTCATCCTGGGCCTGGGCCACACGAGCATCTACTTTTTCTTTCACTCGCTGCTGACGCTGTTGGTAGCGTTCATCACTCATTGGGAAATTCCTGGTTTACGGCCCGGTTGGGCGTCAAAGGTCATGCCGGTTTTACGGCGGCTGTCATCGCCCATCAGCCAGAGGTAGAGCGGCATGATATCAGCGGGTGTTTT‑AAGTTTCTGTGGATCTTCGGTCGGGAAGGCGCTGGCACGCATTGc                                                                                                                                                                                                         <  1:98400/251‑1 (MQ=255)
                                               ataataCCGCGTTCATCCTGGGCCTGGGCCACACGAGCATCTACTTTTTCTTTCACTCGCTGCTGACGCTGTTGGTAGCGTTCATCACTCATTGGGAAATTCCTGGTTTACGGCCCGGTTGGGCGTCAAAGGTCATGCCGGTTTTACGGCGGCTGTCATCGCCCCTCAGCCCGAGGTTGAGCGGCATGATATCAGCGGGTGGTTT‑AAGGTTCTGGGGGTCTTTGGTCGGGGGGGCGCCGGCCCGCCTTGCg                                                                                                                                                                                                        >  2:53274/1‑251 (MQ=255)
                                                  ataCCGCGTTCATCCTGGGCCTGGGCCACACGAGCATCTACTTTTTCTTTCACTCGCTGCTGACGCTGTTGGTAGCGTTCATCACTCATTGGGAAATTCCTGGTTTACGGCCCGGTTGGGCGTCAAAGGTCATGCCGGTTTTACGGCGGCTGTCATCGCCCATCAGCCAGAGGTAGAGCGGCATGATATCAGCGGGGGTTTT‑AAGTTTCTGTGGATCTTCGGTCGGGAAGGCGCTGGCCCGCATTTCgggg                                                                                                                                                                                                     >  1:268175/1‑249 (MQ=255)
                                                                 tgggcctgggccACACGAGCATCTACTTTTTCTTTCACTCGCTGCTGACGCTGTTGGTAGCGTTCATCACTCATTGGGAAATTCCTGGTTTACGGCCCGGTTGGGCGTCAAAGGTCATGCCGGTTTTACGGCGGCTGTCCTCGCCCCTCCGCCAGAGGGAGAGGGGCATGATATTAGCGGGGGGTTT‑TAGTTTTTGTGGGTCTTCGGTCGGGGAGGGGGCGGCACGCCGTGCGGGGCGCGGACCGCCTggg                                                                                                                                                                                      >  2:131161/1‑251 (MQ=255)
                                                                  gggcctgggccACACGAGCATCTACTTTTTCTTTCACTCGCTGCTGACGCTGTTGGTAGCGTTCATCACTCATTGGGAAATTCCTGGTTTACGGCCCGGTTGGGCGTCAAAGGTCATGCCGGTTTTACGGCGGCTGTCATCGCCCATCAGCCAGAGGTAGAGCGGCATGATATCAGCGGGTGTTTT‑AAGTTTCTGTGGATCTTCGGTCGGGAAGGCGCTGGCACGCATTGCGGTGCGCGTACCGCCTgg                                                                                                                                                                                       <  1:41214/249‑1 (MQ=255)
                                                                        gggccACACGAGCATCTACTTGTTCTTTCACTCGCTGCTGACGCTGTTGGTAGCGTTCATCACTCATTGGGAAATTCCTGGTTTACGGCCCGGTTGGGCGTCAAAGGTCATGCCGGTTTTACGGCGGCTGTCATCGCCCATCAGCCAGAGGTAGAGCGGCATGATATCAGCGGGTGTTTT‑AAGTTTCTGTGGATCTTCGGTCGGGAAGGCGCTGGCACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGc                                                                                                                                                                               <  2:489850/251‑1 (MQ=255)
                                                                         ggccACACGAGCATCTACTTTTTCTTTCACTCGCTGCTGACGCTGTTGGTAGCGTTCATCACTCATTGGGAAATTCCTGGTTTACGGCCCGGTTGGGCGTCAAAGGTCATGCCGGTTTTACGGCGGCTGGCATCTCCCATCAGCCCGAGGTAGAGCGGCATGATATCAGCGGGTGTTTT‑TAATTTCTGTGGATCCTCGGGCGGGAAGGCGCTGGTACGCATTGCGGGgcgggggccccccgggtgtatgcc                                                                                                                                                                              >  2:94366/1‑230 (MQ=255)
                                                                               acGAGCATCTACTTTTTCTTTCACTCGCTGCTGACGCTGTTGGTAGCGTTCATCACTCATTGGGAAATTCCTGGTTTACGGCCCGGTTGGGCGTCAAAGGTCATGCCGGTTTTACGGCGGCTGTCATCGCCCATCAGCCAGAGGTAGAGCGGCATGATATCAGCGGGGGTTTT‑TAGTTTCTGTGGATCTTCGGTCGGGGAGGCGCTGGCACGCATTGCGGTGCGCGTACCGCCTGGGGTAATGCAGTTGac                                                                                                                                                                        >  1:186267/1‑251 (MQ=255)
                                                                                    cATCTACTTTTTCTTTCACTCGCTGCTGACGCTGTTGGTAGCGTTCATCACTCATTGGGAAATTCCTGGTTTACGGCCCGGTTGGGCGTCAAAGGTCATGCCGGTTTTACGGCGGCTGTCATCGCCCATCAGCCAGAGGTAGAGCGGCATGATATCAGCGGGTGTTTT‑AAGTTTCTGTGGATCTTCGGTCGGGAAGGCGCTGGCACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCa                                                                                                                                                                   <  1:174310/251‑1 (MQ=255)
                                                                                    cATCTACTTTTTCTTTCACTCGCTGCTGACGCTGTTGGTAGCGTTCATCACTCATTGGGAAATTCCTGGTTTACGGCCCGGTTGGGCGTCAAAGGTCATGCCGGTTTTACGGCGGCTGTCATCGCCCATCAGCCAGAGGTAGAGCGGCATGATATCAGCGGGGGTTTT‑TAGTTTCTGTGGATCTTCGGGCGGGGAGGCGCTGGCCCGCCTTTCGGTGCGCGTACCGCCTGGGGTAATGCAGGTGACACGc                                                                                                                                                                    >  1:232800/1‑250 (MQ=255)
                                                                                        tACTTTTTCTTTCACTCGCTGCTGACGCTGTTGGTAGCGTTCATCACTCATTGGGAAATTCCTGGTTTACGGCCCGGTTGGGCGTCAAAGGTCATGCCGGTTTTACGGCGGCTGTCATCGCCCATCCGCCCGAGGTAGAGCGGCATGATATCAGCGGGGGTTTT‑AAGTTTCTGTGGATCTTCGGTCGGGGAGGGGCTGGCACGCATTGCGGGGCGCGTACCGCCTGGGTTAATGCAGTTGAAACGCAGgcg                                                                                                                                                               >  1:370730/1‑251 (MQ=255)
                                                                                        tACTTTTTCTTTCACTCGCTGCTGACGCTGTTGGTAGCGTTCATCACTAATTGGGAACTTCATGGTTTACGGCCCGGTTGGGCGTCAAAGGTCATGCCGGTTTTACGGCGGCTGTCATCGCCCCTCAGCCCGAGGTAGAGCGGCATGATATCAGCCGGGGTTTT‑AAGTTTCTGTGGGACTTCGGTAGGGGAGGGGCTGGGCCGCCTTGCCGTGCGGGTACCGCCCGGGGTTATGCCGTTGACACGCAGgcg                                                                                                                                                               >  1:565088/1‑251 (MQ=255)
                                                                                                ctttcACTCGCTGCTGACGCTGTTGGTAGCGTTCATCACTCATTGGGAAATTCCTGGTTTACGGCCCGGTTGGGCGTCAAAGGTCATGCCGGTTTTACGGCGGCTGTCATCGCCCATCAGCCCGAGGTAGAGCGGCATGGTATCAGCGGGTGTTTT‑TAGTTTTTGTGGGTATTCGGTCGGGGAAGCGCTGGCACGCACTTCGGGGCGCGGACCGCCCGGGTTAATGCCGTTGACAAGCCGGCGCCGCTGat                                                                                                                                                       >  2:158733/1‑251 (MQ=255)
                                                                                                  ttcACTCGCTGCTGACGCTGTTGGTAGCGTTCATCACTCATTGGGAAATTCCTGGTTTACGGCCCGGTTGGGCGTCAAAGGTCATGCCGGTTTTACGGCGGCTGTCATCGCCCATCAGCCAGAGGGAGAGCGGCATGATATCAGCGGGTGTTTT‑TAGTTTTTGTGGATCTTCGGTCGGGAAGGCGCTGGCACGCATTTCGGGGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCCGatat                                                                                                                                                     >  1:344708/1‑251 (MQ=255)
                                                                                                     aCTCGCTGCTGACGCTGTTGGTAGCGTTCATCACTCATTGGGAAATTCCTGGTTTACGGCCCGGTTGGGCGTCAAAGGTCATGCCGGTTTTACGGCGGCTGTCATCGCCCATCAGCCAGAGGTAGAGCGGCATGATATCAGCGGGTGTTTT‑AAGTTTCTGTGGATCTTCGGTCGGGAAGGCGCTGGCACTCATTGCGGTGCGCGTACCGCCTGGGTTAATTCAGTGGACACGCAGGTGCTGCTGAGATTCa                                                                                                                                                  <  2:259381/251‑1 (MQ=255)
                                                                                                            gctgACGCTGTTGGTAGCGTTCATCACTCATTGGGAAATTCCTGGTTTACGGCCCGGTTGGGCGTCAAAGGTCATGCCGGTTTTACGGCGGCTGTCATCGCCCATCAGCCAGAGGTAGAGCGGCATGATATCAGCGGGTGTTTT‑AAGTTTCTGTGGATCTTCGGTCGGGAAGGCGCTGGCACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCAT‑CGGc                                                                                                                                            <  2:268175/249‑1 (MQ=255)
                                                                                                                  gCTGTTGGTAGCGTTCATCACTCATTGGGAAATTCCTGGTTTACGGCCCGGTTGGGCGTCAAAGGTCATGCCGGTTTTACGGCGGCTGTCATCGCCCATCAGCCAGAGGTAGAGCGGCATGATATCAGCGGGTGTTTT‑AAGTTTCTGTGGATCTTCGGTCGGGAAGGCGCTGGCACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCAT‑CGGCCAGTACCt                                                                                                                                    <  1:399895/251‑1 (MQ=255)
                                                                                                                         gTAGCGTTCATCACTCATTGGGAAATTCCTGGTTTACGGCCCGGTTGGGCGTCAAAGGTCATGCCGGTTTTACGGCGGCTGTCATCGCCCATCAGCCAGAGGTAGAGCGGCATGATATCAGCGGGGGTTTT‑AAGTTTCTGTGGATCTTCGGTCGGGAAGGCGCTGGCACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCAT‑CGGCCCGGACCCGcatcct                                                                                                                             >  1:126182/1‑249 (MQ=255)
                                                                                                                           aGCGTTCATCACTCATTGGGAAATTCCTGGTTTACGGCCCGGTTGGGCGTCAAAGGTCATGCCGGTTTTACGGCGGCTGTCATCGCCCATCAGCCAGAGGTAGAGCGGCATGATATCAGCGGGTGTTTT‑AAGTTTCTGTGGATCTTCGGTCGGGAAGGCGCTGGCACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCAT‑CGGCCAGTACCTGCATCATcc                                                                                                                           <  1:429529/251‑1 (MQ=255)
                                                                                                                                tcatcaCTCATTGGGCAATTCCTGGTTTACGGCCCGGTTGGGCGTCAAAGGTCATGCCGGTTTTACGGCGGCTGTCATCGCCCAACAGCCAGAGGTAGAGCGGCATGATATCAGCGGGTGTTTT‑AAGTTTCTGTGGATCTTCGGTCGGGGAGGCGCTGGCACGCCTTGCGGGGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCGGCTGATATTCAT‑CGGCCAGTACCTGCATCATCCCTTCg                                                                                                                      >  1:511142/1‑251 (MQ=255)
                                                                                                                                   tcaCTCATTGGGAAATTCCTGGTTTACGGCCCGGTTGGGCGTCAAAGGTCATGCCGGTTTTACGGCGGCTGTCATCGCCCATCAGCCAGAGGTAGAGCGGCATGATATCAGCGGGGGGTTT‑TAGTTTTTGTGGGTCTTCGGGCGGGGAGGGGCGGGGACGCATTGCGGGGGGCGGACCGCCTGGGGTAATGCAGGTTACACGCAGGGGCTGCTGATATTCATTCGGCCAGAACC‑GCATCATCCCTTCGGTg                                                                                                                   >  2:90309/1‑251 (MQ=255)
                                                                                                                                       tCATTGGGAAATTCCTGGTTTACGGCCCGGTTGGGCGTCAAAGGTCATGCCGGTTTTACGGCGGCTGTCTTCGCCCATCAGCCCGAGGTAGAGCGGCATGATAACAGCGGGTGGTTT‑AAGTTTCTGTGGATCTTCGGGCGGGAAAGGGCTGGCTCGCAATTCGGTGCGCGGACCGCCTGGGGTTATGGAATTTACACGCAGGGGCTGCTGATATTCAT‑CGGCCAGTAACTGTATCCTCCCTTCggggggc                                                                                                                >  1:318435/1‑245 (MQ=255)
                                                                                                                                        cATTGGGAAATTCCTGGTTTACGGCCCGGTTGGGCGTCAAAGGTCATGCCGGTTTTACGGCGGCTGTCATCGCCCATCAGCCAGAGGTAGAGCGGCATGATATCAGCGGGTGTTTT‑AAGTTTCTGTGGATCTTCGGTCGGGAAGGCGCTGGCACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCAT‑CGGCCAGTACCTGCATCATCCCTTCGGTGG‑Ca                                                                                                               <  1:255737/249‑1 (MQ=255)
                                                                                                                                              gAAATTCCTGGTTTACGGCCCGGTTGGGCGTCAAAGGTCATGCCGGTTTTACGGCGGCTGTCATCGCCCATCAGCCAGAGGTAGAGCGGCATGATATCAGCGGGTGTTTT‑AAGTTTCTGTGGATCTTCGGTCGGGAAGGCGCTGGCACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCAT‑CGGCCAGTACCTGCATCATCCCTTCGGTGG‑CAAATTTCGa                                                                                                       <  2:126182/251‑1 (MQ=255)
                                                                                                                                                             cGGCCCGGTTGGGCGTCAAAGGTCATGCCGGTTTTACGGCGGCTGTCATCGCCCATCAGCCAGAGGTAGAGCGGCCTGATATCAGCGGGGGTTTTTAAGTT‑CTGTGGATCTTCGGTCGGGGAGGCGCTGGCCCGCATTGCGGTGCGCGTAACGCCTGGGGTAATGCAGTTGAAACGCAGGCGCTGCTGATATTCAT‑CGGCCCGTACCCGCATCATCCCCTCGGGGG‑CAAATTTCGACGCTGC‑ATATGCAcc                                                                                       >  1:28389/1‑251 (MQ=255)
                                                                                                                                                             cGGCCCGGTTGGGCGTCAAAGGTCATGCCGGTTTTACGGCGGCTGTCATCGCCCATCAGCCAGAGGTAGAGCGGCATGATATCAGCGGGTGTTTT‑AAGTTTCTGTGGATCTTCGGTCGGGAAGGCGCTGGCACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCAT‑CGGCCAGTACCTGCATCATCCCTTCGGTGG‑CAAATTTCGACGCTGC‑ATATGCAcc                                                                                       <  2:199354/251‑1 (MQ=255)
                                                                                                                                                              ggCCCGGTTGGGCGTCAAAGGTCATGCCGGTTTTACGGCGGCTGTCATCGCCCATCAGCCAGAGGTAGAGCGGCATGATATCAGCGGGTGTTTT‑AAGTTTCTGTGGATCGTCGGTCGGGAAGGCGCTGGCACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCATGCGCTGCTGATATTCAT‑CGGCCAGTACCTGCATCATCCCTGCGGTGG‑CAAATTTCGACGCTGC‑ATATGCAccc                                                                                      <  2:328852/251‑1 (MQ=255)
                                                                                                                                                                 ccGGTTGGGCGTCAAAGGTCATGCCGGTTTTACGGCGGCTGTCATCGCCCATCAGCCAGAGGTAGAGCGGCATGATATCAGCGGGGGGTTT‑AAGTTTCTGTGGATCTTCGGGCGGGAAGGGGCTGGCACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTTACACGCAGGCGCTGCTGATATTCAT‑CGGCCAGTAACTGCATCATCCCTTCGGGGG‑CAAAATTTGACGCTGG‑ATATGCACCCCa                                                                                    >  1:509187/1‑250 (MQ=255)
                                                                                                                                                                 ccGGTTGGGCGTCAAAGGTCATGCCGGTTTTACGGCGGCTGTCACTGCCCATCAGCCAGAGGGGGAGCGGCATGCTATCAGGGGGGGTTTT‑AAAGTTCTGTGGATCCTCGGTCGGGGAGGCGGTGGCACGCAATGCGGTGGGCGTATCTCCTGGGGTAATGCAGTTGACCCGCAGACGCTGGTGATATTAAT‑CGGCCAGGACCTGCCTCATCCCCTCGGGGG‑GCAACTTCGACCCCGC‑CTATGCACCCCac                                                                                   >  2:463445/1‑250 (MQ=255)
                                                                                                                                                                  cGGTTGGGCGTCAAAGGTCATGCCGGTTTTACGGCGGCTGTCATCGCCCATCAGCCAGAGGTAGAGCGGCATGATATCAGCGGGTGTTTT‑AAGTTTCTGTGGATCTTCGGTCGGGAAGGCGCTGGCACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCAT‑CGGCCAGTACCTGCATCATCCCTTCGGTGG‑CAAATTTCGACGCTGC‑ATATGCACCCCAGt                                                                                  <  2:94412/251‑1 (MQ=255)
                                                                                                                                                                    gTTGGGCGTCAAAGGTCATGCCGGTTTTACGGCGGCTGTCATCGCCCATCAGCCAGAGGTAGAGCGGCATGATATCAGCGGGGGGTTT‑AAATTTCTGTGGGTCTTCGGTCGGGAAGGCGCTGGCACGCATTGCGGTTCGCGTACCGCCTGGGTTAATTCAGTTTACACCCAGGCGCTGCTGCTATTCCA‑CGGCCAGTACCTGCCTCCTCCCTTCGGGGGGCAATTTTC‑ACGCCTC‑ATATGCACCCCCGttt                                                                                >  2:515294/1‑250 (MQ=255)
                                                                                                                                                                      tGGGCGTCAAAGGTCATGCCGGTTTTACGGCGGCTGTCATCGCCCATCAGCCAGAGGTAGAGCGGCATGATATCAGCGGGTGTTTT‑AAGTTTCTGTGGATCTTCGGTCGGGAAGGCGCTGGCACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCAT‑CGGCCAGTACCTGCATCATCCCTTCGGTGG‑CAAATTTCGACGCTGC‑ATATGCACCCCAGTTgg                                                                               <  1:111009/250‑1 (MQ=255)
                                                                                                                                                                        ggCGTCACAGGTCATGCCGGTTTTACGGCGGCTGTCATCGCCCATCAGCCAGAGGGAGAGCGGCATGATATCAGCGGGGGGTTT‑AAGGTTCTGTGGATCTTCGGGCGGGAAGGCGCTGGCACGCATTGCGGTGCGCGTACCGCCTGGGGTAATGCAGTTGACACGCAGGCGCTGCTGATATTCCT‑CGGCCAGTACCTGCATCATCCCTTCGGGGG‑CAAATTTTGAcagatcgggagagcacacgtctgga                                                                              >  1:149627/1‑226 (MQ=255)
                                                                                                                                                                        ggCGTCAAAGGTCATGCCGGTTTTACGGCGGCTGTCATCGCCCATCAGCCAGAGGTAGAGCGGCATGATATCAGCGGGTGTTTT‑AAGTTTCTGTGGATCTTCGGTCGGGAAGGCGCTGGCACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCAT‑CGGCCAGTACCTGCATCATCCCTTCGGTGG‑CAAATTTCGAc                                                                                                      <  2:149627/226‑1 (MQ=255)
                                                                                                                                                                             cAAAGGTCATGCCGGTTTTACGGCGGCTGTCATCGCCCATCAGCCAGAGGTAGAGCGGCATGATATCAGCGGGTGTTTC‑AAGTTTCTGTGGATCTTCGGTCGGGAAGGCGCTGGCACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCAT‑CGGCCAGTACCTGCATCATCCCTTCGGTGG‑CAAATTTCGACGCTGC‑ATATGCACCCCAGTTGGCTCGTcc                                                                        <  2:170450/250‑1 (MQ=255)
                                                                                                                                                                              aaaGGTCATGCCGGTTTTACGGCGGCTGTCATCGCCCATCAGCCAGAGGTAGAGCGGCATGATATCAGCGGGTGTTTT‑AAGTTTCTGTGGATCTTCGGTCGGGAAGGCGCTGGCACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCAT‑CGGCCAGTACCTGCATCATCCCTTCGGTGG‑CAAATTTCGACGCTGC‑ATATGCACCCCAGTTGGCTCGTCCCt                                                                      <  1:90309/251‑1 (MQ=255)
                                                                                                                                                                                   tCATGCCGGTTTTACGGCGGCTGTCATCGCCCATCAACCCGAGGTAGAGCGGCATGATATCAGCGGGGGGTTT‑TAGTTTTTTTGGGTCTTCGGTCGGGGAGGCGCTGGCACGCATTGCGGCGCGCGTACCTCCGGGGGTAATTCAGGTGACACGCAGGCGCTGCTGGTATTCAT‑CCGCCCGTACCCGCATCATCCCCTCGGGGG‑GACATTTTGAAGCCGCCAT‑TTCACCCCCGTTGGCTCGTcccccgcgg                                                                 >  2:413964/1‑245 (MQ=255)
                                                                                                                                                                                        ccGGTTTTACGGCGGCTGTCATCGCCCATCAGCCAGAGGTAGAGCGGCATGATATCAGCGGGTGTTTT‑AAGTTTCTGTGGATCTTCGGTCGGGAAGGCGCTGGCACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCAT‑CGGCCAGTACCTGCATCATCCCTTCGGTGG‑CAAATTTCGACGCTGC‑ATATGCACCCCAGTTGGCTCGTCCCTGACGTCCAAc                                                            <  2:56921/251‑1 (MQ=255)
                                                                                                                                                                                           gTTTTACGGCGGCTGTCATCGCCCATCAGCCAGAGGTAGAGCGGCATGATATCAGCGGGTGTTTT‑AAGTTTCTGTGGATCTTCGGTCGGGAAGGCGCTGGCACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCAT‑CGGCCAGTACCTGCATCATCCCTTCGGTGG‑CAAATTTCGACGCTGC‑ATATGCACCCCAGTTGGCTCGTCCCTGACGTCCAACGCt                                                         <  2:65357/251‑1 (MQ=255)
                                                                                                                                                                                                  ggcggcTGTCATCGCCCATCAGCCAGAGGTAGAGCGGCATGATATCAGCGGGTGTTTT‑AAGTTTCTGTGGATCTTCGGTCGGGAAGGCGCTGGCACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCAT‑CGGCCAGTACCTGCATCATCCCTTCGGTGG‑CAAATTTCGACGCTGC‑ATATGCACCCCAGTTGGCTCGTCCCTGACGTCCAACGCTtgatga                                                   >  1:65878/1‑250 (MQ=255)
                                                                                                                                                                                                      gcTGTCATCGCCCATCAGCCAGAGGTAGAGCGGCATGATATCAGCGGGTGTTTT‑AAGTTTCTGTGGATCTTCGGTCGGGAAGGCGCTGGCACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCAT‑CGGCCAGTACCTGCATCATCCCTTCGGTGG‑CAAATTTCGACGCTGC‑ATATGCACCCCAGTTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTaa                                              <  1:224061/251‑1 (MQ=255)
                                                                                                                                                                                                               gCCCATCAGCCAGAGGTAGAGCGGCATGATATCAGCGGGTGTTTT‑AAGTTTCTGTGGATCTTCGGTCGGGAAGGCGCTGGCACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCAT‑CGGCCAGTACCTGCATCATCCCTTCGGTGG‑CAAATTTCGACGCTGC‑ATATGCACCCCAGTTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTg                                     >  1:177606/1‑251 (MQ=255)
                                                                                                                                                                                                                  cATCAGACAGAGGTAGAGCGGCATGATATCAGCGGGGGTTTT‑AAGTTTCTGTGGATCTTCGGTCGGGAAGGCGCTGGCACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCAT‑CGGCCAGTACCTGCATCATCCCTTCGGTGG‑CAAATTTCGACGCTGC‑ATATGCACCCCAGTTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAAc                                  >  2:517688/1‑251 (MQ=255)
                                                                                                                                                                                                                     cAGCCAGAGGTAGAGCGGCATGATATCAGCGGGTGTTTT‑AAGTTTCTGTGGATCTTCGGTCGGGAAGGCGCTGGCACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCAT‑CGGCCAGTACCTGCATCATCCCTTCGGTGG‑CAAATTTCGACGCTGC‑ATATGCACCCCAGTTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCgg                               <  2:264493/251‑1 (MQ=255)
                                                                                                                                                                                                                         cAGAGGTAGAGCGGCATGATATCAGCGGGTGTTTT‑AAGTTTCTGTGGATCTTCGGTCGGGAAGGCGCTGGCACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCAT‑CGGCCAGTACCTGCATCATCCCTTCGGTGG‑CAAATTTCGACGCTGC‑ATATGCACCCCAGTTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTc                           >  1:348976/1‑251 (MQ=255)
                                                                                                                                                                                                                            agGTAGAGCGGCATGATATCAGCGGGTGTTTT‑AAGTTTCTGTGGATCTTCGGTCGGGAAGGCGCTGGCACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCAT‑CGGCCAGTACCTGCATCATCCCTTCGGTGG‑CAAATTTCGACGCTGC‑ATATGCACCCCAGTTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCg                         >  1:455547/1‑250 (MQ=255)
                                                                                                                                                                                                                              gTAGAGCGGCATGATATCAGCGGGTGTTTT‑AAGTTTCTGTGGATCTTCGGTCGGGAAGGCGCTGGCACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCAT‑CGGCCAGTACCTGCATCATCCCTTCGGTGG‑CAAATTTCGACGCTGC‑ATATGCACCCCAGTTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGAtt                      <  2:264965/251‑1 (MQ=255)
                                                                                                                                                                                                                                      gCATGATATCAGCGGGTGTTTT‑AAGTTTCTGTGGATCTTCGGTCGGGAAGGCGCTGGCACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCAT‑CGGCCAGTACCTGCATCATCCCTTCGGTGG‑CAAATTTCGACGCTGC‑ATATGCACCCCAGTTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGtaat              <  1:257975/251‑1 (MQ=255)
                                                                                                                                                                                                                                      gCATGATATCAGCGGGTGTTTT‑AAGTTTCTGTGGATCTTAGGTCGGGAAGGCGCTGGCACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCAT‑CGGCCAGTACCTGCATCATCCCTTCGGTGG‑CAAATTTCGACGCTGC‑ATATGCACCCCAGTTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGGCCGATTTGAGtaat              >  1:431905/1‑251 (MQ=255)
                                                                                                                                                                                                                                        aTGATATCAGCGGGTGTTTT‑AAGTTTCTGTGGATCTTCGGTCGGGAAGGCGCTGGCACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCAT‑CGGCCAGTACCTGCATCATCCCTTCGGTGG‑CAAATTTCGACGCTGC‑ATATGCACCCCAGTTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGtaata             >  1:392522/1‑250 (MQ=255)
                                                                                                                                                                                                                                               cAGCGGGTGTTTT‑AAGTTTCTGTGGATCTTCGGTCGGGAAGGCGCTGGCACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCAT‑CGGCCAGTACCTGCATCATCCCTTCGGTGG‑CAAATTTCGACGCTGC‑ATATGCACCCCAGTTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGc     <  1:515294/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  cGGGTGTTTT‑AAGTTTCTGTGGATCTTCGGTCGGGAAGGCGCTGGCACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCAT‑CGGCCAGTACCTGCATCATCCCTTCGGTGG‑CAAATTTCGACGCTGC‑ATATGCACCCCAGTTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGt  <  1:327811/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                          
CAAATGCCGCGGTGGTTTTGCCTTTTCCATTGCCGGTAAAGACGATGATAATACCGCGTTCATCCTGGGCCTGGGCCACACGAGCATCTACTTTTTCTTTCACTCGCTGCTGACGCTGTTGGTAGCGTTCATCACTCATTGGGAAATTCCTGGTTTACGGCCCGGTTGGGCGTCAAAGGTCATGCCGGTTTTACGGCGGCTGTCATCGCCCATCAGCCAGAGGTAGAGCGGCATGATATCAGCGGGTGTTTT‑AAGTTTCTGTGGATCTTCGGTCGGGAAGGCGCTGGCACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCAT‑CGGCCAGTACCTGCATCATCCCTTCGGTGG‑CAAATTTCGACGCTGC‑ATATGCACCCCAGTTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGT  >  NC_000913/1328219‑1328711

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