Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,457,616 G→A A414T (GCC→ACC)  gspD → general secretory pathway component, cryptic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,457,6160GA100.0% 126.6 / NA 40A414T (GCC→ACC) gspDgeneral secretory pathway component, cryptic
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/19);  total (21/19)

GTTGGCGCACAGCAATTTACCAATACCGGATTACCGATTTTTAACGCTGCGCAAGGTGTGGCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGATAATAAACTCGCGTCCTTCA  >  NC_000913/3457496‑3457751
                                                                                                                        |                                                                                                                                       
gTTGGCGCACAGCAATTTACCAATACCGGATTACCGATTTTTAACGCTGCGGAAGGTGTGGCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCCACAATGGCATGGCCg                                                                                                                       <  2:718372/139‑1 (MQ=255)
            cAATTTACCAATACCGGATTACCGATTTTCAACGCTGAGCAAGATGTGGCTGATTATAAAAAGAATGGTGGGATGACCAGCGCGAATCCTGCCTGGGAGATGTTTAGCACCTACAATGGCATGGCCGCAGGCttcttcc                                                                                                           >  1:718372/1‑138 (MQ=255)
              aTTTACCAATACCGGATTACCGATTTTTAACGCTGCGCAAGGTGTGGCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAAt                                                                                                         <  1:602783/139‑1 (MQ=255)
                    cAATACCGGATTACCGATTTTTAACGCTGCGCAAGGTGTGGCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGAc                                                                                                   >  2:260850/1‑139 (MQ=255)
                      aTACCGGATTACCGATTTTTAACGCTGCGCAAGGTGTGGCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCa                                                                                                                     >  2:406858/1‑119 (MQ=255)
                      aTACCGGATTACCGATTTTTAACGCTGCGCAAGGTGTGGCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCa                                                                                                                     <  1:406858/119‑1 (MQ=255)
                                          aaCGCTGCGCAAGGTGTGGCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCgcgc                                                                             <  1:639717/139‑1 (MQ=255)
                                               tGCGCAAGGTGTGGCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCg                                                                                                                       <  2:1067498/92‑1 (MQ=255)
                                               tGCGCAAGGTGTGGCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCg                                                                                                                       >  1:1067498/1‑92 (MQ=255)
                                               tGCGCAAGGTGTGGCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGc                                                                                                                  >  1:703660/1‑97 (MQ=255)
                                               tGCGCAAGGTGTGGCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGcc                                                                        >  2:31676/1‑139 (MQ=255)
                                               tGCGCAAGGTGTGGCTGATTATAAAAAGAATGGTAGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGACGCAGGc                                                                                                                  <  2:703660/97‑1 (MQ=255)
                                                      ggTGTGGCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACa                                                                 <  1:31676/139‑1 (MQ=255)
                                                                aTTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTggg                                                                                               >  2:652081/1‑99 (MQ=255)
                                                                aTTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTggg                                                                                               <  1:652081/99‑1 (MQ=255)
                                                                  tataAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGc                                                                                                                  >  1:18/1‑78 (MQ=255)
                                                                  tataAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGc                                                                                                                  <  2:18/78‑1 (MQ=255)
                                                                  tataAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACa                                                     >  1:421137/1‑139 (MQ=255)
                                                                           aaTGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTa                                                                                   <  1:390743/100‑1 (MQ=255)
                                                                           aaTGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTa                                                                                   >  2:390743/1‑100 (MQ=255)
                                                                            aTGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTC‑‑CAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGc                                              <  2:85212/134‑1 (MQ=255)
                                                                            aTGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTC‑‑CAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGc                                              >  1:85212/1‑134 (MQ=255)
                                                                                   gATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGc                                    >  2:174998/1‑139 (MQ=255)
                                                                                     tCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCAt                                  >  1:926593/1‑139 (MQ=255)
                                                                                     tCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGACATGGCCg                                                                                                                       <  2:611760/54‑1 (MQ=255)
                                                                                     tCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGACATGGCCg                                                                                                                       >  1:611760/1‑54 (MQ=255)
                                                                                       aCCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGa                                                                                             <  2:442475/78‑1 (MQ=255)
                                                                                       aCCAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGa                                                                                             >  1:442475/1‑78 (MQ=255)
                                                                                         cAGCGCGAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTa                              >  2:162446/1‑139 (MQ=255)
                                                                                             gcgAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTggg                                                                                               <  1:724808/70‑1 (MQ=255)
                                                                                             gcgAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTggg                                                                                               >  2:724808/1‑70 (MQ=255)
                                                                                              cgAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGa                                                                                                    >  2:226559/1‑64 (MQ=255)
                                                                                              cgAATCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGa                                                                                                    <  1:226559/64‑1 (MQ=255)
                                                                                               caaTCCTTACTGGGAGATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTGACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTg                        <  2:63067/138‑1 (MQ=255)
                                                                                                 aTCCTGCCTGGGATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGa                      >  2:46684/1‑139 (MQ=255)
                                                                                                           ggATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGa                                                                                                    <  2:274463/51‑1 (MQ=255)
                                                                                                           ggATATGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGa                                                                                                    >  1:274463/1‑51 (MQ=255)
                                                                                                              tatGTTTAGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGATAATAAACTcgcg         <  1:174998/139‑1 (MQ=255)
                                                                                                                     aGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACg                           <  2:519726/114‑1 (MQ=255)
                                                                                                                     aGCACCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGATAATAAACTCGCGTCCTTCa  >  1:578933/1‑139 (MQ=255)
                                                                                                                        |                                                                                                                                       
GTTGGCGCACAGCAATTTACCAATACCGGATTACCGATTTTTAACGCTGCGCAAGGTGTGGCTGATTATAAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGATAATAAACTCGCGTCCTTCA  >  NC_000913/3457496‑3457751

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