Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,100,717 G→T 76.5% intergenic (‑192/‑93) rhaT ← / → sodA rhamnose/lyxose:H(+) symporter/superoxide dismutase (Mn)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,100,7170GT76.5% 85.8 / 26.4 49intergenic (‑192/‑93)rhaT/sodArhamnose/lyxose:H(+) symporter/superoxide dismutase (Mn)
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base T (16/20);  minor base A (5/6);  total (22/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.22e-01

TAAAAATCGTGCTGTCGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCATTTTCCTGCAAAACCATACCCTTACGAAAAGTACGGCATTGATAATCATTTTCAATATCATTTAATTAACTATAATGAACCAACTGCTTACGCGGCATTAACAATCGGCCGCCCGACAATACTGGAGATGAATATGAGCTATACCCTGCCATCCCTGCCGTATGCTTACGATGCCC  >  NC_000913/4100577‑4100852
                                                                                                                                            |                                                                                                                                       
tAAAAATCGTGCTGTCGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCATTTTCCTGCAACACCATCCCCTTCCGAAAAGTACGGCATTTAt                                                                                                                                       >  1:212313/1‑143 (MQ=255)
 aaaaaTCGTGCTGTCGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCATTTTCCTGCAAAGCCATACCCTTACGAAAAGTACGGCATTtata                                                                                                                                      >  2:559828/1‑143 (MQ=255)
 aaaaaTCGTGCTGTCGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCATTTTCCTGCAAAACCATACCCTTACGAAAAGTACGGCATTtata                                                                                                                                      <  2:112671/143‑1 (MQ=255)
         tGCTGTCGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCATTTTCCTGCAAAACCATACCCTTACGAAAAGTACGGCATTTATAATCAtttt                                                                                                                              <  2:625429/143‑1 (MQ=255)
           cTGTCGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCATTCTCCTGCAAAACCATACCCTTCCGAACAGTCCGCCATACATAATCATCTTCa                                                                                                                            >  2:837490/1‑143 (MQ=255)
           cTGGCGATGAACCTTTCTCCTGTTGCATCCTGTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGATTTGCTTCACATCTCTGGCATTTTCCTGCAAAACCATACCCTTACGAAAAGTACGGCATTTATAATCATTTTCa                                                                                                                            <  1:119920/143‑1 (MQ=255)
            tgacGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGAGTTACTGGCAATCACGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCATTTTCCTGCAAAACCATACCCTTACGAAAAGTACGGCATTTATAATCATTTTCaa                                                                                                                           <  2:597593/140‑1 (MQ=255)
            tGTCGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCATTTTCCTGCAAAACCATACCCTTACGAAAAGTACGGCATTTATAATCATTTTCaa                                                                                                                           <  1:262329/143‑1 (MQ=255)
                    aaCCTTTCGCCTCTTGCCGCCGTTGTCGATTTCCTGCCAATCACGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCATTTTCCTGCAAAACCATACCCTTACGAAAAGTACGGCATTTATAATCATTTCAATATCATtta                                                                                                                  <  1:447862/143‑1 (MQ=255)
                     aCCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCATTTTCCTGCAAAACCATACCCTTACGAAAAGTACGGCATTTATAATCATTTTCAATATCATtta                                                                                                                  >  1:100686/1‑143 (MQ=255)
                      ccTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGATTTGATTCACATCTCGGGCATTTTCCTGCAAAACCATATCCTTACGAAAAGTCGGGCATTTCTAATCATTTTTAACATGAAttaa                                                                                                                 >  2:577956/1‑143 (MQ=255)
                          ccGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCATTTTCCTGCAAAACCATACCCTTACGAAAAGTACGGCATTTATAATCATTTTCAATATCATttaattaa                                                                                                             <  2:85735/142‑1 (MQ=255)
                                gTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCATTTTCCTGCACAACCATACCCTTAGGACAGCTTCGGCATTTATAATCATTATCAATATTATTTAATTAAATATaa                                                                                                       >  1:510294/1‑143 (MQ=255)
                                gTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCATTTTCCTGCAAAACCATACCCTTACGAAAAGTACGGCATTTATAATCATTTTCAATATCATTTAATTAACTATaa                                                                                                       <  2:74603/143‑1 (MQ=255)
                                    ccgccgTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCATTTTCCTGCAAAACCATACCCTTACGAAAAGTACGGCATTAATAATCATTTTCAATATCATTTAATTAACTATAATGaa                                                                                                   >  1:680512/1‑143 (MQ=255)
                                    ccgccgTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCATTTTCCTGCAAAACAATACCCTTACGAAAAGTACGGCATTTATAATCATTTTCAATATCATTTAATTAACTATAATGaa                                                                                                   >  1:180374/1‑143 (MQ=255)
                                        cgTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCATTTTCCTGCAAAACCATACCCTTACGAAAAGTACGGCATTTATAATCATTTTCAATATCATTTAATTAACTATAATGAACCaa                                                                                               <  1:332998/143‑1 (MQ=255)
                                            gTCGATTTACTGGCAATCACGGCATTAAGTGGGTGAGTTGCTTCACATCTCGGGCATTTTCCTGCAAAACCATACCCTTACGAAAAGTACGGCATTTATAATCATTTTCAATATCATTTAATTAACTATAATGAACCAACTGc                                                                                           <  1:143767/143‑1 (MQ=255)
                                               gATTTACTGGCAATCACGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCATTTTCCTGCAAAACCATACCCTTACGAAAAGTATGGCATTTATAATCATTTTCAATATCATTTAATTAACTATAATGAACCACATGCTTa                                                                                        >  2:449531/1‑143 (MQ=255)
                                               gATTTACTGGCAATCACGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCATTTTCCTGCAAAACCATACCCTTACGAAAAGTACGGCATTAATAATCATTTTCAATATCATTTAATTAACTATAATGAACCAACTGCTTa                                                                                        <  2:862916/143‑1 (MQ=255)
                                                    aCTGGCAATCACGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCATTTTCCTGCAAAACCATACCCTTACGAAAAGTACGGCATTTATAATCATTTTCGATATCATTTAATTAACTATAATGAACCAACTGCTTACGCgg                                                                                   <  2:141260/143‑1 (MQ=255)
                                                     cTGGCAATCACGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCATTTTCCTGCAAAACCATACCCTTACGAAAAGTACGGCATTTATAATCATTTTCAATATCATTTAATTAACTATAATGAACCAACTGCTTACGCGGc                                                                                  <  1:32996/143‑1 (MQ=255)
                                                      tGGCAATCACGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCATTTTCCTGCAAAACCATACCCTTACGAAAAGTACGGCATTTATAATCATTTTCAATATCATTTAATTAACTATAATGAACCAACTGCTTACGCGGCa                                                                                 <  1:397787/143‑1 (MQ=255)
                                                        gCAATCACGGCATTAAGTGGGAGATTTGCTTCACATCTCGGGCATTTTCCTGCAAAACCATGCCCTTACGAAAAGTACGGCATTAATAATCATTTTCAATATCATTTACTTAACTATAATGAACAAGATCCTTACGCGGCAtt                                                                               >  2:211903/1‑143 (MQ=255)
                                                              aCGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCATTTTCCTGCAAAACCATACCCTTACGAAAAGTACGGCATTTATAATCATTTTCAATATCATTTAATTAACTATAATGAACCAACTGCTTACGCGGCATTAACAAt                                                                         >  1:96145/1‑143 (MQ=255)
                                                                   aTTAAGTGGGTGATTTACTTCACATCTCGCGCATTTTCCTGCAAAACCATACCCTTACGAAAAGTACGGCATTAATAATCATTTTCAATATCATTTAATTAACTATAATTAACCAACTGCTTACGCGGCATTAACAATCGgcc                                                                    <  1:322293/143‑1 (MQ=255)
                                                                      aaGTGGGTGGTTTGCTTCACATCTCGGGCATTTTCCTGCAAAACCATACCCTTACGAAAAGTACGGCATTAATAATCATTTTCAATATCATTTAATTAACTATAATGAACCAACTGCTTACGCGGCATTAACAATCGgccgcc                                                                 <  1:614569/143‑1 (MQ=255)
                                                                       agggggTGATTTGCTTCACATCTCGGGCATTTGCCTGCAAAACCATACCCTTACGAAAAGTACGGCATTAATAATCATTTTCAATATCATTTAATTAACTATAATGAACCAACTGCTTACGCGGCATTAACAATCGGCCGccc                                                                <  1:643228/140‑1 (MQ=255)
                                                                          gggTGATTTGCTTCACATCTCGGGCATTTTCCTGCAAAACCATACCCTTACGAAAAGTACGGCATTAATAATCATTTTCAATATCATTTAATTAACTATAATGAACCAACTGCTTACGCGGCATTAACAATCGGCCGCCCGAc                                                             >  1:457510/1‑143 (MQ=255)
                                                                                tttGCTTCACATCTCGGGCATTTTCCTGCAAAACCATACCCTTACGAAAAGTACGGCATTTATAATCATTTTCAATATCATTTAATTAACTATAATGAACCAACTGCTTACGCGGCATTAACAATCGGCCGCCCGACAATACt                                                       >  1:198327/1‑143 (MQ=255)
                                                                                tttGCTTCACATCTCGGGCATTTTCCTGCAAAACCATACACTTACGAAAAGTACGGCATTTGTAATCATTCTCCATTTTTTTTAATTAACTATAATGAACCAACTGCTTACGCGGCATTAACAATCGGCCGCCCGACAATACt                                                       >  1:878670/1‑143 (MQ=255)
                                                                                  tGCTTCACATCTCGGGCATTTTCCTGCAAAACCATACCCTTACGAAAAGTACGGCATTTATAATCATTTTCAATATCATTTAATTAACTATAATGAACCAACTGCTTACGCGGCATTAACAATCGGCCGCCCGACAATACTgg                                                     <  1:754837/143‑1 (MQ=255)
                                                                                   gCTTCACATCTCGGGCATTTTCCTGCAAAACCATACCCTTACGAAAAGTACGGCATTAATAATCATTTTCAATATCATTTAATTAACTATAATGAACCAACTGCTTACGCGGCATTAACAATCGGCCGCCCGACAATACTGga                                                    <  1:98860/143‑1 (MQ=255)
                                                                                   ccatCACATCTCGGGCATTTTCCTGCAAAACCATACCCTTACGAAAAGTACGGCATTTATAATCATTTTCAATATCATTTAATTAACTATAATGAACCAACTGCTTACGCGGCATTAACAATCGGCCGCCCGACAATACTGga                                                    <  1:866826/140‑1 (MQ=255)
                                                                                          aTCTCGGGCATTGTCCTGCAAAACCAGACCCTTACGAAAAGTACGGCATTGATAATCAGGTTCAATATCAGTTAATTAACTATAATGAACCAACTGCTTACGCGGCATTAACAATCGGCCGCCCGACAAGACTGGAGATGAat                                             <  1:16847/143‑1 (MQ=255)
                                                                                            cacGGGCATTGTCCTGCAAAACCATACCCTTACGAAAAGTACGGCATTTATAATCATTTTCAATATCATTTAATTAACTATAATGAACCAACTGCTTACGCGGCATTAACAATCGGCCGCCCGACAATACTGGAGATGAatat                                           <  1:281266/141‑1 (MQ=255)
                                                                                              cGGGCATTTTCCTGCAAAACCATACCCTTACGAAAAGTACGGCATTTATAATCATTTTCAATATCATTTAATTAACTATAATGAACCAACTGCTTACGCGGCATTAACAATCGGCCGCCCGACAATACTGGAGATGAATATGa                                         <  1:252097/143‑1 (MQ=255)
                                                                                                  cATTTTCCTGCAAAACCATACCCTTACGAAAAGTACGGCATTTATAATCATTTTCAATATCATTTAATTAACTATAATGAACCAACTGCTTACGCGGCATTAACAATCGGCCGCCCGACAATACTGGAGATGAATATGAGCta                                     >  2:354252/1‑143 (MQ=255)
                                                                                                      ttCCTGCAACACCATACCCTTACGAAAAGTACGGCATTAATAATCATTTTCAATATCATTTAATTAACTATAATGAACCAACTGCTTACGCGGCATTAACAATCGGCCGCCCGACAATACTGGAGATGAATATGAGCTATAcc                                 <  1:529637/143‑1 (MQ=255)
                                                                                                        ccTGCAAAACCATACCCTTACGAAAAGTACGGCATTTATAATCATTTTCAATATCATTTAATTAACTATAATGAACCAACTGCTTACGCGGCATTAACACTCGGCCGCCCGACAATACTGGAGATGAATATGAGCTATAcacg                               >  1:337964/1‑140 (MQ=255)
                                                                                                             aaaaCCATACCCTTACGAAAAGTACGGCATTTATAATCATTTTCAATATCATTTAATTAACTATAATGAACCAACTGCTTACGCGGCATTAACAATCGGCCGCCCGACAATACTGGAGATGAATATGAGCTATACCCTGCCAt                          <  1:859032/143‑1 (MQ=255)
                                                                                                              aaaCCATACCATTACGAAAAGTACGGCATTAATAATCATTTTCAATATCATTTAATTAACTAAAATGAACCAACTGCTTACGCAGCATTAACAATCGGCCGCCCGACAATACTGCAGACGAATATGAGATACACCCGCCCATc                         >  2:195323/1‑143 (MQ=255)
                                                                                                                  cATACCCTTACGAAAAGTACGGCATTTATAATCATTTTCAATATCATTTAATTAACTATAATGAACCAACTGCTTACGCGGCATTAACAATCGGCCGCCCGACAATACTGGAGATGAATATGAGCTATCCCGCGCCGTCCCTg                     >  1:116205/1‑143 (MQ=255)
                                                                                                                  cATACCCTTACGAAAAGTACGGCATTTATAATCATTTTCAATATCATTTAATTAACTATAATGAACCAACTGCTTACGCGGCATTAACAATCGGCCGCCCGACAATACTGGAGATGAATATGAGCTATACCCTGCCATCCCTg                     <  1:566314/143‑1 (MQ=255)
                                                                                                                        cTTACGAAAAGTACGGCATTTATAATCATTTTCAATATCATTTAATTAACTATAATGAACCAACTGCTTACGCGGCATTAACAATCGGCCGCCCGACAATACTGGAGATGAATATGAGCTATACCCTGCCATCCCTGCCGTAt               >  1:481038/1‑143 (MQ=255)
                                                                                                                        cTTACGAAAAGTACGGCATTAATAATCATTTTCAATATCATTTAATTCACTATAATGCACCAACTGCTTACACGACATTAACAATCGGCCGCCCGACAATACTGGAGATGAATATGAGCTATACCCTGCCATCCCTGCCGTAt               >  1:665737/1‑143 (MQ=255)
                                                                                                                            cGAAAAGTACGGCATTTATAATCATTTTCAATATCATTTAATTAACTATAATGAACCAACTGCTTACGCGGCATTAACAATCGGCCGCCCGACAATACTGGAGATGAATATGAGCTATACCCTGCCATCCCTGCCGTATGCtt           <  1:648279/143‑1 (MQ=255)
                                                                                                                             gAAAAGTACGGCATTTATAATCATTTTCAATATCATTTAATCAACTATAATGAACCAACTGCTTACGCGGCATTAACAATCGACCACCCGACAATACTGGAGATGAATATGAGCTATACCCTGCCATCCCTGCCGTATGCTTa          >  2:597880/1‑143 (MQ=255)
                                                                                                                              aaaaGTACGGCATTTATAAGCAGTTTCAATATGATTGAATTAACTATAATGAACCAACTGCTTACGCGGCATTAACAATCGGCCGCCCGACAATACTGGAGATGAATATGAGCTATACCCTGCCATCCCTGCCGTATGCTTAc         <  1:104160/143‑1 (MQ=255)
                                                                                                                                aaGTACGGCATTTATAATCATTTTCAATATAATTTAATTAACGATAAGTAACCAACTGCTTACGCGGCATTAACAATCGGCCGCCCGACAATACTGGAGATGAATATGAGCTATACCCTGCCATCCCTGCCGTATGCTTACGa       >  1:168381/1‑143 (MQ=255)
                                                                                                                                     cGGCATTTATAATCATTTTCAATATCATTTAATTAACTATAATGAACCAACTGCTTACGCGGCATTAACAATCGGCCGCCCGACAATACTGGAGATGAATATGAGCTATACCCTGCCATCCCTGCCGTATGCTTACGATGccc  >  1:573436/1‑143 (MQ=255)
                                                                                                                                            |                                                                                                                                       
TAAAAATCGTGCTGTCGATTAACCTTTCGCCTGTTGCCGCCGTTGTCGATTTACTGGCAATCACGGCATTAAGTGGGTGATTTGCTTCACATCTCGGGCATTTTCCTGCAAAACCATACCCTTACGAAAAGTACGGCATTGATAATCATTTTCAATATCATTTAATTAACTATAATGAACCAACTGCTTACGCGGCATTAACAATCGGCCGCCCGACAATACTGGAGATGAATATGAGCTATACCCTGCCATCCCTGCCGTATGCTTACGATGCCC  >  NC_000913/4100577‑4100852

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: