Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,908,514 A→G intergenic (‑148/+35) pstB ← / ← pstA phosphate ABC transporter ATPase/phosphate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,908,5140AG100.0% 151.0 / NA 44intergenic (‑148/+35)pstB/pstAphosphate ABC transporter ATPase/phosphate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (23/21);  total (23/21)

ATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGC  >  NC_000913/3908382‑3908638
                                                                                                                                    |                                                                                                                            
aTCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTa                                                                                                                        >  1:293474/1‑139 (MQ=255)
                   gCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCgtg                                                                                                     <  2:2282/139‑1 (MQ=255)
                            gCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCa                                                                                                  >  2:9350/1‑133 (MQ=255)
                            gCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCa                                                                                                  <  1:9350/133‑1 (MQ=255)
                            gCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAata                                                                                            >  2:499324/1‑139 (MQ=255)
                              ttACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATc                                                                                          >  2:102710/1‑139 (MQ=255)
                                       gCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGtt                                                                                 >  1:352340/1‑139 (MQ=255)
                                       gCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGtt                                                                                 >  1:146/1‑139 (MQ=255)
                                                 gCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGc                                                                       >  1:533020/1‑139 (MQ=255)
                                                          ccgcATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAaacaac                                                              <  2:208515/139‑1 (MQ=255)
                                                          ccgcATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAaacaac                                                              >  2:492104/1‑139 (MQ=255)
                                                           cgcATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGc                                                                                                <  2:523873/104‑1 (MQ=255)
                                                           cgcATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGc                                                                                                >  1:523873/1‑104 (MQ=255)
                                                           cgcATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcg                                                             <  1:243235/139‑1 (MQ=255)
                                                                    cAAAGTGTTAACCGTGTACAGAGTCAAAAGGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTGAGGGTGAGGGGAAACCCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGa                                                    <  2:392116/139‑1 (MQ=255)
                                                                    cAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGa                                                    <  2:146/139‑1 (MQ=255)
                                                                    cAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCATCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGa                                                    <  2:352340/139‑1 (MQ=255)
                                                                                  tgtACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagcagcag                                      >  2:651313/1‑139 (MQ=255)
                                                                                   gtACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGt                                                                                                                    <  1:646254/60‑1 (MQ=255)
                                                                                   gtACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGt                                                                                                                    >  2:646254/1‑60 (MQ=255)
                                                                                   gtACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGt                                     >  2:423788/1‑139 (MQ=255)
                                                                                        gagTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAAc                                                                                       <  1:503794/84‑1 (MQ=255)
                                                                                        gagTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAAc                                                                                       >  2:503794/1‑84 (MQ=255)
                                                                                        gagTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGca                                <  2:533020/139‑1 (MQ=255)
                                                                                         agTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGcac                               >  2:193089/1‑139 (MQ=255)
                                                                                         agTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGcac                               >  1:424053/1‑139 (MQ=255)
                                                                                         agTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGcac                               >  2:543020/1‑139 (MQ=255)
                                                                                         agTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGcac                               >  2:525213/1‑139 (MQ=255)
                                                                                          gTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCg                                                                        <  1:343487/97‑1 (MQ=255)
                                                                                          gTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCg                                                                        >  2:343487/1‑97 (MQ=255)
                                                                                          gTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGcaca                              <  2:293474/139‑1 (MQ=255)
                                                                                               aaaGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTa                         <  1:499324/139‑1 (MQ=255)
                                                                                               aaaGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTa                         <  2:424053/139‑1 (MQ=255)
                                                                                               aaaGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTa                         <  1:256688/139‑1 (MQ=255)
                                                                                                     gcTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATc                   >  1:606424/1‑139 (MQ=255)
                                                                                                     gcTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATc                   >  1:284057/1‑139 (MQ=255)
                                                                                                       tcggtcggTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCaa                 <  2:645479/139‑1 (MQ=255)
                                                                                                                   cTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCa                                                                                                  >  2:693106/1‑46 (MQ=255)
                                                                                                                   cTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCa                                                                                                  <  1:693106/46‑1 (MQ=255)
                                                                                                                      gCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGc  <  1:651313/139‑1 (MQ=255)
                                                                                                                         ccTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAatat                                                                                           <  2:507027/47‑1 (MQ=255)
                                                                                                                         ccTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAatat                                                                                           >  1:507027/1‑47 (MQ=255)
                                                                                                                         ccTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAat                                                                                             <  2:555195/45‑1 (MQ=255)
                                                                                                                         ccTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAat                                                                                             >  1:555195/1‑45 (MQ=255)
                                                                                                                                    |                                                                                                                            
ATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGC  >  NC_000913/3908382‑3908638

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: