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

TGGAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGGAA  >  NC_000913/4101377‑4101648
                                                                                                                                         |                                                                                                                                       
tGGAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTtgtg                                                                                                                                       <  1:378373/140‑1 (MQ=255)
       tggtgAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAAtt                                                                                                                                >  1:15477/1‑140 (MQ=255)
       tggtgAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAAtt                                                                                                                                >  1:791495/1‑140 (MQ=255)
         aggAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTc                                                                                                                              <  1:49821/138‑1 (MQ=255)
                gggACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATaa                                                                                                                      <  1:727933/140‑1 (MQ=255)
                   aCGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTTTGATCAATTTCaaa                                                                                                                          >  1:1143685/1‑133 (MQ=255)
                   aCGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCaaa                                                                                                                          <  2:1143685/133‑1 (MQ=255)
                      aaGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGTCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTTAAAATAAAAACaa                                                                                                                >  1:733289/1‑140 (MQ=255)
                      aaGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTTAAAATAAAAACaa                                                                                                                >  2:909060/1‑140 (MQ=255)
                      aaGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAAATTTaaaaaaaaaaact                                                                                                                >  2:49821/1‑137 (MQ=255)
                       agcagcGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAAt                                                                                                               <  2:17052/140‑1 (MQ=255)
                       agcagcGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAAt                                                                                                               <  1:659125/140‑1 (MQ=255)
                           gcGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTTAAAATAAAAAAAATGATc                                                                                                           >  1:260435/1‑140 (MQ=255)
                            cGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATaaa                                                                                                                     <  2:108296/129‑1 (MQ=255)
                            cGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATaaa                                                                                                                     >  1:108296/1‑129 (MQ=255)
                               cACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGat                                                                                                       <  2:563190/140‑1 (MQ=255)
                                  gTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTTAAAATAAAAAAAAAGATCCGataaa                                                                                                    >  2:31795/1‑140 (MQ=255)
                                            gAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAg                                                                                          <  1:909060/140‑1 (MQ=255)
                                               aaaaTAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGt                                                                                       <  2:15477/140‑1 (MQ=255)
                                                   taatCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCa                                                                                   <  1:1563829/140‑1 (MQ=255)
                                                     atCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAAt                                                                                 <  1:1177515/140‑1 (MQ=255)
                                                          ttGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGAtg                                                                            <  1:31795/140‑1 (MQ=255)
                                                          ttGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGAtg                                                                            <  2:518868/140‑1 (MQ=255)
                                                          ttGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGAtg                                                                            <  2:171544/140‑1 (MQ=255)
                                                          ttGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGAtg                                                                            <  2:746013/140‑1 (MQ=255)
                                                          ttGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGAtg                                                                            <  1:1492140/140‑1 (MQ=255)
                                                           tGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGGTCAATTTTAAAAAAAAAAAAAAGATCCGATAAAAATAAAAAAGCGGTTCAATTGAtgg                                                                           >  2:1563829/1‑139 (MQ=255)
                                                            gccgccTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTTAAAATAAAAAAAATGATCCGATAAAAATAAAAAAACGTTTCAATTGAtgtg                                                                          >  1:1142273/1‑140 (MQ=255)
                                                             ccgccTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTTAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTgg                                                                         >  1:561098/1‑140 (MQ=255)
                                                                    ctgcAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTATTTTTGATCAATTTTAAAATAAAAAAAAAGATCCGGTAAAAATAAAAAAGCGTTTCAATTTATGTGGGTTTGAt                                                                  >  2:1429249/1‑140 (MQ=255)
                                                                      gagacagggCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCa                                                                <  1:162362/133‑1 (MQ=255)
                                                                       cAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCAc                                                               <  2:912475/140‑1 (MQ=255)
                                                                          tGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTTTGATCAATTTTAAAATAAAAAAAAAAATCCCATAAAAATAAAAAAACGTTTCAATTTATTTGGGTTTTATCACttt                                                            >  2:943232/1‑140 (MQ=255)
                                                                          tGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTTAAAAAAAAAAAAAAAATTCGATAAAAAAAAAAAAACGGTTCAATTGATGTGGTTTTTATTACttt                                                            >  2:774348/1‑140 (MQ=255)
                                                                            gggCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTa                                                          <  1:1557221/139‑1 (MQ=255)
                                                                            aGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTa                                                          <  1:959925/140‑1 (MQ=255)
                                                                            aGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTa                                                          <  1:945986/140‑1 (MQ=255)
                                                                            aGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTa                                                          <  2:1485083/140‑1 (MQ=255)
                                                                               cGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTGTGATCAATTTTAAAATAAAAACAAAGATCCGATAAAAATAAAACAGCGGTTCAATTGATGTGGTTTTTATCACTTTTAttt                                                       >  1:644804/1‑139 (MQ=255)
                                                                                 tataGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAAATTTAAAATAAAAAAAATGATCCGGTAAAAATAAAAAAGCGTTTTAATTTATGTGGTTTTTATTACTTTTTTTGAt                                                     >  2:35217/1‑140 (MQ=255)
                                                                                  ataGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGAtt                                                    <  1:832459/140‑1 (MQ=255)
                                                                                      gCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTTAAAATAAAAACAATGATCCGATAAAAATAAAAAAGCGTTTCAATTGATGTGGTTTTGATCACTTTTTTTGATTAaat                                                >  1:91709/1‑140 (MQ=255)
                                                                                      gCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGGGGGTTTTTTTTGATCAATTTCAAAATAAAAACAAAGATCCGGTAAAAAAAAAAAAGCGTTTTAATTTATGGGGTTTTGATCACTTTTATTTATTaaaa                                                >  2:959328/1‑139 (MQ=255)
                                                                                        cGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTTAAAAAAAAAAAAAAGATCCGATAAAAAAAAAAAAGCGTTTTAATTTATGTGGGTTTGATTACTTTTATTTATTAaatgc                                              >  2:1557221/1‑139 (MQ=255)
                                                                                        cGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTCAAAAAAAAAACAAAGATCCGATAAAAAAAAAACAACGTTTCAATTGATGTGGTTTTGATTACTTttatttatctgtctc                                              >  2:162362/1‑133 (MQ=255)
                                                                                          cATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTTAAAATAAAAAAAAAGATCCGATAAAAATAAAAAAACGTTTTAATTGATGTGGTTTTTATCACTTTTTTTTATTAaatgaat                                            >  2:83122/1‑140 (MQ=255)
                                                                                          cATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAaatgaat                                            <  2:597399/140‑1 (MQ=255)
                                                                                           atatCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTCAAAATAAAAAAAATGATCCGATAAAAATAAAAAAGCGTTTCAATTGATGTGGTTTTGATCACTTTTTTTTATTAaatgaatg                                           >  1:182368/1‑140 (MQ=255)
                                                                                                  cTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATTAATTTTAAAATAAAAAAAAAGATCCGAAAAAAAAAAAACAGCGTTTTAATTTATGTGGGTTTTATTACTTTTTTTTATTAAATTAATGTCTATCt                                    >  2:1124223/1‑140 (MQ=255)
                                                                                                       aaaaTGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGtt                               <  1:35217/140‑1 (MQ=255)
                                                                                                       aaaaTGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGtt                               <  1:943232/140‑1 (MQ=255)
                                                                                                       aaaaTGAACCATCGCCAACGGCGGTGGGTTTTTTGTGATCAATTTTAAAAAAAAAAAAAAGATCCGATAAAAAAAAAAAAACGGTTTAATTTATGTGGGTTTGATCACTTTTATTGATTAAAAGAAAGTCTATCTTCGtt                               >  2:741801/1‑140 (MQ=255)
                                                                                                            gAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTTAAAATAAAAAAAATGATCCGATAAAAATAAAAAAGCGTTTCAATTTATGTGGGTTTGATTACTTTTATTGATTAAATGAATGTCTATCTTTGTTTCCAt                          >  2:881526/1‑140 (MQ=255)
                                                                                                            gAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCAt                          <  2:1326605/140‑1 (MQ=255)
                                                                                                                cATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCaa                       <  2:672771/140‑1 (MQ=255)
                                                                                                                  tCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAacac                    <  2:91709/140‑1 (MQ=255)
                                                                                                                            cggTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCa          <  1:1429249/140‑1 (MQ=255)
                                                                                                                               tggtTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTg       <  1:959328/140‑1 (MQ=255)
                                                                                                                                 gtTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAg     <  2:1142273/140‑1 (MQ=255)
                                                                                                                                 gtTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAg     <  1:83122/140‑1 (MQ=255)
                                                                                                                                    ttttttTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGGaa  >  2:773235/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                       
TGGAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGGAA  >  NC_000913/4101377‑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: