Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913505,4400CA31.2% 38.9 / 10.6 32A176E (GCA→GAA) ushA5'‑nucleotidase/UDP‑sugar hydrolase
Reads supporting (aligned to +/- strand):  ref base C (13/9);  new base A (0/10);  total (13/19)
Fisher's exact test for biased strand distribution p-value = 1.65e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACTTACAGGATGCCGAACCTGATTTTCGCGGTATGAATCTGGTGGGCTATGACGCGATGGCGATCGGTAATCATGAATTTGATAATCCGCTCACCGTATTACGCCAGCAGGAAAAGTGGGCCAAGTTCCCGTTGCTTTCCGCGAATATCTACCAGAAAAGTACTGGCGAGCGCCTGTTTAAACCGTGGGCGCTGTTTAAGCGTCAGGATCTGAAAATTG‑CCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAA‑TACTTCACTGATATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTACGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCGATGATCGTCGGTGGTCACTCGCAAGATCCGGTCT  >  NC_000913/505191‑505685
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                      
aCTTACAGGATGCCGAACCTGATTTTCGCGGTATGAATCTGGTGGGCTATGACGCGATGGCGATCGGTAATCATGAATTTGATAATCCGCTCACCGTATTACGCCAGCAGGAAAAGTGGGCCAAGTTCCCGTTGCTTTCCGCGAATATCTACCAGAAAAGTACTGGCGAGCGCATGTTTAAACCGTGGGCGCTGTTTAAGCGTCAGGATCTGAAAAGTG‑CCGTTATTGGGCTGACAACCGATGACACAGCa                                                                                                                                                                                                                                                       <  2:270819/251‑1 (MQ=255)
   tACAGGATGCCGAACCTGATTTTCGCGGTATGAATCTGGTGGGCTATGACGCGATGGCGATCGGTAATCATGAATTTGATAATCCGCTCACCGTATTACGCCAGCAGGAAAAGTGGGCCAAGTTCCCGTTGCTTTCCGCGAATATCTACCAGAAAAGTACTGGCGAGCGCCTGTTTAAACCGTGGGCGCTGTTTAAGCGTCAGGATCTGAAAATTG‑CCGTTATTGGGCTGACAACCGATGACACAGCaa                                                                                                                                                                                                                                                      <  2:93629/249‑1 (MQ=255)
             cGAACCTTAGTTTCGCGGTATGAATCTGGTGTGCGATGACGCGATGTCGATCTGTAAACATGAATTTTCTAATCCGCTCACCGTATTACGCCAGCAGGAAAAGTTGGCCAAGTTCCCGTTTCTTTCCGCGAATATCTACCACAAAATGACTCGCGAGCGCCTGTTTAAACCGTGGGCGCTGTTTAAGCTTCAGGATCTGAAAATGG‑CCGTGAGTGGGCGGACAACCGAGGACACAGCAAAAATTGGTAAc                                                                                                                                                                                                                                           <  2:369041/250‑1 (MQ=255)
              gAACCTGATTTTCGCGGTATGAATCTGGTGGGCTATGACGCGATGGCGATCGGTAATCATGAATTTGATAATCCGCTCACCGTATTACGCCAGCAGGAAAAGTGGGCCAAGTTCCCGTTGCTTTCCGCGAATATCTACCAGAAAAGTACTGGCGAGCGCCTGTTTAAACCGTGGGCGCTGTTTAAGCGTCAGGATCTGAAAATTG‑CCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAAccc                                                                                                                                                                                                                                         >  2:365392/1‑251 (MQ=255)
                aCCTGATTTTCGCGGTATGAATCTGGTGGGCTATGACGCGATGGCGATCGGTAATCATGAATTTGATAATCCGCTCACCGTATTACGCCAGCAGGAAAAGTGGGCCAAGTTCCCGTTGCTTTCCGCGAATATCTACCAGAAAAGTACTGGCGAGCGCCTGTTTAAACCGTGGGCGCTGTTTAAGCGTCAGGATCTGAAAATTG‑CCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCgg                                                                                                                                                                                                                                       >  1:247986/1‑251 (MQ=255)
                         tCGCGGTATGATTCTGGTGGGCTATGACGCGATGGCGATCGTTAATCATGAATTTGATAATCCGCTCACCGTATTACGCCAGCAGGAAAAGTGGGCCAAGTTCCCGTTGCTTTCCGCGAATATCTACCAGAAAAGTACTGGCGAGCGCCTGTTTAAACCGTGGGCGCTGTTTAAGCGTCAGGATCTGAAAATTG‑CCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAA‑TACTTc                                                                                                                                                                                                                              <  1:221443/250‑1 (MQ=255)
                          cgcgGTATGAATCTGGTGGGCTATGACGCGATGGCGATCGGTAATCATGAATTTGATAATCCGCTCACCGTATTACGCCAGCAGGAAAAGTGGGCCAAGTGCCCGTTGCTTTCCGCGAATATCTACCAGAAAAGTACTGGCGAGCGCCTGTTTAAACCGTGGGCGCTTTTTAAGCGTCAGGATCTGAAAATTG‑CCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAA‑TACTTCa                                                                                                                                                                                                                             <  1:371123/250‑1 (MQ=255)
                          cgccgTATGAATCTGGTGGGCTATGACGCGATGGCGATCGGTAATCATGAATTTGATAAGCCGCTCACCGTATTACGCCAGCAGGAAAAGTGGGCCAAGTTCCCGTTGCTTTCCGCTAATATCTACCAGAAAAGTACTGGCGAGCGCCTGTTAAAACCGTGGGCTCTGTTTAAGCGTCAGGATCTTAAAATTG‑CCGTTATTTTGCGGACAACCGATGACACAGAAAAAATTGGTAACCCGGAA‑TGCTTCa                                                                                                                                                                                                                             <  2:208718/246‑1 (MQ=255)
                                           gtgCTATGGCGCGATCGCGATCGGTAATCATGAATTTGATATTCCGCTCACCGTATTACGCCAGCATGAAAAGTTGGCTAAGTGCCCGTTGCTTTCCGCGAATATCTACCAGAAAAGGACTGGCGATCGACTGTTTACAACGTGGTCTCTGTGTAAGCGTCATGATCGGAAAATGG‑CCGTTATTGAGCTGACAACCGATGACACATCAAAAATTGGTAATCCGTAA‑TTCTTCACTGATATCGAATTTCGTa                                                                                                                                                                                                           <  2:139344/249‑1 (MQ=255)
                                                         tGGCGATCGGTAATCATGAATTTGATAATCCGCTCACCGTATTACGCCAGCAGGAAAAGTGGGCCAAGTTCCCGTTGCTTTCCGCGAATATCTACCAGAAAAGTACTGGCGAGCGCCTGTTTAAACCGTGGGCGCTGTTTAAGCGTCAGGATCTGAAAATTG‑CCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAA‑TACTTCACTGATATCGAATTTCGTAAGCCCGCCGATga                                                                                                                                                                                              >  1:129830/1‑250 (MQ=255)
                                                            cGATCGGTAATCATGAATTTGATAATCCGCTCACCGTATTACGCCAGCAGGAAAAGTGGGCCAAGTTCCCGTTGCTTTCCGCGAATATCTACCAGAAAAGTACTGGCGAGCGCCTGTTTAAACCGTGGGCGCTGTTTAAGCGTCAGGATCTGAAAATTG‑CCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAA‑TACTTCACTGATATCGAATTTCGTAAGCCCGCCGATgaagcg                                                                                                                                                                                          >  1:46274/1‑251 (MQ=255)
                                                                 ggTAATCATGAATTTGATAATCCGCTCACCGTATTACGCCAGCAGGAAAAGTGGGCCAAGTTCCCGTTGCTTTCCGCGAATATCTACCAGAAAAGTACTGGCGAGCGCCTGTTTAAACCGTGGGCGCTGTTTAAGCGTCAGGATCTGAAAATTG‑CCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAA‑TACTTCACTGATATCGAATTTCGTAAGCCCGCAGATGAAGCGAAGCt                                                                                                                                                                                     >  1:342026/1‑251 (MQ=255)
                                                                       aaTGAATGTGATAACCCTCTCACCGTATTACGCCAGCAGGAAAAGGGGGCCAATTCCCCGTTGCTTTCCGCGAATATCTACCAGAAAAGTGCTGGGGAGCGCCTGTTAAACCCGTGGGCGTTGTTAAAGCGTCAGGATCTGAAAATTG‑CCGTTATGGGGCTGACAACCGATGACACATAAAAAATTGGTAACCCGAAA‑TACTTCACTGATATCGATGTTCTTAAGCCCGCCGATGAAGCGAAGCTGGTGa                                                                                                                                                                                <  1:321240/249‑1 (MQ=255)
                                                                                   aaTCCGCTCACCGTATTACGCCAGCAGGAAAAGTGGGCCAAGTTCCCGTTGCTTTCCGCGAATATCTACCAGAAAAGTACTGGCGAGCGCCTGTTTAAACCGTGGGCGCTGTTTAAGCGTCAGGATCTGAAAATTG‑CAGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAA‑TACTTCACTGATATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGc                                                                                                                                                                    >  1:412163/1‑250 (MQ=255)
                                                                                                          gcagAAAAAGGGGGCAAAGTCCCCGTTTCTTCCCGAGAATATCTACAAGAAAGGTTCGGGGGAGCGCCTTGTTAAACGGGGGGCGTTGTTTAAGCGTGAGGATCTGAAAATGG‑CCTTAATTGGGGGAACACCCGGGGAAACAGAAAAAATTGGAGCCCCGGAA‑GACTTCACGGATATGGAATTTAGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACatta                                                                                                                                            <  1:247354/251‑1 (MQ=255)
                                                                                                            agaaaaggtgggccaatcccccGTTGCTTCCCGCAAATATCTACCAAAAAAGTACTGGCGAGCCCTTTTTTAAACCGGGGGCGTTGTTTAAGCGTCAGGTTCT‑AAAATTGGCCTTTATTGGGTTGCCAACCTAGCAACCAAAAAAAATTGTTCACCCGAAA‑TACTCCTTTGTTTTCAAATTTCTAAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACattatt                                                                                                                                          <  1:36020/232‑1 (MQ=255)
                                                                                                                   gggggCCAAGTTCCCGTTGCTTCCCGAGAATATCTACAAGAAAAGTACTGGCGGGCGCCTTTTTAAACCGGGGGCGTTGTTAAGGCGTCAGGATCTGAAAATTG‑CCGTTATTGGGCTGACACCCGATGACACAGCAAAAATTGGTACCCCGGAA‑TACTTCACTGATATCGAATTGCGAAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCgcgg                                                                                                                                   <  1:319783/249‑1 (MQ=255)
                                                                                                                                 ggttgttTTCCGCGAATTTCTCCCAGAAAATTACTGGCGAGCGCCTTTTTAACCGGTGGGCGCTTTTTAAGCGTCAGGTTCTAAAAATTG‑CCGTTATGGGGCTGACACCCGATGACACAGCAAAAATTGGTACCCCGGAA‑TCCTTCACTGATATCGATTTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATgggg                                                                                                                     <  1:44475/250‑1 (MQ=255)
                                                                                                                                  gTTGCTTCCCGAGAATATCTACCAGAAAGGTACTGGGGAGCGCTTGTTAAAACCGTGGGCGTTGTTTAAGCGTCAGTATCTAAAAATTG‑CCGTTATTGGTCTGACACCCAATGACACAGCAAAAATTGGTAACCCGGAA‑TACTTCACTGATACCAAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATggg                                                                                                                      <  2:2604/249‑1 (MQ=255)
                                                                                                                                      cTTTCCGCGAATATCTACCAGAAAAGTACTGGCGAGCGCCTGTTTAAACCGTGGGCGCTGTTTAAGCGTCAGGATCTGAAAATTG‑CCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAA‑TACTTCACTGATATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCAtt                                                                                                                 >  1:144042/1‑250 (MQ=255)
                                                                                                                                                atatcaaaaaaaaaaGTACTGGCGAGCGCTTGTTTAAACGGTGGGGGTTGTTTAAGCTACAGAATCGAAAAATTG‑CCGTAATGGGGCGGACAACCAATGAAACAGAAAAAATTGTAAACCCGAAA‑TACTTCACTAATAGAAAATTTCGAAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGTAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTACGATAATGGt                                                                                                      <  1:384643/240‑1 (MQ=255)
                                                                                                                                                                   gggCGACCGCCTTTTTAACCCGTGGGCGTTGTTTAAGCGTAAGGATCTAAAAATTG‑CCTTTATTGGGCTGACAACCGATGACACAGAAAAAATTGTTACCCCGAAA‑TACTTCACTGATATCGAATTCCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTACGATAATGGTGAGCACGGCTCTAACGCAc                                                                                   <  1:365392/250‑1 (MQ=255)
                                                                                                                                                                          cgcCTGTGTAAACCGTGGGCGCTGTTTAAGAGTCAGGATCTGACACTTG‑CCGTTATTGGGCTGACAACCGCTGACACAGCACAAATTGGTAACACGGAA‑TACTTCACTGATATCGAATTTCGTAAGCACGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTACGATAATGGTGCGCCAGGCTCTAACGCACCGGGCg                                                                             >  2:58190/1‑250 (MQ=255)
                                                                                                                                                                           gcCTTTTTAAACCGGGGGCGTTGTTTAAGCGGCAGAATATAAAAATGG‑CCGTTATTGGGCTGACACCCGATGACACAGAAAAAATTGGAACCCCGGAA‑TACTCCACGGATATCAAATTTCGAAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATgggg                                                                                                                     <  1:56627/209‑1 (MQ=255)
                                                                                                                                                                           gcCTGTTTAAACCGTGGGCGCTGTTTAAGCGTCAGGATCTGAAAATTG‑CCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAA‑TACTTCACTGATATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATgggg                                                                                                                     >  2:56627/1‑209 (MQ=255)
                                                                                                                                                                             cTGTTTAAACCGTGGGCGCTGTTTAAGCGTCAGGATCTGAAAATTG‑CCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAA‑TACTTCACTGATATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTACGAt                                                                                                            >  1:235996/1‑216 (MQ=255)
                                                                                                                                                                                            gggctgttaaagcgcgagaatcaaaaaaTGG‑CCGTTATGGGGCGGAAACCCGATGACACAGAAAAAATGGATAACCCGAAA‑TACTCCACTGATATCTAATTGCGAAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTACGATAATTGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCAcgcgc                                                          <  1:154567/227‑1 (MQ=255)
                                                                                                                                                                                                      aGCGTCAGGATCTGAAAATTG‑CCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAA‑TACTTCACTGATATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCAGGTGATTCAGGAGATGCAACAGCCAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTACGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGctgcctacc                                                >  1:418397/1‑248 (MQ=255)
                                                                                                                                                                                                          acagaaTCTGAAAATGG‑CCGTTATTGGGATGACACCCGATGACACAGAAAAAATTGGTA‑CCCGGAAATACTTCACTGATATCGATTTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTACGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGAt                                            <  1:246401/246‑1 (MQ=255)
                                                                                                                                                                                                                atttgaaaattgctttttgtGGGCTGACACCCGATGACACAGAAAAAATTGGAAACCCGGAA‑TACTTCACTGATATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTACGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTg                                       <  1:194672/232‑1 (MQ=255)
                                                                                                                                                                                                                                            aaCCGATGACACAGCAAAAATTGGTAACCCGGAA‑TACTTCACTGATATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTACGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCGATGATCGTCGGTGGTCACTCGCa             >  2:230142/1‑249 (MQ=255)
                                                                                                                                                                                                                                                     cacaGCAAAAATTGGTAACCCGGAA‑TACTTCACTGATATCGAATTTCGTAAGCCAGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCAGCGACCCATATGGGGCATTACGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGAGATGCTCGTCGGTGGTCACTCGCAAGATCCGGTCt  >  2:339894/1‑251 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                      
ACTTACAGGATGCCGAACCTGATTTTCGCGGTATGAATCTGGTGGGCTATGACGCGATGGCGATCGGTAATCATGAATTTGATAATCCGCTCACCGTATTACGCCAGCAGGAAAAGTGGGCCAAGTTCCCGTTGCTTTCCGCGAATATCTACCAGAAAAGTACTGGCGAGCGCCTGTTTAAACCGTGGGCGCTGTTTAAGCGTCAGGATCTGAAAATTG‑CCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAA‑TACTTCACTGATATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTACGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCGATGATCGTCGGTGGTCACTCGCAAGATCCGGTCT  >  NC_000913/505191‑505685

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