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

GTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGGAA  >  NC_000913/4101386‑4101648
                                                                                                                                |                                                                                                                                      
gtgAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAAtttt                                                                                                                             >  2:1162570/1‑139 (MQ=255)
  gAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCaa                                                                                                                           <  1:347522/140‑1 (MQ=255)
       gggACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATaa                                                                                                                      >  1:131463/1‑140 (MQ=255)
        ggACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATaaa                                                                                                                     <  1:1162570/140‑1 (MQ=255)
        ggACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATaaa                                                                                                                     >  1:1167172/1‑140 (MQ=255)
             aaGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGCGGTTTTTTTGTGATCAATTTCAAAATAAAAACaa                                                                                                                >  1:928018/1‑140 (MQ=255)
              agcagcGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAAt                                                                                                               >  1:96675/1‑140 (MQ=255)
                 agcGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGAt                                                                                                            <  2:131463/140‑1 (MQ=255)
                   cGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTTAAAATaaaaaaaaaaatac                                                                                                          >  2:1214255/1‑134 (MQ=255)
                      cACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGat                                                                                                       <  2:844237/140‑1 (MQ=255)
                      cACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCAACGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGat                                                                                                       <  2:183202/140‑1 (MQ=255)
                        cGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTTTGATCAATTTCAAAATAAAAACAATGATCCGataa                                                                                                     >  1:27728/1‑140 (MQ=255)
                          ttttACGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTTAAAATAAAAAAAAAGATCCGataaaa                                                                                                   >  2:2670/1‑140 (MQ=255)
                                      aaaaTAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGt                                                                                       <  2:1167172/140‑1 (MQ=255)
                                       aaaTAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGtt                                                                                      <  1:779123/140‑1 (MQ=255)
                                        aataatCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGttt                                                                                     <  2:781143/140‑1 (MQ=255)
                                           aatCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCaa                                                                                  <  2:27728/140‑1 (MQ=255)
                                            atCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAAt                                                                                 <  1:1171375/140‑1 (MQ=255)
                                            atCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAAt                                                                                 <  2:819949/140‑1 (MQ=255)
                                                tttGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGAt                                                                             <  2:236668/140‑1 (MQ=255)
                                                 ttGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGAtg                                                                            <  2:1214297/140‑1 (MQ=255)
                                                 ttGCCGCCTGCTGCAATGAGGCGTATAGGCAGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGGTCAATTTTAAAATAAAAAAAAAGATCCGATAAAAAAAAAAAAGCGTTTCAATTTAtg                                                                            >  2:633803/1‑140 (MQ=255)
                                                    ccgccTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTgg                                                                         <  1:240139/140‑1 (MQ=255)
                                                           ctgcAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGAt                                                                  <  2:1234891/140‑1 (MQ=255)
                                                           ctgcAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGAt                                                                  <  1:373125/140‑1 (MQ=255)
                                                           ctgcAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTTAAAATAAAAAAAAAGATCCCATAAAAAAAAAAAAACGGTTCAATTGATGTGGGTTTTAt                                                                  >  1:552169/1‑140 (MQ=255)
                                                                      cGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTTTGGTCAATTTTAAAATAAAAAAAATGATCCGATAAAAATAAAAAAGCGTTTCAATTGATGTGGTTTTTATCACtttgtttt                                                       >  2:1159435/1‑135 (MQ=255)
                                                                      cGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGGTCAATTTTAAAATAAAAAAAATGATCCGATAAAAATAAAACAGCGTTTTAATTGATGTGGTTTTTATCACTTTTAttt                                                       >  2:1061202/1‑139 (MQ=255)
                                                                        tataGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTTTGATCAATTTTAAAATAAAAAAAAAGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGAt                                                     >  2:752091/1‑140 (MQ=255)
                                                                        tataGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGGTCCATTTCTAAATAAAAAAAATGATTCCATAAAAAAAAAACACCGTTTTCATTTATGTGGGTTTTATCACTTTTATTGAt                                                     >  2:736825/1‑140 (MQ=255)
                                                                         ataGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGAtt                                                    <  2:898973/140‑1 (MQ=255)
                                                                             gCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAaat                                                <  2:726716/140‑1 (MQ=255)
                                                                             gCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAaat                                                <  1:716710/140‑1 (MQ=255)
                                                                             gCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGGTCAATTTTAAAAAAAAAAAAAATATCCGATAAAAAAAAAAAAACGTTTTAATTTATGTGGGTTTTATCACTTTTTTTTATTaaaa                                                >  2:601122/1‑139 (MQ=255)
                                                                               cGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAaatga                                              <  2:774843/140‑1 (MQ=255)
                                                                               cGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTTAAAATAAAAAAAATGATCCGATAAAAATAAAACAGCGTTTTAATTGATGTGGTTTTGATCACTTTTTTTGATTAaatga                                              >  1:975293/1‑140 (MQ=255)
                                                                                gCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAaatgaa                                             <  2:96675/140‑1 (MQ=255)
                                                                                gCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAaatgaa                                             <  1:1061202/140‑1 (MQ=255)
                                                                                 cATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAaat                                                >  1:77686/1‑136 (MQ=255)
                                                                                 cATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAaat                                                <  2:77686/136‑1 (MQ=255)
                                                                                  atatCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTTAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTTATTAaatgaatg                                           >  2:1014662/1‑140 (MQ=255)
                                                                                       aGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTTTGATCAATTTTAAAAAAAAAACAAAGATCCGATAAAAAAAAAAAAGCGGTTCAATTGATGTGGTTTTGATTACTTTTTTTGATTAAATAAATGTCttt                                      >  1:1209335/1‑138 (MQ=255)
                                                                                       aGCTTAAAAAATGAACCATCGCCAACGGCGGCGGTTTTTTTTGGATCAAATTCAAAATAAAAACAATGATCCGATAAAAAAAAAACAGCGTTTCAATTGATGTGGTTTTGATCACATTTATTTTTTAAATGAATGTCTAt                                      >  2:716710/1‑140 (MQ=255)
                                                                                         cTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCt                                    <  1:146313/140‑1 (MQ=255)
                                                                                           tAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAAATTTAAAAAAAAAAAAAAGATCCGATAAAAATAAAACAGCGTTTTAATTTATGTGGGTTTTATCACTTTTATTTATTAAATTAATGTCTATCttt                                  >  2:146313/1‑139 (MQ=255)
                                                                                              aaaaTGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGtt                               <  1:1214255/140‑1 (MQ=255)
                                                                                                   gAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCAt                          <  2:1209335/140‑1 (MQ=255)
                                                                                                   gAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCAt                          <  2:552169/140‑1 (MQ=255)
                                                                                                     aCCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCa                        <  2:728753/140‑1 (MQ=255)
                                                                                                      ccATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCaa                       <  1:736825/140‑1 (MQ=255)
                                                                                                         tCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAacac                    <  1:758537/140‑1 (MQ=255)
                                                                                                          cGCCAACGGCGGTGGGTTTTTTTTGATCAATTTCAAAATAAAAAAAAAGATCCGATAAAAATAAAAAAGCGTTTTAATTTATGTGGTTTTTATCACTTTTATTTATTAAATGAATGTCTATCTTTGTTTTCATCAACACt                   >  2:758537/1‑140 (MQ=255)
                                                                                                                cggcggTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCt             <  2:155856/140‑1 (MQ=255)
                                                                                                                cggcggTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCt             <  1:1159435/140‑1 (MQ=255)
                                                                                                                 ggcggTGGGTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAAAAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTc            >  1:123222/1‑140 (MQ=255)
                                                                                                                   cggTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCa          <  1:2670/140‑1 (MQ=255)
                                                                                                                           tttttGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGGaa  >  2:573007/1‑140 (MQ=255)
                                                                                                                               tgtgATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTg       <  2:1152800/131‑1 (MQ=255)
                                                                                                                               tgtgATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTg       >  1:1152800/1‑131 (MQ=255)
                                                                                                                                |                                                                                                                                      
GTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGGAA  >  NC_000913/4101386‑4101648

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