Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,101,5140GT28.3% 106.5 / 15.9 46intergenic (+84/‑176)sodA/kdgTsuperoxide dismutase (Mn)/2‑dehydro‑3‑deoxy‑D‑gluconate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base G (4/29);  new base T (13/0);  total (17/29)
Fisher's exact test for biased strand distribution p-value = 2.34e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CTGGAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGGAATTA  >  NC_000913/4101376‑4101651
                                                                                                                                          |                                                                                                                                         
cTGGAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTtgt                                                                                                                                          >  1:238379/1‑140 (MQ=255)
 tGGAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTtgtg                                                                                                                                         <  1:303627/140‑1 (MQ=255)
        tggtgAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTtgtg                                                                                                                                         >  2:990883/1‑133 (MQ=255)
        tggtgAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTtgtg                                                                                                                                         <  1:990883/133‑1 (MQ=255)
           tgAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATGAAtttta                                                                                                                               >  2:434580/1‑138 (MQ=255)
             aaCTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTCaaa                                                                                                                             >  2:206379/1‑140 (MQ=255)
                  ggACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTTCGATCAATTTCAAAATaaa                                                                                                                        >  1:771949/1‑140 (MQ=255)
                  ggACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATaaa                                                                                                                        >  1:1255573/1‑140 (MQ=255)
                       aaGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTTTGATTAATTTTAAAATaaaaaaaa                                                                                                                   >  2:870872/1‑140 (MQ=255)
                       aaGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATTAATTTCAaaaaaaaaaaaa                                                                                                                   >  1:745387/1‑140 (MQ=255)
                            gcGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGGTTAATTTTAAAATaaaaaaaaaaatc                                                                                                              >  2:1097288/1‑135 (MQ=255)
                            gcGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTCAaaaaaaaaaaaaaaaac                                                                                                              >  2:727638/1‑135 (MQ=255)
                                cACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGat                                                                                                          <  2:457482/140‑1 (MQ=255)
                                     tttGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGataaaaa                                                                                                     <  1:1347268/140‑1 (MQ=255)
                                           gcgAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAAc                                                                                               <  2:978614/140‑1 (MQ=255)
                                           gcgAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAAc                                                                                               <  1:713495/140‑1 (MQ=255)
                                             gAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTTAAAATAAAAAAAATGATCCGataaaaaaaaaaaaa                                                                                             >  2:50961/1‑137 (MQ=255)
                                                      atCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAAt                                                                                    <  1:782421/140‑1 (MQ=255)
                                                      atCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTCGTCAATTTTAAAATAAAAAAAAAGATCCGATAAAAAAAAAAAAACGTTTTAAt                                                                                    >  2:1334883/1‑140 (MQ=255)
                                                           ttGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGAtg                                                                               <  1:1449954/140‑1 (MQ=255)
                                                           ttGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGAtg                                                                               <  1:1492386/140‑1 (MQ=255)
                                                            tGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTTAAAAAAAAAACAAAGAACCGATAAAAAAAAAAAAACGTTTCAAATTAtgg                                                                              >  2:1492386/1‑139 (MQ=255)
                                                             gccgccTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGAtgtg                                                                             <  1:1355442/140‑1 (MQ=255)
                                                                     ctgcAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGAt                                                                     <  1:1097288/140‑1 (MQ=255)
                                                                     ctgcAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGAt                                                                     <  1:434580/140‑1 (MQ=255)
                                                                        cAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGCGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCAc                                                                  <  2:771949/140‑1 (MQ=255)
                                                                             aGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTa                                                             <  1:350032/140‑1 (MQ=255)
                                                                             aGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTa                                                             <  1:211976/140‑1 (MQ=255)
                                                                                cGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTTTGATCAATTTTAAAAAAAAAAAAATGATCCGATAAAAAAAAAAAAGCGGTTCAATTGATGTGGTTTTTATCACTTTTAttt                                                          >  2:211976/1‑139 (MQ=255)
                                                                                cGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTCAAAATAAAAAAAAAGATCCGATAAAAATAAAACAACGTTTCAATTTATGTGGTTTTTATCACTTTTAttt                                                          >  2:135924/1‑139 (MQ=255)
                                                                                 gTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGa                                                         <  1:451952/140‑1 (MQ=255)
                                                                                 gTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGa                                                         <  1:135924/140‑1 (MQ=255)
                                                                                     aGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTaa                                                     <  2:709157/140‑1 (MQ=255)
                                                                                       gCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAaat                                                   <  2:834513/140‑1 (MQ=255)
                                                                                        ccGCATATCAGCTTAAAAAATGAACCATCGCCAAAGGCGGTGGGTTTTTTGTTGTCAATTTTAAAAAAAAAACAAAGATCCGGTAAAAATAAAAAAACGTTTCAATTTGTGTGAGTTTTATCACTTTTTTTGAtgtaata                                                  >  2:1357061/1‑134 (MQ=255)
                                                                                         cGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTTTGATCAATTTTAAAATAAAAAAAAAGATCCGATAAAAATAAAACAACGTTTCAATTGATGTGGTTTTGATCACTTTTTTTGATTAaatga                                                 >  2:1504392/1‑140 (MQ=255)
                                                                                         cGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAaatga                                                 <  2:1255573/140‑1 (MQ=255)
                                                                                         cGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTTTAAAAAAAAACAAAGATCCGATAAAAAAAAAAAAGCGTTTCAATTGATGTGGGTTTTATCACTTTTTTTTATTaaaaaa                                                 >  2:451952/1‑137 (MQ=255)
                                                                                           cATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGGGGGTTTTTTTTGATCAAATTTAAAAAAAAAACAAAGATCCGATAAAAAAAAAAAAACGGTTCAATTGATGTGGTTTTGATTACTTTTATTTATTAaatgaaa                                               >  2:161153/1‑139 (MQ=255)
                                                                                            atatCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGGTCAATTTTAAAATAAAAAAAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTTATCACTTTTATTGATTAaatgaatg                                              >  1:1428668/1‑140 (MQ=255)
                                                                                             tatCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTGTGATCAATTTTAAAATAAAAAAAATGATCCGATAAAAATAAAAAAGCGTTTCAATTGATGTGGGTTTTATCACTTTTATTTATTAAATGAATGt                                             >  1:1466206/1‑140 (MQ=255)
                                                                                              atCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTc                                            <  1:161153/140‑1 (MQ=255)
                                                                                                        aaaaTGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGtt                                  <  1:1357061/140‑1 (MQ=255)
                                                                                                        aaaaTGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGtt                                  <  1:50961/140‑1 (MQ=255)
                                                                                                           aTGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTcc                               <  1:870872/140‑1 (MQ=255)
                                                                                                             gAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCAt                             <  2:238379/140‑1 (MQ=255)
                                                                                                             gAACCATCGCCAACGGCGGTGGGTTTTTTGTGATCAATTTTAAAATAAAAACAATGATCCCATAAAAAAAAAAAAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCAt                             >  1:798215/1‑140 (MQ=255)
                                                                                                                       cAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGAt                   <  1:130102/140‑1 (MQ=255)
                                                                                                                        aaCGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATg                  <  2:1428668/140‑1 (MQ=255)
                                                                                                                          cggcggTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCt                <  2:1466206/140‑1 (MQ=255)
                                                                                                                                  gtTTTTTTGTGATCAATTTCAAAATAAAAACAGTGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAg        <  2:745387/140‑1 (MQ=255)
                                                                                                                                   tttttttGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAgg       <  1:1436728/140‑1 (MQ=255)
                                                                                                                                        ttGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGGAATTa  <  1:203415/140‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                         
CTGGAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGGAATTA  >  NC_000913/4101376‑4101651

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