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

GATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACGCCACGGTATGGGTTTTGA  >  NC_000913/1842764‑1843062
                                                                                                                                                      |                                                                                                                                                       
gATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGga                                                                                                                                                         <  1:476511/151‑1 (MQ=255)
gATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTCTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGga                                                                                                                                                         <  1:182388/151‑1 (MQ=255)
   ttCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGgatg                                                                                                                                                      <  2:490421/150‑1 (MQ=255)
          cGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCAACTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGAGATGaaaa                                                                                                                                               >  1:611617/1‑150 (MQ=255)
              aGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGACGaaaaa                                                                                                                                              >  1:486405/1‑147 (MQ=255)
                  aaaaGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGTTGgaaaaaagctcc                                                                                                                                       >  2:387952/1‑143 (MQ=255)
                  aaaaGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGGTGGAAAAAATCTcc                                                                                                                                       >  1:371767/1‑150 (MQ=255)
                  aaaaGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCca                                                                                                                                      <  1:288101/151‑1 (MQ=255)
                        gAAGGTGAAGTGCTGCGATTGTTCCAGGCGCTGATGACTGAAATGTATCGCCACCTGGGCGCGGATACCGACGTTTCGGCAGGTGATATCGAGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTATCcaacaaga                                                                                                                                <  2:633279/151‑3 (MQ=255)
                        gAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGATGGAAAAAATTTCcaaaca                                                                                                                                  >  2:102121/1‑146 (MQ=255)
                         aaGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGATGAAAAAGCTCTCcaacaaaa                                                                                                                                >  1:516841/1‑148 (MQ=255)
                         aaGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATAc                                                                                                                               <  1:521713/151‑1 (MQ=255)
                                                         atgaCTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGGTTGAAAAAAACTACCAAAAAAACCCCCTGGGTCTTACCCCGGAGGGGCCttcat                                                                                              >  2:333053/1‑148 (MQ=255)
                                                          tgaCTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGTAAGGGCCTTTCAt                                                                                             <  1:580087/150‑1 (MQ=255)
                                                           gaCTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGTAAGGGCCTTTCAttt                                                                                           <  2:334992/151‑1 (MQ=255)
                                                               gAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGTAAGGGCCTTTCATTTggcg                                                                                       <  2:599595/151‑1 (MQ=255)
                                                                aaCTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGTTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGTAAGGGCCTTTCATTTggcg                                                                                       <  1:594772/150‑1 (MQ=255)
                                                                 aCTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGAGGATTAAAAAAACCTCCACCAAAACCGCCCGCGGCTTCACCGGGAAGGGCCCTTTAtttttcggg                                                                                     >  1:256389/1‑145 (MQ=255)
                                                                       tCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGATGAAAAAACTCTCCAACAATACCGCCCGCGTCTTCACCGGGAAGGGGCTTTCCTTTTGCGGCCGTCtt                                                                               >  2:594866/1‑151 (MQ=255)
                                                                        cGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCtt                                                                               <  2:50697/150‑1 (MQ=255)
                                                                               tGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCTTATTCg                                                                          <  1:387082/148‑1 (MQ=255)
                                                                                    gcgGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGaa                                                                   <  1:557078/150‑1 (MQ=255)
                                                                                           ccGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGGTGAAAAAGAACTCCAACAAAACCCCCCGCCTCTTTCCCCCGAAGGGCGTTTCATTTTTCGGCCGGCTTTTTTCCCCCGCAGCAAcc                                                             >  2:80216/1‑149 (MQ=255)
                                                                                            cGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGATGAAAAAACTCTCCAAAAAAAACGCCCGCGTCTTTACCCGGAAGGGGCTTTTTTTTTGGGGCCGGCTTTTTTTCCCCGGAGCTAcccg                                                           >  1:40663/1‑148 (MQ=255)
                                                                                                   ccGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGt                                                    <  1:57782/150‑1 (MQ=255)
                                                                                                     ggCAGGTGATATCGGGGTGGGTGGTCGTGAAGCCGGCTTTAGGGCGGGGATGATGAAAAACCTCCCCAACAATTCCGCCTGCGTCTTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTATgg                                                     <  2:448700/147‑1 (MQ=255)
                                                                                                          gTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGGTGGAAAAAATCTCCCACAAAACCCCCCGCGTCTTCTCCCGTAAGGGCCTTTCTTTTTGCGGCGGTCTTATTTTCCCCGAAGCTACCGGCCACGGGCTCGttt                                             >  1:62184/1‑150 (MQ=255)
                                                                                                            gATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGATGAAAAAACTCTCCAACAAAAACCCCTTCGTCTTTACCCGTTAGGGCCTTTTCTTTTGCGGGCGTCTTTTTCTCCCCGCAGCCACCCGCCACCGTCTGGGTtttt                                          >  1:639915/1‑151 (MQ=255)
                                                                                                             atatCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGATGAAAAAAATCTCCCACAAAACCGCCTGCGGCTTCACCCGGAAGGGGCCTTTATTTTGCGGGCGGCTTTTTTTCCCCGGAGCTACCCGCCACGGGCTGGGtttttt                                         >  1:501881/1‑151 (MQ=255)
                                                                                                                       tggtggTCGTGAAGTCGGCTTTATGGCGGGGGGGATGGAAAAAATCTCCCACAAAAACCCCCGCGTCTTTAACCGGAAAGGGCCTTTTTTTTGCGGCAGGCTTTTTTTCCCCGAAGGCACCGGCCACGGTCTTGTttttttttccaaa                                  >  2:507411/1‑141 (MQ=255)
                                                                                                                           ggtCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGc                           <  1:458764/151‑1 (MQ=255)
                                                                                                                            gtCGTGAAGTCGGCTTTATGGCGGGGGGGGTGAAAAAAATCTCCCACAAAACCCCCCGCGGCCTCCCCCGGAAAGGGCCTTTCTTTTGCGGCCGTCTTATTTCCCCCGGAGCTACCGGGTACCGGCTGGGTTATTTTTCAAAAGCaa                              >  1:42242/1‑147 (MQ=255)
                                                                                                                                    gTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACGcc                  <  1:440292/151‑1 (MQ=255)
                                                                                                                                      cGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACGCCAc                <  2:607016/151‑1 (MQ=255)
                                                                                                                                             aTGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACGCCACGGTATg          <  2:469522/150‑1 (MQ=255)
                                                                                                                                                      atgatgAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACGCCACGGTATGGGTTTTGa  >  1:84006/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                                       
GATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACGCCACGGTATGGGTTTTGA  >  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: