Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,045,2840A.62.5% 5.4 / 10.6 25coding (618/744 nt)ygfFputative NAD(P)‑binding oxidoreductase with NAD(P)‑binding Rossmann‑fold domain
Reads supporting (aligned to +/- strand):  ref base A (5/4);  new base . (0/15);  total (5/20)
Fisher's exact test for biased strand distribution p-value = 2.96e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTG‑TTTTATTT‑CCCGCCCGCCAAATCGAT‑AAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGC‑ATGGGGA‑TG‑TTCGACTTAACGCGATCGACGCGTCCAGGCTCG‑CCGCCGCTGGCGTGCATTTCGGTATAAATAAACCCTGGCCGCACGCAGTTAACGCGGATCCCCTGCGCGGCGACTTCCAGCGATAGTCCGGTGGTT‑AACGTATCAATCGCCCCTTTCGATGCCGCGTAATCAACATATTCCCCTGGCGAA‑CCCAACCGTGAGGCCACCGAAGAGACATTGACGATAGCGCCGCCACTGCCACCA‑‑TTTTTAAGCGCCA‑TGCGTTTTACCGCCTCGCGGCAGCAGAGAAAATAT‑CCCGTCACGTTGGTGGAAAGTACTCGGTTGATTCGCTCTGCGGTAAGGTTTTCAACGGTGCACTGGGTAAACAAGATCC  >  NC_000913/3045000‑3045639
                                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                                                                                                          
gtgAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACACCTT‑CCTG‑TTTTATTTCCCCGCCCGCCAAATCGAT‑AAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGC‑ATGGGGA‑TG‑TTCGACTT‑ACGCGATCGACGCGTCCAGGCTCG‑CCGCCGCTGGCGTGCATTTCGGTataaata                                                                                                                                                                                                                                                                                                                     <  1:510377/338‑1 (MQ=255)
                           gATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTTATTTCCCCGCCCGCCAAATCGATAAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGC‑ATGGGGA‑TG‑TTCGACTT‑ACGCGATCGACGCGTCCAGGCTCG‑CCGCCGCTGGCGTGCATTTCGGTATAAATAAACCCTGGCCGCACGCAGTTAACGCGGATCCCCTGCGCGGCGACTTCCAGCGATAGTCCggt                                                                                                                                                                                                                                                       <  1:510380/376‑1 (MQ=255)
                                              agaagaGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATC‑‑GTCAGACAACTTCCCTG‑TTTTATTT‑CCCGCCCGCCAAATCGAT‑AAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGC‑ATGGGGA‑TG‑TTCGACTTAACGCGATCGACGCGTCCAGGCTcg‑ccgccg                                                                                                                                                                                                                                                                                                                                             <  1:1390051/267‑1 (MQ=255)
                                                   agCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAG‑AGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTGTTTTTATTT‑CCCGCCCGCCAAATCGAT‑AAAACTTCCCGTGACG‑AAGAGGCTTTATCACTTAGTAGCCAGACAA‑GGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGC‑ATGGGGA‑TG‑TTCGACTT‑ACGCGATCGACGCGTCCAGGCTCG‑CCGCCGCTGGCGTGCATTTCGGTATAAATAAACCCTGGCCGCACGCAG‑TAACGCGGAt                                                                                                                                                                                                                                                                                        <  1:510382/313‑1 (MQ=255)
                                                                                                             aaGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTT‑CCTG‑‑TTTATTT‑CCCGCCCGCCAAATCGATAAAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGC‑ATGGGGA‑TG‑TTCGACTT‑ACGCGATCGACGCGTCCAGGCTCG‑CCGCCGCTGGCGTGCATTTCGGTATAAATAAACCCTGGCCGCACGCAGTTAACGCGGATCCCCTGCGCGGCGACTTCCAGCGATAGTc                                                                                                                                                                                                                                                           <  1:510386/286‑1 (MQ=255)
                                                                                                                                  ttgCATCCGGTCAGACAACTTCCCTG‑TTTTATTT‑CCCGCCCGCCAAATCGAT‑AAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGC‑ATGGGGA‑TG‑TTCGACTT‑ACGCGATCGACGCGTCCAGGCTCG‑CCGCCGCTGGCGTGCATTTCGGtata                                                                                                                                                                                                                                                                                                                         <  1:1390055/204‑1 (MQ=255)
                                                                                                                                                          tG‑TTTTATTT‑CCCGCCCGCCAAATCGAT‑AAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGC‑ATGGGGA‑TG‑TTCGACTTAACGCGATCGACGCGTCCAGGCTCG‑CCGCCGCTGGCGTGCATTTCGGTATAAATAAACCCTGGCCGCACGCAGTTAAcgc                                                                                                                                                                                                                                                                                            <  1:1390057/210‑1 (MQ=255)
                                                                                                                                                                  ttt‑CCCGCCCGCCAAATCGAT‑AAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAA‑GGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGC‑ATGGGGA‑TG‑TTCGACTT‑ACGCGATCGACGCGTCCAGGCTCG‑CCGCCGCTGGCGTGCATTTCGGTATAAATAAACCCTGGCCGCACGCAGTTAACGCGGATCCCCTGCGCGGCGACTTCCAGCGATAGTCCGGTGGTT‑AACGTATCAATCGCCCCTTTCGATGCCGCGTAATCAACATATTCCCCTGGCGAA‑CCCAACCGTGAGGCCACCGAAGAGACATTGACGATAGCGCCGCCACTGcca                                                                                                                                        <  1:510390/347‑1 (MQ=255)
                                                                                                                                                                           ccgccAAATCGAT‑AAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGC‑ATGGGGA‑TG‑TTCGACTTAACGCGATCGACGCGTCCAGGCTCG‑CCGCCGCTGGCGTGCATTTCGGTATAAATAAACCCTGGCCGCACGCAGTTAACGCGGATCCCCTGCgcggc                                                                                                                                                                                                                                                                            >  1:1390061/1‑211 (MQ=255)
                                                                                                                                                                             gccAAATCGAT‑AAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAA‑GGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGC‑ATGGGGA‑TG‑TTCGACTT‑ACGCGATCGACGCGTCCAGGCTCG‑CCGCCGCTGGCGTGCATTTCGGTATAAATAAACCCTGGCCGCACGCAGTTAACGCGGATCCCCTGCGCGGCGACTTCCAGCGATAGTCCGGTGGTT‑AACGTATCAATCGCCCCTTTCGATGCCGCGTAATCAACATATTCCCCTGGCGAA‑CCCAACCGTGAGGCCACCGAAGAGACATTGACGATAGCGCCGCCACTGCCACCATTTTTTTAAGCGCCA‑TGCGTTTTACCGCCTCGCGgcag                                                                                              <  1:510393/378‑1 (MQ=255)
                                                                                                                                                                                  aTCGAT‑AAAAC‑TCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGC‑ATGGGGA‑TG‑TTCGACTT‑ACGCGATCGACGCGTCCAGGCTCG‑CCGCCGCTGGCGTGCATTTCGGTATAAATAAACCCTGGCCGCACGCAGTTAACGCGGATCCCCTGCGCGGCGACTTCCAGCGATAGTCCGGTGGTT‑AACGTATCAATCGCCCCTTTCGATGCCGCGTAATCAACATATTCCCCTGGCGAA‑CCCAACCGt                                                                                                                                                                                  <  1:1390058/290‑1 (MQ=255)
                                                                                                                                                                                    cGAT‑AAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGC‑ATGGGGA‑TG‑TTCGACTTAACGCGATCGACGCGTCCAGGCTcg‑ccgccg                                                                                                                                                                                                                                                                                                                                             <  1:1390059/137‑1 (MQ=255)
                                                                                                                                                                                          aaaCTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGC‑ATGGGGA‑TG‑TTCGA‑TCTACGCGATCGACGCGTCCAGGCTCG‑CCGCCGCTGGCGTGCATTTCGGTATAAATAAACCCTGGCCGCACGCAGTTAACGCGGATCCCCTGCGCGGCGACTTCCAGCGATAGTCCGGTGGTT‑AACGTATCa                                                                                                                                                                                                                                         <  1:510395/230‑1 (MQ=255)
                                                                                                                                                                                           aaCTTCCCGTGACG‑AAGAGGCTTTATCAC‑TAGTAGCCAGACAA‑GGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGC‑ATGGGGA‑TGTTTCGACT‑‑ACGCGATCGACGCGTCCAGGCTCG‑CCGCCGCTGGCGTGCATTTCGGTATAAAT‑AACCCTGGCCGCACGCAGT‑‑ACGCGGATCCCCTGCGCGGCGACTTCCAGCGATAGTCCGGTGGTT‑AACGTATCAATCGCCCCTTTCGATGCCGCGTAATCAACATATTCCCCTGGCGAA‑CCCAACCGTGAGGCCACCGAAGAGACATTGACGATAGCGCCGCCACTGCCACCA‑‑TTTTTAAGCGCCA‑TGCGTTTTACCGCCTCGCGGCAGCAGAGAAAATATCCCCGTCACGtttgggtgga                                                              <  1:510394/390‑9 (MQ=255)
                                                                                                                                                                                                                     tCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGC‑ATGGGGA‑TG‑TTCGACTTAACGCGATCGACGCGTCCAGGCTCG‑CCGCCGCTGGCGTGCATTTCGGTATAAATAAACCCTGGCCGCACGCAGTTAACGCGGATCCCCTGCGCGGCGACTTCCAGCGATAGTCCGGTGGTT‑AACGTATCAATCGCCCCTTTCGATGCCGCGTAATCaa                                                                                                                                                                                                             >  1:510397/1‑232 (MQ=255)
                                                                                                                                                                                                                     tCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGC‑ATGGGGA‑TG‑TTCGACTTAACGCGATCGACGCGTCCAGGCTCG‑CCGCCGCTGGCGTGCATTTCGGTATAAATAAACCCTGGCCGCACGCAGTTAACGCGGATCCCCTGCGCGGCGACTTCCAGCGATAGTCCGGTGGTT‑AACGTATCAATCGCCCCTTTCGATGCCGCGTAATCAACATATTcccc                                                                                                                                                                                                   >  1:1390066/1‑242 (MQ=255)
                                                                                                                                                                                                                     tCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGC‑ATGGGGA‑TG‑TTCGACTT‑ACGCGATCGACGCGTCCAGGCTCG‑CCGCCGCTGGCGTGCATTTCGGTATAAATAAACCCTGGCCGCACGCAGTTAACGCGGATCCCCTGCGCGGCGACTTCCAGCGATAGTCCGGTGGTT‑AACGTATCAATCGCCCCTTTCGATGCCGCG‑AATCAACATATTCCCCTGGCGAA‑CCCAACCGTGAGGCCACCGAAGAGACATTGAc                                                                                                                                                           <  1:1390064/279‑1 (MQ=255)
                                                                                                                                                                                                                                    aCAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGC‑ATGGGGA‑TG‑TTCGACTT‑ACGCGATCGACGCGTCCAGGCTCG‑CCGCCGCTGGCGTGCATTTCGGTATAAATAAACCCTGGCCGCACGCAGTTAACGCGGATCCCCTGCGCGGCGACTTCCAGCGATAGTCCGGTGGTT‑AACGTATc                                                                                                                                                                                                                                          <  1:1390065/187‑1 (MQ=255)
                                                                                                                                                                                                                                             tGCGCGACCTCTTCTGCCTGTCCACCACGCTGCGATGGGGA‑TG‑TTCGACT‑‑ACGCGATCGACGCGTCCAGGCTCG‑CCGCCGCTGGCGTGCATTTCGGTATAAATAAACCCTGGCCGCACGCAGTTAACGCGGATCCCCTGCGCGGCGACTTCCAGCGATAGTCCGGTGGTTAAACGTATCAATCGCCCCTTTCGATGCCGCGTAATCAACATATTCCCCTGGCGAACCCCAACCGt                                                                                                                                                                                  <  1:510396/235‑1 (MQ=255)
                                                                                                                                                                                                                                               cgcgACCTCTTCTGCCTGTCCACCACGCTGC‑ATGGGGA‑TG‑TTCGACTT‑ACGCGATCGACGCGTCCAGGCTCG‑CCGCCGCTGGCGTGCATTTCGGTATAAATAAACCCTGGCCGCACGCAGTTAACGCGGATCCCCTGCGCGGCGACTTCCAGCGATAGTCCGGTGGTT‑AACGTATCAATCGCCCCTTTCGATGCCGCGTAATCAACATATTCCCCTGGCGAA‑CCCAACCGTGAGGCCACCGAAGAGACATTGACGATAGCGCCGCCACTGCCACCA‑‑TTTTTa                                                                                                                             <  1:1390067/282‑1 (MQ=255)
                                                                                                                                                                                                                                                                  ccaccaCGCTGCGATGGGGATTG‑TTCGACTT‑ACGCGATCGACGCGTCCAGGCTCG‑CCGCCGCTGGCGTGCATTTCGGTATAAAT‑AACCCTGGCCGCACGCAG‑TAACGCGGATCCCCTGCGCGGCGACTTCCAGCGATAGTCCGGTGGTT‑AACGTATCAATCGCCCCTTTCGATGCCGCGTAATCAACATATTCCCCTGGCGAA‑CCCAACCGTGAGGCCACCGAAGAGACATTGACGATAGCGCCGCCACTGCCACCA‑‑TTTTTAAGCGCCATTGCGTTTTACCGCCTCGCGGCAGCAGAGAAAATAT‑CCCGTCACGTTGGTGGAAAGTACTCGGTTGATTCGCTCTGCGGTAAGGTTTTCAACGGTGCACTGGGTAAACAAGATcc  <  1:510398/385‑1 (MQ=255)
                                                                                                                                                                                                                                                                       acGCTGC‑ATGGGGA‑TG‑TTCGACTTAACGCGATCGa                                                                                                                                                                                                                                                                                                                                                                  >  1:510400/1‑35 (MQ=255)
                                                                                                                                                                                                                                                                                  gggA‑TG‑TTCGACTTAACGCGATCGACGCGTCCAGGCTCG‑CCGCCGCTGGCGTGCATTTCGGTATAAATAAACCCTGGCCGCACGCAGTTAACGCGGATCCCCTGCGCGGCGACTTCCAGCGATAGTCCGGTGGTT‑AACGTATCAATCGCCCCTTTCGATGCCGCGTAATCaa                                                                                                                                                                                                             <  1:1390069/172‑1 (MQ=255)
                                                                                                                                                                                                                                                                                    gA‑TG‑TTCGACTT‑ACGCGATCGACGCGTCCAGGCTCG‑CCGCCGCTGGCGTGCATTTCGGTATAAATAAACCCTGGCCGCACGCAGTTAACGCGGATCCCCTGCGCGGCGACTTCCAGCGATAGTCCGGTGGTT‑AACGTATCAATCGCCCCTTTCGATGCCGCGTAATCAACATATTCCCCTGGCGAA‑CCCAACCGTGAGGCCACCGAAGAGACATTGACGATAGCGCCGCCACTGCCACCA‑‑TTTTTAAGCGCCA‑TGCGTTTTACCGCCTCGCGGCAGCAGAGaaaa                                                                                     <  1:1390070/285‑1 (MQ=255)
                                                                                                                                                                                                                                                                                            cGACTTAACGCGATCGACGCGTCCAGGCTCGCCCGCCGCTGGCGTGCATTTCGGTATAAATAAACCCTGGCCGCACGCAGTTAACGCGGATCCCCTg                                                                                                                                                                                                                                                                                  >  1:1390071/1‑97 (MQ=255)
                                                                                                                                                                                                                                                                                                  aaCGCGATCGACGCGTCCAGGCTCG‑CCGCCGCTGGCGTGCATTTCGGTATAAATAAACCCTGGCCGCACGCAGTTAACGCGGATcccc                                                                                                                                                                                                                                                                                    >  1:1390072/1‑88 (MQ=255)
                                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                                                                                                          
GTGAATAAACATGACGACGGTACTGGCGATATGTCGCGTTCACGTAAGAAGAGCTTTAACTGGTACAAAGAGGTGATTGCCAGCAACGGCGAGAAGCTTTAAGTCGATGAAGTACCGGATGCAATACTTGTTGCATCCGGTCAGACAACTTCCCTG‑TTTTATTT‑CCCGCCCGCCAAATCGAT‑AAAACTTCCCGTGACGTAAGAGGCTTTATCACTTAGTAGCCAGACAATGGCCTGCGCGACCTCTTCTGCCTGTCCACCACGCTGC‑ATGGGGA‑TG‑TTCGACTTAACGCGATCGACGCGTCCAGGCTCG‑CCGCCGCTGGCGTGCATTTCGGTATAAATAAACCCTGGCCGCACGCAGTTAACGCGGATCCCCTGCGCGGCGACTTCCAGCGATAGTCCGGTGGTT‑AACGTATCAATCGCCCCTTTCGATGCCGCGTAATCAACATATTCCCCTGGCGAA‑CCCAACCGTGAGGCCACCGAAGAGACATTGACGATAGCGCCGCCACTGCCACCA‑‑TTTTTAAGCGCCA‑TGCGTTTTACCGCCTCGCGGCAGCAGAGAAAATAT‑CCCGTCACGTTGGTGGAAAGTACTCGGTTGATTCGCTCTGCGGTAAGGTTTTCAACGGTGCACTGGGTAAACAAGATCC  >  NC_000913/3045000‑3045639

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