Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,457,094 A→G I240V (ATC→GTC)  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,0940AG100.0% 137.7 / NA 37I240V (ATC→GTC) gspDgeneral secretory pathway component, cryptic
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/22);  total (15/22)

CTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTATCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTGAAATATGCTAAAGCCACGAATCTGGTGGAA  >  NC_000913/3456965‑3457225
                                                                                                                                 |                                                                                                                                   
cTGGAATATGCGTCAGCAGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTg                                                                                                                            >  2:449222/1‑139 (MQ=255)
  ggAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGa                                                                                                                          <  2:180327/139‑1 (MQ=255)
         gCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGaaa                                                                                                                   <  1:449222/139‑1 (MQ=255)
              aGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCa                                                                                                              >  2:506812/1‑139 (MQ=255)
              aGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCa                                                                                                              >  1:7232/1‑139 (MQ=255)
                 ggAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGc                                                                                                           <  1:506812/139‑1 (MQ=255)
                  gAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGcc                                                                                                          <  2:273780/139‑1 (MQ=255)
                              gagaTTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGTCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACtt                                                                                              <  2:513702/139‑1 (MQ=255)
                               agaTTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACtt                                                                                              <  2:389302/138‑1 (MQ=255)
                               agaTTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACtt                                                                                              >  1:389302/1‑138 (MQ=255)
                                 aTTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGc                                                                                                           >  2:557968/1‑123 (MQ=255)
                                 aTTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGc                                                                                                           <  1:557968/123‑1 (MQ=255)
                                             ttAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGcc                                                                               >  1:712798/1‑139 (MQ=255)
                                                aTCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTg                                                                            <  1:113515/139‑1 (MQ=255)
                                                       aaaGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGa                                                                     <  2:7232/139‑1 (MQ=255)
                                                        aaGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGaa                                                                    <  2:571663/139‑1 (MQ=255)
                                                               ggTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGa                                                                                                                          >  2:223947/1‑78 (MQ=255)
                                                               ggTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGa                                                                                                                          <  1:223947/78‑1 (MQ=255)
                                                               ggTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCg                                                             <  1:187395/139‑1 (MQ=255)
                                                                  aaaaGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGagg                                                          <  2:715361/139‑1 (MQ=255)
                                                                           aTGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCGGCGCAt                                                                                                   <  1:67408/89‑1 (MQ=255)
                                                                           aTGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCAt                                                                                                   >  2:67408/1‑89 (MQ=255)
                                                                           aTGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATa                                                 >  2:197872/1‑139 (MQ=255)
                                                                                         ctcCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTAc                                   <  2:309792/139‑1 (MQ=255)
                                                                                            cgcgAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAg                                                                                                                >  2:317613/1‑59 (MQ=255)
                                                                                            cgcgAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAg                                                                                                                <  1:317613/59‑1 (MQ=255)
                                                                                                   aTTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGttt                                        <  2:109374/124‑1 (MQ=255)
                                                                                                   aTTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGttt                                        >  1:109374/1‑124 (MQ=255)
                                                                                                     tgtgGCGGCTAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTGAAATATGCt                       <  1:91230/139‑1 (MQ=255)
                                                                                                     tgtgGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTGAAATATGCt                       >  2:290277/1‑139 (MQ=255)
                                                                                                              tAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTGAAATATGCTAAAGCCACg              <  2:712798/139‑1 (MQ=255)
                                                                                                                       cAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGaa                                                                    <  2:105943/76‑1 (MQ=255)
                                                                                                                       cAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGaa                                                                    >  1:105943/1‑76 (MQ=255)
                                                                                                                       cAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTGAAATATGCTAAAGCCACGAATCtggtg     >  2:437683/1‑139 (MQ=255)
                                                                                                                         aCTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTAcc                                  <  1:341921/108‑1 (MQ=255)
                                                                                                                         aCTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTAcc                                  >  2:341921/1‑108 (MQ=255)
                                                                                                                          ctctctTGTCATCAGTGGACCGGACAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTGAAATATGCTAAAGCCACGAATCTGGTGGaa  <  1:437683/139‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                   
CTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTATCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTGAAATATGCTAAAGCCACGAATCTGGTGGAA  >  NC_000913/3456965‑3457225

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