Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,741,2320CT37.8% 12.5 / 33.6 37intergenic (+469/+284)gltA/yeiWcitrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (19/4);  new base T (6/8);  total (25/12)
Fisher's exact test for biased strand distribution p-value = 2.74e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.77e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGAACCTCCGCCAACTGCAATTCGTGAACCAAGCAAAAGGCCCCGAAAAGGCCCTCCGCAGCCCCCCCAACCTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGC  >  NC_002947/4741185‑4741338
                                                 |                                                                                                            
ttGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgc                                  >  2:3518499/1‑123 (MQ=255)
   aaCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgcggcg                              >  1:3518414/1‑124 (MQ=255)
   aaCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgcg                                 <  2:3518359/121‑1 (MQ=255)
    aCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTGTGCGACTTCCTTTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGgcgc                             >  2:3518590/1‑124 (MQ=255)
    aCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGgcgc                             >  2:3518475/1‑124 (MQ=255)
      tGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTAGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgcggc                               >  2:3518473/1‑120 (MQ=255)
       gggTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCt                           <  2:3518360/123‑1 (MQ=255)
        ggTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCAGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCt                           <  2:3518361/122‑1 (MQ=255)
        ggTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGtcctcc                                                          <  2:3518363/91‑1 (MQ=255)
        ggTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCtt                          <  2:3518362/123‑1 (MQ=255)
          tataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGAt                       <  1:3518364/124‑1 (MQ=255)
           ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGAt                       >  2:3518681/1‑123 (MQ=255)
            taGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATc                      <  2:3518365/123‑1 (MQ=255)
            taGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCa                     >  2:3518525/1‑124 (MQ=255)
            taGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCa                     >  2:3518497/1‑124 (MQ=255)
              ggTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACg                   >  1:3518804/1‑124 (MQ=255)
              ggTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACg                   >  1:3518582/1‑124 (MQ=255)
                tCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCa                     <  2:3518367/120‑1 (MQ=255)
                tCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAg                                   >  2:3518565/1‑106 (MQ=255)
                  cAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTca               >  1:3518544/1‑124 (MQ=255)
                  cAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTca               >  2:3518649/1‑124 (MQ=255)
                   aGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTca               >  1:3518529/1‑123 (MQ=255)
                    gCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGCTCACg                   >  1:3518588/1‑118 (MQ=255)
                     ctctGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCacac            <  2:3518374/124‑1 (MQ=255)
                     ctctGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTcaca             <  2:3518390/123‑1 (MQ=255)
                       ctGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcg          >  1:3518619/1‑124 (MQ=255)
                        tGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc         >  2:3518566/1‑124 (MQ=255)
                        tGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcg          <  2:3518414/123‑1 (MQ=255)
                         ggTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc         <  2:3518416/123‑1 (MQ=255)
                           tCACCTTCGTGCTTGAGCAGCTTCTACGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgcg        >  2:3518501/1‑122 (MQ=255)
                             aCCTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGa       >  1:3518656/1‑121 (MQ=255)
                             aCCTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAg    >  1:3518609/1‑124 (MQ=255)
                              ccTTCGTGCTTGAGCAGCTTCTTCGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAgg   >  2:3518615/1‑124 (MQ=255)
                              ccTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAgg   >  1:3518571/1‑124 (MQ=255)
                              ccTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAgg   >  2:3518587/1‑124 (MQ=255)
                              ccTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAgg   >  2:3518600/1‑124 (MQ=255)
                                ttCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGc  >  2:3518560/1‑123 (MQ=255)
                                                 |                                                                                                            
TTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGAACCTCCGCCAACTGCAATTCGTGAACCAAGCAAAAGGCCCCGAAAAGGCCCTCCGCAGCCCCCCCAACCTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGC  >  NC_002947/4741185‑4741338

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