Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,911,819 T→G 9.1% intergenic (‑294/+20) pstS ← / ← glmS phosphate ABC transporter periplasmic binding protein/L‑glutamine‑‑D‑fructose‑6‑phosphate aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,911,8190TG 9.1% 114.2 / 2.3 43intergenic (‑294/+20)pstS/glmSphosphate ABC transporter periplasmic binding protein/L‑glutamine‑‑D‑fructose‑6‑phosphate aminotransferase
Reads supporting (aligned to +/- strand):  ref base T (14/24);  new base G (2/3);  total (16/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 = 9.96e-01

TCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCAC  >  NC_000913/3911677‑3911958
                                                                                                                                              |                                                                                                                                           
tCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCt                                                                                                                                             >  2:73494/1‑143 (MQ=255)
       cAAACATTAATAACGAAGAGATGACGAACAAATTTACATTCTGTGGTAGAGACATAGTAGGTGAACTCGGTGGATTGTCACATGGAGTTGGCAGGATGTTTGACAGAAATTGTGTCAGGAAGAAAAATGCCCCGCTTACGCAg                                                                                                                                      <  2:449426/143‑1 (MQ=255)
        aaaCATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAgg                                                                                                                                     >  2:652973/1‑143 (MQ=255)
        aaaCATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAgg                                                                                                                                     <  1:265960/143‑1 (MQ=255)
                 taaCGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCTGCGTAAGCGGGGCAtttttc                                                                                                                            <  2:855013/143‑6 (MQ=255)
                  aaCGAAGAGACGACAGAAAAATTTTCATGCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCAttt                                                                                                                           <  2:501789/143‑1 (MQ=255)
                    ggAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATttat                                                                                                                         <  1:677518/142‑1 (MQ=255)
                      aaGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATttatta                                                                                                                       <  2:464602/143‑1 (MQ=255)
                       agggaggACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTGTGATTAACAACATAACAGGAAGAAAAATGCCCCGCTTAGGCAGGGCATCCATTTATTAc                                                                                                                      <  2:479310/137‑1 (MQ=255)
                        gagaTGACAGAAAAATTTTCAGTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTTACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACt                                                                                                                     <  2:561961/143‑1 (MQ=255)
                           aagACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATAACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCaa                                                                                                                  <  1:500899/141‑1 (MQ=255)
                                  aaaaaTTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAGACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTCTTACTCAACCGTAAc                                                                                                           >  1:669628/1‑143 (MQ=255)
                                    aaaTTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCg                                                                                                         >  1:36627/1‑143 (MQ=255)
                                       gtgtCATTCTCTGACAGGTGAAAAGTAGCCGAAGATGACGGTTTGTCACAGGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGAtt                                                                                                      <  2:910467/140‑1 (MQ=255)
                                              tCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAg                                                                                               <  1:226644/143‑1 (MQ=255)
                                                tgtgACAGAGAAGAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGt                                                                                             <  2:299937/143‑1 (MQ=255)
                                                     cAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTAcgc                                                                                        <  1:440080/143‑1 (MQ=255)
                                                      agagAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACAGATTTTACCAGATTAcaag                                                                                       >  1:806976/1‑140 (MQ=255)
                                                         gaAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCt                                                                                    <  2:769851/143‑1 (MQ=255)
                                                         gaAAAAGAAGCCGAAGATTACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCt                                                                                    <  2:528150/143‑1 (MQ=255)
                                                          aaaaaGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGCTCTTGCCAGGTTACGCGGCTg                                                                                   >  2:421075/1‑143 (MQ=255)
                                                               gTAGCCGACGATGGCGGTCTGTCACATGGAGTTGGCAGCGTGTTTGATTAAAAACATAACAGGAACAAATAGGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCaa                                                                              <  2:465443/143‑1 (MQ=255)
                                                                tAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCATGTCTTCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAAc                                                                             <  1:823748/143‑1 (MQ=255)
                                                                    cggcggtGACGGTTTGTCGCATGGAGTGGGGGGGGTGTGTGGTGAAGGGCGTAACGGGAAGAAAAAGGCCCGGCTTGCGCAGGGCATCCATTTATTACTCGACCGTAGCCGGTTTTGCCAGGTTACGCGGCTGGTCAACGTCg                                                                         <  2:244295/137‑1 (MQ=255)
                                                                    cGAAGGTGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCACCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCg                                                                         <  2:628430/143‑1 (MQ=255)
                                                                     gAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATCTATTACTCAACCGTAACCGATTTTGCCCGGTTACGCGGCTGGTCAACGTCgg                                                                        >  2:547801/1‑143 (MQ=255)
                                                                               gTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGAt                                                              >  1:656306/1‑143 (MQ=255)
                                                                                  tGTCACAGGAAGTTGCCAGGAGGTTTGATAAAAAACATACCAGGAGGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAg                                                           <  1:891057/143‑1 (MQ=255)
                                                                                        aTGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGAc                                                     >  1:5888/1‑143 (MQ=255)
                                                                                         tGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACa                                                    >  2:754319/1‑143 (MQ=255)
                                                                                                  cAGAATGTTTGAGTAAAAACATAACAGGAAGCAAAATGCCCCGCGTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGcc                                           <  1:228352/143‑1 (MQ=255)
                                                                                                   aGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCca                                          >  2:535027/1‑143 (MQ=255)
                                                                                                              ttAAAAACATAACAGGAAGAGACATGCCCCGCTTACGCAGGGCATCCATTTATGACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTATGCCAGCAGCTGCAGc                               <  2:774904/143‑1 (MQ=255)
                                                                                                                 aaaaCATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTCAGCCACCAGCTGCAGCgca                            >  1:788467/1‑141 (MQ=255)
                                                                                                                           aGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGgaga                  >  1:228356/1‑143 (MQ=255)
                                                                                                                            ggAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGGCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGCACGGTGTagaa                 >  1:128673/1‑143 (MQ=255)
                                                                                                                            ggAAGAAAAAGGCCCCGCGTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTagaa                 <  1:272795/143‑1 (MQ=255)
                                                                                                                                 aacaaTGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCg            <  2:416670/140‑1 (MQ=255)
                                                                                                                                 aaaaaTGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAGCGGTGTAGAAGATCg            <  2:5888/143‑1 (MQ=255)
                                                                                                                                   aaaTGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGt          <  1:73494/143‑1 (MQ=255)
                                                                                                                                    gaTGCCCTGCGTAAGCGGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAACAGCTGCAGCGGAACGGTGTAGAAGATAgggg         >  2:66243/2‑141 (MQ=255)
                                                                                                                                     aTGCCCTGCGTAAGCGGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCACCTGCAGCGGAACGGTGGAGAAGCTCgtcgc        >  2:768695/1‑139 (MQ=255)
                                                                                                                                     aTGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTGGGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGc        <  2:656306/143‑1 (MQ=255)
                                                                                                                                           cGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGAGCAATCAc  >  2:711805/1‑143 (MQ=255)
                                                                                                                                              |                                                                                                                                           
TCAATTACAAACATTAATAACGAAGAGATGACAGAAAAATTTTCATTCTGTGACAGAGAAAAAGTAGCCGAAGATGACGGTTTGTCACATGGAGTTGGCAGGATGTTTGATTAAAAACATAACAGGAAGAAAAATGCCCCGCTTACGCAGGGCATCCATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCAC  >  NC_000913/3911677‑3911958

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: