Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,842,9140AG46.9% 23.5 / 10.4 32M182V (ATG→GTG) gdhAglutamate dehydrogenase, NADP‑specific
Reads supporting (aligned to +/- strand):  ref base A (1/16);  new base G (15/0);  total (16/16)
Fisher's exact test for biased strand distribution p-value = 5.66e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGAT‑GATGAAAAAGCTCTCCAACAATACCGCCTGC‑GTC‑TTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACGCCACGGTATGGGTTTTGA  >  NC_000913/1842764‑1843062
                                                                                                                                                      |                                                                                                                                                       
gATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGga                                                                                                                                                         <  1:476511/151‑1 (MQ=255)
gATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTCTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGga                                                                                                                                                         <  1:182388/151‑1 (MQ=255)
   ttCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGgat‑g                                                                                                                                                      <  2:490421/150‑1 (MQ=255)
          cGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCAACTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGA‑GATGaaaa                                                                                                                                               >  1:611617/1‑150 (MQ=255)
              aGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGT‑GACGaaaaa                                                                                                                                              >  1:486405/1‑147 (MQ=255)
                  aaaaGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGT‑GTTGgaaaaaagctcc                                                                                                                                       >  2:387952/1‑143 (MQ=255)
                  aaaaGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGT‑GGTGGAAAAAATCTcc                                                                                                                                       >  1:371767/1‑150 (MQ=255)
                  aaaaGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGAT‑GATGAAAAAGCTCTCca                                                                                                                                      <  1:288101/151‑1 (MQ=255)
                        gAAGGTGAAGTGCTGCGATTGTTCCAGGCGCTGATGACTGAAATGTATCGCCACCTGGGCGCGGATACCGACGTTTCGGCAGGTGATATCGAGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGAT‑GATGAAAAAGCTATCcaacaaga                                                                                                                                <  2:633279/151‑3 (MQ=255)
                        gAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGG‑GATGGAAAAAATTTCcaaaca                                                                                                                                  >  2:102121/1‑146 (MQ=255)
                         aaGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGG‑GATGAAAAAGCTCTCcaacaaaa                                                                                                                                >  1:516841/1‑148 (MQ=255)
                         aaGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGAT‑GATGAAAAAGCTCTCCAACAATAc                                                                                                                               <  1:521713/151‑1 (MQ=255)
                                                         atgaCTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGGTTGAAAAAAACTACCAAAAAAACCCCCTGG‑GTC‑TTACCCCGGAGGGGCCttcat                                                                                              >  2:333053/1‑148 (MQ=255)
                                                          tgaCTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGAT‑GATGAAAAAGCTCTCCAACAATACCGCCTGC‑GTC‑TTCACCGGTAAGGGCCTTTCAt                                                                                             <  1:580087/150‑1 (MQ=255)
                                                           gaCTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGAT‑GATGAAAAAGCTCTCCAACAATACCGCCTGC‑GTC‑TTCACCGGTAAGGGCCTTTCAttt                                                                                           <  2:334992/151‑1 (MQ=255)
                                                               gAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGAT‑GATGAAAAAGCTCTCCAACAATACCGCCTGC‑GTC‑TTCACCGGTAAGGGCCTTTCATTTggcg                                                                                       <  2:599595/151‑1 (MQ=255)
                                                                aaCTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGTTCGTGAAGTCGGCTTTATGGCGGGGAT‑GATGAAAAAGCTCTCCAACAATACCGCCTGC‑GTC‑TTCACCGGTAAGGGCCTTTCATTTggcg                                                                                       <  1:594772/150‑1 (MQ=255)
                                                                 aCTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGAG‑GATTAAAAAAACCTCCACCAAAACCGCCCGC‑GGC‑TTCACCGGGAAGGGCCCTTTAtttttcggg                                                                                     >  1:256389/1‑145 (MQ=255)
                                                                       tCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGT‑GATGAAAAAACTCTCCAACAATACCGCCCGC‑GTC‑TTCACCGGGAAGGGGCTTTCCTTTTGCGGCCGTCtt                                                                               >  2:594866/1‑151 (MQ=255)
                                                                        cGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGAT‑GATGAAAAAGCTCTCCAACAATACCGCCTGC‑GTC‑TTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCtt                                                                               <  2:50697/150‑1 (MQ=255)
                                                                               tGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGAT‑GATGAAAAAGCTCTCCAACAATACCGCCTGC‑GTC‑TTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCTTATTCg                                                                          <  1:387082/148‑1 (MQ=255)
                                                                                    gcgGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGAT‑GATGAAAAAGCTCTCCAACAATACCGCCTGC‑GTC‑TTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGaa                                                                   <  1:557078/150‑1 (MQ=255)
                                                                                           ccGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGG‑GGTGAAAAAGAACTCCAACAAAACCCCCCGC‑CTC‑TTTCCCCCGAAGGGCGTTTCATTTTTCGGCCGGCTTTTTTCCCCCGCAGCAAcc                                                             >  2:80216/1‑149 (MQ=255)
                                                                                            cGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGT‑GATGAAAAAACTCTCCAAAAAAAACGCCCGC‑GTC‑TTTACCCGGAAGGGGCTTTTTTTTTGGGGCCGGCTTTTTTTCCCCGGAGCTAcccg                                                           >  1:40663/1‑148 (MQ=255)
                                                                                                   ccGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGAT‑GATGAAAAAGCTCTCCAACAATACCGCCTGC‑GTC‑TTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGt                                                    <  1:57782/150‑1 (MQ=255)
                                                                                                     ggCAGGTGATATCGGGGTGGGTGGTCGTGAAGCCGGCTTTAGGGCGGGGAT‑GATGAAAAACCTCCCCAACAATTCCGCCTGC‑GTC‑TTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTATgg                                                     <  2:448700/147‑1 (MQ=255)
                                                                                                          gTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGT‑GGTGGAAAAAATCTCCCACAAAACCCCCCGC‑GTC‑TTCTCCCGTAAGGGCCTTTCTTTTTGCGGCGGTCTTATTTTCCCCGAAGCTACCGGCCACGGGCTCGttt                                             >  1:62184/1‑150 (MQ=255)
                                                                                                            gATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGT‑GATGAAAAAACTCTCCAACAAAAACCCCTTC‑GTC‑TTTACCCGTTAGGGCCTTTTCTTTTGCGGGCGTCTTTTTCTCCCCGCAGCCACCCGCCACCGTCTGGGTtttt                                          >  1:639915/1‑151 (MQ=255)
                                                                                                             atatCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGT‑GATGAAAAAAATCTCCCACAAAACCGCCTGC‑GGC‑TTCACCCGGAAGGGGCCTTTATTTTGCGGGCGGCTTTTTTTCCCCGGAGCTACCCGCCACGGGCTGGGtttttt                                         >  1:501881/1‑151 (MQ=255)
                                                                                                                       tggtggTCGTGAAGTCGGCTTTATGGCGGGGGG‑GATGGAAAAAATCTCCCACAAAAACCCCCGC‑GTCTTTAACCGGAAAGGGCCTTT‑TTTTTGCGGCAGGCTTTTTTTCCCCGAAGGCACCGGCCACGGTCTTGTttttttttccaaa                                  >  2:507411/1‑141 (MQ=255)
                                                                                                                           ggtCGTGAAGTCGGCTTTATGGCGGGGAT‑GATGAAAAAGCTCTCCAACAATACCGCCTGC‑GTC‑TTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGc                           <  1:458764/151‑1 (MQ=255)
                                                                                                                            gtCGTGAAGTCGGCTTTATGGCGGGGGG‑GGTGAAAAAAATCTCCCACAAAACCCCCCGCGGCC‑TCCCCCGGAAAGGGCCTTTC‑TTTTGCGGCCGTCTTATTTCCCCCGGAGCTACCGGGTACCGGCTGGGTTATTTTTCAAAAGCaa                              >  1:42242/1‑147 (MQ=255)
                                                                                                                                    gTCGGCTTTATGGCGGGGAT‑GATGAAAAAGCTCTCCAACAATACCGCCTGC‑GTC‑TTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACGcc                  <  1:440292/151‑1 (MQ=255)
                                                                                                                                      cGGCTTTATGGCGGGGAT‑GATGAAAAAGCTCTCCAACAATACCGCCTGC‑GTC‑TTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACGCCAc                <  2:607016/151‑1 (MQ=255)
                                                                                                                                             aTGGCGGGGAT‑GATGAAAAAGCTCTCCAACAATACCGCCTGC‑GTC‑TTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACGCCACGGTATg          <  2:469522/150‑1 (MQ=255)
                                                                                                                                                      at‑gatgAAAAAGCTCTCCAACAATACCGCCTGC‑GTC‑TTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACGCCACGGTATGGGTTTTGa  >  1:84006/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                                       
GATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGAT‑GATGAAAAAGCTCTCCAACAATACCGCCTGC‑GTC‑TTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACGCCACGGTATGGGTTTTGA  >  NC_000913/1842764‑1843062

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