Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077792,225,7000TG42.2% 64.1 / 16.7 45V61G (GTG→GGG) dldD‑lactate dehydrogenase, FAD‑binding, NADH independent
Reads supporting (aligned to +/- strand):  ref base T (2/24);  new base G (18/1);  total (20/25)
Fisher's exact test for biased strand distribution p-value = 1.96e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CAAAAACGGCCCGCTATCGCAAGGGCTTCCGTTCTGGTCAGGGCGACGCGCTGGCTGTCGTTTTCCCTGGCTCACTACTAGAATTGTGGCGGGTGCTGAAAGCCTGCGTCACCGCCGACAAAATTATTCTGATGCAGGCCGCCAATACAGGCCTGACCGAAGGATCGACGCCAAACGGTAACGATTATGATC  >  NC_007779/2225607‑2225798
                                                                                              |                                                                                                  
cAAAAACGGCCCGCTATCGCAAGGGCTTCCGTTCTGGTCAGGGCGACGCGCTGGCTGTCGTTTTCCCTGGCTCACTACTAGAATTGTGGCGGGGGCTGgaa                                                                                             >  3:804380/1‑101 (MQ=255)
 aaaaaCGGCCCGCTATCGCAAGGGCTTCCGTTCTGGTCAGGGCGACGCGCTGGCTGTCGTTTTCCCTGGCTCACTACTAGAATTGTGGCGGGTGCTGAAAg                                                                                            <  3:667235/101‑1 (MQ=255)
  aaaaCGGCCCGCTATCGCAAGGGCTTCCGTTCTGGTCAGGGCGACGCGCTGGCTGTCGTTTTCCCTGGCTCACTACTAGAATTGTGGCGGGTGCTGAAAGc                                                                                           <  1:1283238/101‑1 (MQ=255)
            gCTATCGCAAGGGCTTCCGTTCTGGTCAGGGCGACGCGCTGGCTGTCGTTTTCCCTGGCTCACTACTAGAATTGTGGCGGGTGCTGAAAGCCTGCGTCAcc                                                                                 <  3:1265109/101‑1 (MQ=255)
              tATCGCAAGGGCTTCCGTTCTGGTCAGGGCGACGCGCTGGCTGTCGTTTTCCCTGGCTCACTACTAGAATTGTGGCGGGTTGTGTAAACACGCGCCAccgc                                                                               >  3:294865/1‑101 (MQ=255)
              tATCGCAAGGGCTTCCGTTCTGGTCAGGGCGACGCGCTGGCTGTCGTTTTCCCTGGCTCACTACTAGAATTGTGGCGGGTGCTGAAAGCCTGCGTCAccgc                                                                               <  1:774172/101‑1 (MQ=255)
               aTCGCAAGGGCTTCCGTTCTGGTCAGGGCGACGCGCTGGCTGTCGTTTTCCCTGGCTCACTACTAGAATTGTGGCGGGTGCTGAAAGCCTGCGCCAccgcc                                                                              <  1:1255735/101‑1 (MQ=255)
                      gggCTTCCGTTCTGGTCAGGGCGACGCGCTGGCTGTCGTTTTCCCTGGCTCACTACTAGAATTGTGGCGGGTGCTGAAAGCCTGCGTCACCGCCGACaaaa                                                                       <  1:1097115/101‑1 (MQ=255)
                          ttCCGTTCTGGTCAGGGCGACGCGCTGGCTGTCGTTTTCCCTGGCTCACTACTAGAATTGTGGCGGGGGCTGAAAAACCGCGGCACCGCCCCCAAAaatat                                                                   >  1:163039/1‑101 (MQ=255)
                          ttCCGTTCTGGTCAGGGCGACGCGCTGGCTGTCGTTTTCCCTGGCTCACTACTAGAATTGTGGCGGGGGCGGGAAGCCCGCGTCACCGCCGACAAAaatat                                                                   >  1:1044516/1‑101 (MQ=255)
                            ccGTTCTGGTCAGGGCGACGCGCTGGCTGTCGTTTTCCCTGGCTCACTACTAGAATTGTGGCGGGGGGTGAAAGCCTGCGTCACCCCCGACAAAaatattt                                                                 >  3:953004/1‑100 (MQ=255)
                            ccGTTCTGGTCAGGGCGACGCGCTGGCTGTCGTTTTCCCTGGCTCACTACTAGAATTGTGGCGGGGGCTGGAAGCCTGCGTCACCCCCCGCAAAaatattt                                                                 >  3:24770/1‑100 (MQ=255)
                            ccGTTCTGGTCAGGGCGACGCGCTGGCTGTCGTTTTCCCTGGCTCACTACTAGAATTGTGGCGGGGGCTGGAAAACCGCGGCCCCCCCCACAAAaatattt                                                                 >  3:332623/1‑100 (MQ=255)
                            ccGTTCTGGTCAGGGCGACGCGCTGGCTGTCGTTTTCCCTGGCTCACTACTAGAATTGTGGCGGGGGCTGAAAGCCCGCGGCGCCCCCGACAAAaatattt                                                                 >  1:403592/1‑100 (MQ=255)
                             cGTTCTGGTCAGGGCGACGCGCTGGCTGTCGTTTTCCCTGGCTCACTACTAGAATTGTGGCGGGGGCTGAAAGACCGCGTCACCGCCCCCAAAAATATTCt                                                                >  1:675998/1‑101 (MQ=255)
                              gTTCTGGTCAGGGCGACGCGCTGGCTGTCGTTTTCCCTGGCTCACTACTAGAATTGTGGCGGGGGGGGGAAAGCCGCGTCACCCCCGACAAAAATATTCtt                                                               >  3:946730/1‑100 (MQ=255)
                               ttCTGGTCAGGGCGACGCGCTGGCTGTCGTTTTCCCTGGCTCACTACTAGAATTGTGGCGGGTGCTGAAAGCCTGCGTCACCGCCGACAAAATTATTCTGa                                                              <  3:337875/101‑1 (MQ=255)
                                tCTGGTCAGGGCGACGCGCTGGCTGTCGTTTTCCCTGGCTCACTACTAGAATTGTGGCGGGGGCTGAAAGCCCGCGGCACCCCCGGCAAAattaattttat                                                             >  1:732258/1‑94 (MQ=255)
                                 cTGGTCAGGGCGACGCGCTGGCTGTCGTTTTCCCTGGCTCACTACTAGAATTGTGGCGGGGGCTGAAAGACCGCGGCACCGCCGCCAAAAATATTCTGATg                                                            >  4:30394/1‑101 (MQ=255)
                                     tCAGGGCGACGCGCTGGCTGTCGTTTTCCCTGGCTCACTACTAGAATTGTGGCGGGTGCTGAAAGCCTGCGTCACCGCCGACAAAATTATTCTGATGCAgg                                                        <  1:712614/101‑1 (MQ=255)
                                     tCAGGGCGACGCGCTGGCTGTCGTTTTCCCTGGCTCACTACTAGAATTGTGGCGGGTGCTGAAAGCCTGCGTCACCGCCGACAAAATTATTCTGATGCAgg                                                        <  1:208036/101‑1 (MQ=255)
                                     tCAGGGCGACGCGCTGGCTGTCGTTTTCCCTGGCTCACTACTAGAATTGTGGCGGGGGGTGAAAACCCGCGTCACCCCCCCCAAAAATATTCTTATGCAgg                                                        >  3:333775/1‑101 (MQ=255)
                                           cGACGCGCTGGCTGTCGTTTTCCCTGGCTCACTACTAGAATTGTGGCGGGTGCTGAAAGCCTGCGTCACCGCCGACAAAATTATTCTGATGCAGGCCGCCa                                                  <  3:93749/101‑1 (MQ=255)
                                               gcgcTGGCTGTCGTTTTCCCTGGCTCACTACTAGAATTGTGGCGGGGGCCGCAAGACTCGGTCACCCCCAACCAAATTTTTCTTCAGAAGGCCCCCACTAc                                              >  3:216474/1‑101 (MQ=255)
                                                    ggctgTCGTTTTCCCTGGCTCACTACTAGAATTGTGGCGGGGGGTGAAAACCTGCGGCACCGCCGACAAAAATATTCTTATGCAGGCCGCCacaaaaagcc                                         >  1:575416/1‑92 (MQ=255)
                                                          cGTTTTCCCTGGCTCACTACTAGAATTGTGGCGGGTGCTGAAAGCCTGCGTCACCGCCGACAAAATTATTCTGATGCAGGCCGCCAATACAGGCCTGACCg                                   <  3:108170/101‑1 (MQ=255)
                                                            ttttCCCTGGCTCACTACTAGAATTGTGGCGGGGGCTGAAAAACCGCGTCACCGCCCACAAAAATATTCTGATGCAGGCCGCCAATACAGGCCCGACCGaa                                 >  1:556481/1‑101 (MQ=255)
                                                               tCCCTGGCTCACTACTAGAATTGTGGCGGGTGCTGAAAGCCTGCGTCACCGCCGACAAAATTATTCTGATGCAGGCCGCCAATACAGGCCTGACCGAAGGa                              <  1:653226/101‑1 (MQ=255)
                                                               tCCCTGGCTCACTACTAGAATTGTGGCGGGTGCTGAAAGCCTGCGTCACCGCCGACAAAATTATTCTGATGCAGGCCGCCAATACAGGCCTGACCGAAGGa                              <  1:1503321/101‑1 (MQ=255)
                                                               tCCCTCGCTCACTACTAGAATTGTGGCGGGTGCTGAAAGCCTGCGTCACCGCCGACAAAATTATTCTGATGCAGGCCGCCAATACAGGCCTGACCGAAGGa                              <  1:1038300/101‑1 (MQ=255)
                                                                    ggCTCACTACTAGAATTGTGGCGGGTGCTGAAAGCCTGCGTCACCGCCGACAAAATTATTCTGATGCAGGCCGCCAATACAGGCCTGACCGAAGGATCGAc                         <  2:6633/101‑1 (MQ=255)
                                                                     gCTCACTACTAGAATTGTGGCGGGTGCTGAAAGCCTGCGTCACCGCCGACAAAATTATTCTGATGCGGGCCGCCAATACAGGCCTGACCGAAGGATCGCCg                        <  1:373150/101‑1 (MQ=255)
                                                                     gCTCACTACTAGAATTGTGGCGGGTGCTGAAAGCCTGCGTCACCGCCGACAAAATTATTCTGATGCAGGCCGCCAATACAGGCCTGACCGAAGGATCGACg                        <  1:330364/101‑1 (MQ=255)
                                                                     gCTCACTACTAGAATTGTGGCGGGTGCTGAAAGCCTGCGTCACCGCCGACAAAATTATTCTGATGCAGGCCGCCAATACAGGCCTGACCGAAGGATCGACg                        <  1:24472/101‑1 (MQ=255)
                                                                      cTCACTACTAGAATTGTGGCGGGTGCTGAAAGCCTGCGTCACCGCCGACAAAATTATTCTGATGCAGGCCGCCAATACAGGCCTGACCGAAGGATCGACGc                       <  3:1507373/101‑1 (MQ=255)
                                                                      cTCACTACTAGAATTGTGGCGGGTGCTGAAAGCCTGCGTCACCGCCGACAAAATTATTCTGATGCAGGCCGCCAATACAGGCCTGACCGAAGGATCGACGc                       <  3:1400739/101‑1 (MQ=255)
                                                                         actactCGCATTGTGGAGGGGGGTGAAAGCCCGCGGCACCGCCGACAAAATTATTCTGATGCAGGCCGCCAATACAGGCCTGACCGAAGGATCGACGCCaa                    >  4:128570/1‑101 (MQ=255)
                                                                         actactAGAATTGTGGCGGGGGCTGAAAGCCCGCGGCACCGCCGACAAAATTATTCTGATGCAGGCCGCCAATACAGGCCTGACCGAAGGATCGACGCCaa                    >  4:47524/1‑101 (MQ=255)
                                                                           tactaGAATTGTGGCGGGTGCTGAAAGCCTGCGTCACCGCCGACAAAATTATTCTGATGCAGGCCGCCAATACAGGCCTGACCGAAGGATCGACGCCAAAc                  <  3:682710/101‑1 (MQ=255)
                                                                            actaGAATTGTGGCGGGTGCTGAAAGCCTGCGTCACCGCCGACAAAATTATTCTGATGCAGGCCGCCAATACAGGCCTGACCGAAGGATCGACGCCAAACg                 <  3:958486/101‑1 (MQ=255)
                                                                            actaGAATTGTGGCGGGTGCTGAAAGCCTGCGTCACCGCCGACAAAATTATTCTGATGCAGGCCGCCAATACAGGCCTGACCGAAGGATCGACGCCAAACg                 <  2:30413/101‑1 (MQ=255)
                                                                               aGAATTGTGGCGGGGGCTGAAAACCTGCGTCACCGCCGACAAAATTATTCTGATGCAGGCCGCCAATACAGGCCTGACCCAAGGATCGACGCCCAACGGTa              >  3:436807/1‑101 (MQ=255)
                                                                                     gtgtCGGGGGCTGAAAGCCGGCGTCACCGCCGACAAAATTATTCTGATGCAGGCCGCAAATACAGGCCTGACCGAAGGATCGACGCCAAACGGTAACGAtt        <  1:1373464/101‑1 (MQ=255)
                                                                                           gggTGCTGAAAGCCTGCGTCACCGCCGACAAAATTATTCTGATGCAGGCCGCCAATACAGGCCTGACCGAAGGATCGACGCCAAACGGTAACGATTATGAt   >  1:569276/1‑101 (MQ=255)
                                                                                            ggTGCTGAAAGCCTGCGTCACCGCCGACAAAATTATTCTGATGCAGGCCGCCAATACAGGCCTGACCGAAGGATCGACGCCAAACGGTAACGATTATGATc  <  1:937351/101‑1 (MQ=255)
                                                                                              |                                                                                                  
CAAAAACGGCCCGCTATCGCAAGGGCTTCCGTTCTGGTCAGGGCGACGCGCTGGCTGTCGTTTTCCCTGGCTCACTACTAGAATTGTGGCGGGTGCTGAAAGCCTGCGTCACCGCCGACAAAATTATTCTGATGCAGGCCGCCAATACAGGCCTGACCGAAGGATCGACGCCAAACGGTAACGATTATGATC  >  NC_007779/2225607‑2225798

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: