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

TCTGGAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGG  >  NC_000913/4101375‑4101646
                                                                                                                                           |                                                                                                                                    
tCTGGAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTtg                                                                                                                                      <  2:181845/140‑1 (MQ=255)
         tggtgAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTTTGATTAAtt                                                                                                                             >  2:142950/1‑140 (MQ=255)
         tggtgAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTTTGATCAAtt                                                                                                                             >  1:1257212/1‑140 (MQ=255)
              aaCTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCaaa                                                                                                                        <  2:695073/140‑1 (MQ=255)
              aaCTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTTaaa                                                                                                                        >  1:334288/1‑140 (MQ=255)
               aCTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTTaaaa                                                                                                                       <  1:1046350/140‑1 (MQ=255)
                   ggACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATaaa                                                                                                                   <  2:506384/140‑1 (MQ=255)
                   ggACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTGTGATCAATTTCAAAATaaa                                                                                                                   >  1:1392849/1‑140 (MQ=255)
                                 cACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGat                                                                                                     <  1:69514/140‑1 (MQ=255)
                                     ttttGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGataaaa                                                                                                 <  1:295802/140‑1 (MQ=255)
                                          cggcgAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGG‑GTTTTTGTGGTTAATTTTAAAATAAAAAAAATGAACCGataaaaataaaa                                                                                           >  2:172871/1‑140 (MQ=255)
                                            gcgAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTTTGATCAATTTTAAAATAAAAAAAATGATCCGataaaaataaaaa                                                                                          >  2:107810/1‑139 (MQ=255)
                                              gAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAg                                                                                        <  2:597719/140‑1 (MQ=255)
                                               aaaaaaTAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTTAAAATAAAAAAAATGATCCGataaaaaaaaaaaaac                                                                                       >  2:469001/1‑136 (MQ=255)
                                                       atCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAAt                                                                               <  1:914504/140‑1 (MQ=255)
                                                        tCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAAtt                                                                              <  1:461200/140‑1 (MQ=255)
                                                            ttGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGAtg                                                                          <  1:478425/140‑1 (MQ=255)
                                                            ttGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGAtg                                                                          <  2:606941/140‑1 (MQ=255)
                                                            ttGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGAtg                                                                          >  1:1347359/1‑140 (MQ=255)
                                                             tGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATTAAATTTAAAATAAAAACAAAGATCCGATAAAAATAAAAAAACGTTTTAATTTAtgt                                                                         >  2:22397/1‑140 (MQ=255)
                                                             tGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTTAAAATAAAAAAAAAGATCCGATAAAAATAAAAAAGCGTTTTAATTTAtgt                                                                         >  2:1035960/1‑140 (MQ=255)
                                                                cgccTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGt                                                                      <  2:1227174/140‑1 (MQ=255)
                                                                    tgctgcAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTg                                                                  <  2:4732/140‑1 (MQ=255)
                                                                      ctgcAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGAt                                                                <  1:1288599/140‑1 (MQ=255)
                                                                      ctgcAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGAt                                                                <  1:469001/140‑1 (MQ=255)
                                                                      ctgcAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGAt                                                                <  1:218205/140‑1 (MQ=255)
                                                                      ctgcAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGGTCAATTTCAAAAAAAAAAAAATGATCCGATAAAAAAAAAAAAACGGTTCAATTGATGTGGgttttgt                                                                >  2:235466/1‑137 (MQ=255)
                                                                        gcAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCa                                                              <  1:107810/140‑1 (MQ=255)
                                                                         cAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCAc                                                             <  1:837039/140‑1 (MQ=255)
                                                                         cAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCAc                                                             <  1:41475/140‑1 (MQ=255)
                                                                              aGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTTATCACttttt                                                        >  1:294284/1‑139 (MQ=255)
                                                                                gCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGGGGGTTTTTTTTGATCAATTTCAAAAAAAAAAAAAAGAACCGATAAAAAAAAAAAAGCGTTTCAATTTATGTGGGTTTGTTCACTtttttt                                                      >  2:1348953/1‑140 (MQ=255)
                                                                                  gTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGa                                                    <  2:1001915/140‑1 (MQ=255)
                                                                                   tataGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTTTGATCAATTTCAAAATAAAAAAAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGAt                                                   >  1:1338044/1‑140 (MQ=255)
                                                                                   tataGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTg                                                                  <  2:1238978/125‑1 (MQ=255)
                                                                                   tataGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTg                                                                  >  1:1238978/1‑125 (MQ=255)
                                                                                   tataGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGAt                                                   <  1:551558/140‑1 (MQ=255)
                                                                                          cGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAaatga                                            <  1:923503/140‑1 (MQ=255)
                                                                                          cGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATTAATTTTAAAAAAAAAACAATGATCCGATAAAAATAAAAAAGCGTTTTAATTTATGTGGTTTTTATCACTTTTTTTGATTAaatga                                            >  2:923503/1‑140 (MQ=255)
                                                                                          cGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTTAAAATAAAAAAAATGATCCGATAAAAAAAAAAAAACGGTTTAAATTATGTGGGTTTTATCAATTTTATTGATTAaaaga                                            >  2:307865/1‑140 (MQ=255)
                                                                                             atatCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTTTGATCAATTTCAAAAAAAAAAAAATGATCCGATAAAAATAAAACAGCGTTTTAATTGATGTGGTTTTTATCACTTTTTTTTATTAaatgaatg                                         >  2:565155/1‑140 (MQ=255)
                                                                                                       aaaaaaTGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTTAAAATAAAAAAAATGATCCGGTAAAAATAAAAAAGCGTTTTAATTTATGTGGGTTTTATTACTTTTTTTTATTTAATGAATGTCTATCTTct                               >  2:1304160/1‑139 (MQ=255)
                                                                                                         aaaaTGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGtt                             <  2:1392849/140‑1 (MQ=255)
                                                                                                                 ccATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCaa                     <  2:1013049/140‑1 (MQ=255)
                                                                                                                    tCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAacac                  <  2:20441/140‑1 (MQ=255)
                                                                                                                         aaCGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATg             <  1:142950/140‑1 (MQ=255)
                                                                                                                           cggcggTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCt           <  2:1458799/140‑1 (MQ=255)
                                                                                                                            ggcggTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTc          <  1:1240053/140‑1 (MQ=255)
                                                                                                                            ggcggTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTc          <  1:1035960/140‑1 (MQ=255)
                                                                                                                               ggtggtTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCAt       >  1:573650/1‑140 (MQ=255)
                                                                                                                                gtggtTTTTTTGTGATCAATTTCAAAATAAAAACAAAGATCCGATAAAAATAAAAAAGCGGTTTAATTGATGTGGGTTTTATCACTTTTTTTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCAtt      >  2:377243/1‑140 (MQ=255)
                                                                                                                                   gtTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAg   <  2:1338044/140‑1 (MQ=255)
                                                                                                                                    tttttttGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAgg  <  2:872373/140‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                    
TCTGGAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGG  >  NC_000913/4101375‑4101646

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: ‑