Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,099,504 G→A A21V (GCA→GTA)  yggM ← DUF1202 family putative secreted protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,099,5040GA100.0% 149.7 / NA 43A21V (GCA→GTA) yggMDUF1202 family putative secreted protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/24);  total (19/24)

GACCAGGTGGCCTGATTACCCTTAGCATCAAGATTTTTTAAGGTGATAGCATCAAGCTTGAGAATGCCGTGATACTGCTTTTTGAATTGATCTTTTAAGATATTTTCAGTTGGCGGTTCGCCATCTGCCCAGGCATTACCAGTCATCAACATAAGCAGTGCCGTGCCAACAATCCATTGCTTTTTCACCCGAAAAACCTCACGCGGAATAATTTATTCCGGTGATTATCATTAGGGTAATTAAAATTGGCAAT  >  NC_000913/3099378‑3099630
                                                                                                                              |                                                                                                                              
gACCAGGTGGCCTGATTACCCTTAGCATCAAGATTTTTTAAGGTGATAGCATCAAGCTTGAGAATGCCGTGATACTGCTTTTTGAATTGATCTTTTAAGATATTTTCAGTTGGCGGTTCGCCATCTACCCAGGCATTAc                                                                                                                    <  2:386548/139‑1 (MQ=255)
 aCCAGGTGGCCTGATTACCCTTAGCATCAAGATTTTTTAAGGTGGTAGCATCAAGCTTGAGAATGCCGTGATACTGCTTTTTGAATTGATCTTTTAAGATATTTTCAGTTGGCGGTTCGCCATCTACCCAGGCATTAcc                                                                                                                   <  1:632440/139‑1 (MQ=255)
      gTGGCCTGATTACCCTTAGCATCAAGATTTTTTAAGGTGATAGCATCAAGCTTGAGAATGCCGTGATACTGCTTTTTGAATTGATCTTTTAAGATATTTTCAGTTGGCGGTTCGCCATCTACCCAGGCATTACCAGtca                                                                                                              >  1:16015/1‑139 (MQ=255)
        ggCCTGATTACCCTTAGCATCAAGATTTTTTAAGGTGATAGCATCAAGCTTGAGAATGCCGTGATACTGCTTTTTGAATTGATCTTTTAAGATATTTTCAGTTGGCGGTTCGCCATCTACCCAGGCATTACCAGtcatc                                                                                                            <  2:47313/139‑1 (MQ=255)
        ggCCTGATTACCCTTAGCATCAAGATTTTTTAAGGTGATAGCATCAAGCTTGAGAATGCCGTGATACTGCTTTTTGAATTGATCTTTTAAGATATTTTCAGTTGGCGGTTCGCCATCTACCCAGGCATTACCAGtcatc                                                                                                            <  1:239705/139‑1 (MQ=255)
         gCCTGATTACCCTTAGCATCAAGATTTTTTAAGGTGATAGCATCAAGCTTGAGAATGCCGTGATACTGCTTTTTGAATTGATCTTTTAAGATATTTTCAGTTGGCGGTTCGCCATCTACCCAGGCATTACCAGtcatca                                                                                                           >  2:396099/1‑139 (MQ=255)
                        gCATCAAGATTTTTTAAGGTGATAGCATCAAGCTTGAGAATGCCGTGATACTGCTTTTTGAATTGATCTTTTAAGATATTTTCAGTTGGCGGTTCGCCATCTACCCAGGCATTACCAGTCATCAACATAAGCAgtgccg                                                                                            >  2:841570/1‑139 (MQ=255)
                            cAAGATTTTTTAAGGTGATAGCATCAAGCTTGAGAATGCCGTGATACTGCTTTTTGAATTGATCTTTTAAGATATTTTCAGTTGGCGGTTCGCCATCTACCCAGGCATTACCAGTCATCAAc                                                                                                         <  1:358545/122‑1 (MQ=255)
                            cAAGATTTTTTAAGGTGATAGCATCAAGCTTGAGAATGCCGTGATACTGCTTTTTGAATTGATCTTTTAAGATATTTTCAGTTGGCGGTTCGCCATCTACCCAGGCATTACCAGTCATCAAc                                                                                                         >  2:358545/1‑122 (MQ=255)
                             aaGATTTTTTAAGGTGATAGCATCAAGCTTGAGAATGCCGTGATACTGCTTTTTGAATTGATCTTTTAAGATATTTTCAGTTGGCGGTTCGCCATCTACCCAGGCATTACCAGTCATCAACATAAGCAGTGCCGTGCca                                                                                       >  2:717006/1‑139 (MQ=255)
                                             aTAGCATCAAGCTTGAGAATGCCGTGATACTGCTTTTTGAATTGATCTTTTAAGATATTTTCAGTTGGCGGTTCGCCATCTACCCAGGCATTACCAGTCATCAACATAAGCAGTGCCGTGCCAACAATCCATTGCtttt                                                                       >  1:761748/1‑139 (MQ=255)
                                               aGCATCAAGCTTGAGAATGCCGTGATACTGCTTTTTGAATTGATCTTTTAAGATATTTTCAGTTGGCGGTTCGCCATCTACCCAGGCATTACCAGTCATCAACATAAGCAGTGCCGTGCCAACAATCCATTGCTTTTTc                                                                     <  1:841570/139‑1 (MQ=255)
                                                gCATCAAGCTTGAGAATGCCGTGATACTGCTTTTTGAATTGATCTTTTAAGATATTTTCAGTTGGCGGTTCGCCATCTACCCAGGCATTACCAGTCATCAACATAAGCAGTGCCGTGCCAACAATCCATTGCTTTTTCa                                                                    <  2:16015/139‑1 (MQ=255)
                                                            agaATGCCGTGATACTGCTTTTTGAATTGATCTTTTAAGATATTTTCAGTTGGCGGTTCGCCATCTACCCAGGCATTACCAGTCATCAACATAAGCAGTGCCGTGCCAACAATCCATTGCTTTTTCACCCGAAAAACCt                                                        >  2:542395/1‑139 (MQ=255)
                                                                tGCCGTGATACTGCTTTTTGAATTGATCTTTTAAGATATTTTCAGTTGGCGGTTCGCCATCTACCCAGGCATTACCAGTCATCAACATAAGCAGTGCCGTGCCAACAATCCATTGCTTTTTCACCCGAAAAACCTCAcg                                                    <  1:906270/139‑1 (MQ=255)
                                                                tGCCGTGATACTGCTTTTTGAATTGATCTTTTAAGATATTTTCAGTTGGCGGTTCGCCATCTACCCAGGCATTACCAGTCATCAACATAAGCAGTGCCGTGCCAACAATCCATTGCTTTTTCACCCGAAAAACCTCAcg                                                    >  1:128784/1‑139 (MQ=255)
                                                                    gTGATACTGCTTTTTGAATTGATCTTTTAAGATATTTTCAGTTGGCGGTTCGCCATCTACCCAGGCATTACCAGTCATCAACATAAGCAGTGCCGTGCCAACAATCCATTGCTTTTTCACCCGAAAAACCTCACGCGGa                                                <  2:761748/139‑1 (MQ=255)
                                                                       aTACTGCTTTTTGAATTGATCTTTTAAGATATTTTCAGTTGGCGGTTCGCCATCTACCCAGGCATTACCAGTCATCAACAt                                                                                                       >  2:394076/1‑81 (MQ=255)
                                                                       aTACTGCTTTTTGAATTGATCTTTTAAGATATTTTCAGTTGGCGGTTCGCCATCTACCCAGGCATTACCAGTCATCAACAt                                                                                                       <  1:394076/81‑1 (MQ=255)
                                                                           tGCTTTTTGAATTGATCTTTTAAGATATTTTCAGTTGGCGGTTCGCCATCTACCCAGGCATTACCAGTCATCAACTTAAGCAGTGCCGTGCCAACAATCCATTGCTTTTTCACCCGAAAAACCTCACGCGGAATAAttt                                         >  2:746296/1‑139 (MQ=255)
                                                                                tttGAATTGATCTTTTAAGATATTTTCAGTTGGCGGTTCGCCATCTACCCAGGCATTACCAGTCATCAACATAAGCAGTGCCGTGCCAACAATCCATTGCTTTTTCACCCGAAAAACCTCACGCGGAATAATTTATTc                                     <  2:723739/138‑1 (MQ=255)
                                                                                tttGAATTGATCTTTTAAGATATTTTCAGTTGGCGGTTCGCCATCTACCCAGGCATTACCAGTCATCAACATAAGCAGTGCCGTGCCAACAATCCATTGCTTTTTCACCCGAAAAACCTCACGCGGAATAATTTATTc                                     >  1:723739/1‑138 (MQ=255)
                                                                                  tGAATTGATCTTTTAAGATATTTTCAGTTGGCGGTTCGCCATCTACCCAGGCATTACCAGTCATCAACATAAGCAGTGCCGTGCCAACAATCCATTGCTTTTTCACCCGAAAAACCTCACGCGGAATAATTTATTCCgg                                  <  2:725471/139‑1 (MQ=255)
                                                                                  tGAATTGATCTTTTAAGATATTTTCAGTTGGCGGTTCGCCATCTACCCAGGCATTACCAGTCATCAACATAAGCAGTGCCGTGCCAACAATCCATTGCTTTTTCACCCGAAAAACCTCACGCGGAATAATTTATTCCgg                                  <  1:717006/139‑1 (MQ=255)
                                                                                      ttGATCTTTTAAGATATTTTCAGTTGGCGGTTCGCCATCTACCCAGGCATTACCAGTCATCAACATAAGCAGTGCCGTGCCAACAATCCATTGCTTTTTCACCCGAAAAACCTCACGCGGAATAATTTATTCCGGTGAt                              <  2:128784/139‑1 (MQ=255)
                                                                                      ttGATCTTTTAAGATATTTTCAGTTGGCGGTTCGCCATCTACCCAGGCATTACCAGTCATCAACATAAGCAGTGCCGTGCCAACAATCCATTGCTTTTTCACCCGAAAAACCTCACGCGGAATAATTTATTCCGGTGAt                              <  1:396099/139‑1 (MQ=255)
                                                                                            ttttAAGATATTTTCAGTTGGCGGTTCGCCATCTACCCAGGCATTACCAGTCATCAACTTAAGCAGTGCCGTGCCAACAATCCATTGCTTTTTCACCCGAAAAACCTCACGCGGAATAATTTATTCCGGTGATTATCAt                        <  1:746296/139‑1 (MQ=255)
                                                                                              ttAAGATATTTTCAGTTGGCGGTTCGCCATCTACCCAGGCATTACCAGTCATCAACATAAGCAGTGCCGTGCCAACAATCCATTGCTTTTTCACCCGAAAAACCTCACGCGGAATAATTTATTCCGGTGATTATCATTa                      >  2:607623/1‑139 (MQ=255)
                                                                                                     atTTTCAGTTGGCGGTTCGCCATCTACCCAGGCATTAcc                                                                                                                   >  1:629697/1‑39 (MQ=255)
                                                                                                     atTTTCAGTTGGCGGTTCGCCATCTACCCAGGCATTAcc                                                                                                                   <  2:629697/39‑1 (MQ=255)
                                                                                                        ttCAGTTGGCGGTTCGCCGTCTACCCAGGCATTACCAGTCATCAACATAAGCAGTGCCGTGCCAACAATCCATTGCTTTTTCACCCGAAAAACCTCACGCGGAATAATTTATTCCGGTGATTATCATTAGGGTAATTaa            >  1:507073/1‑139 (MQ=255)
                                                                                                        ttCAGTTGGCGGTTCGCCATCTACCCAGGCATTACCAGtca                                                                                                              <  1:5788/41‑1 (MQ=255)
                                                                                                        ttCAGTTGGCGGTTCGCCATCTACCCAGGCATTACCAGtca                                                                                                              >  2:5788/1‑41 (MQ=255)
                                                                                                        ttCAGTTGGCGGTTCGCCATCTACCCAGGCATTACCAGTCATCAACATAAGCAGTGCCGTGCCAACAATCCATTGCTTTTTCACCCGAAAAACCTCACGCGGa                                                >  2:853861/1‑103 (MQ=255)
                                                                                                        ttCAGTTGGCGGTTCGCCATCTACCCAGGCATTACCAGTCATCAACATAAGCAGTGCCGTGCCAACAATCCATTGCTTTTTCACCCGAAAAACCTCACGCGGa                                                <  1:853861/103‑1 (MQ=255)
                                                                                                          cAGTTGGCTGTTCGCCATCTACCCAGGCATTACCAGTCATCAACATAAGCAGTGCCGTGCCAACAATCCATTGCTTTTTCACCCGAAAAACCTCACGCGGAATAATTTATTCCGGTGATTATCATTAGGGTAATTaaaa          <  1:542395/139‑1 (MQ=255)
                                                                                                             ttGGCGGTTCGCCATCTACCCAGGCATTACCAGTCATCAACATAAGCAGTGCCGTGCCAACAATCCATTGCTTTTTCACCCGAAAAACCTCACGCGGAATAATTTATTCCGGTGAt                              >  2:520559/1‑116 (MQ=255)
                                                                                                             ttGGCGGTTCGCCATCTACCCAGGCATTACCAGTCATCAACATAAGCAGTGCCGTGCCAACAATCCATTGCTTTTTCACCCGAAAAACCTCACGCGGAATAATTTATTCCGGTGAt                              <  1:520559/116‑1 (MQ=255)
                                                                                                               ggCGGTTCGCCATCTACCCAGGCATTACCAGTCATCAACATAAGCAGTGCCGTGCCAACAATCCATTGCTTTTTCACCCGAAAAACCTCACGCGGAATAATTTATTc                                     <  1:352300/107‑1 (MQ=255)
                                                                                                               ggCGGTTCGCCATCTACCCAGGCATTACCAGTCATCAACATAAGCAGTGCCGTGCCAACAATCCATTGCTTTTTCACCCGAAAAACCTCACGCGGAATAATTTATTc                                     >  2:352300/1‑107 (MQ=255)
                                                                                                                gCGGTTCGCCATCTACCCAGGCATTACCAGTCATCAACATAAGCAGTGCCGTGCCAACAATCCATTGCTTTTTCACCCGAAAAACCTCACGCGGAATAATTTATTCCGGTGAt                              <  2:211534/113‑1 (MQ=255)
                                                                                                                gCGGTTCGCCATCTACCCAGGCATTACCAGTCATCAACATAAGCAGTGCCGTGCCAACAATCCATTGCTTTTTCACCCGAAAAACCTCACGCGGAATAATTTATTCCGGTGAt                              >  1:211534/1‑113 (MQ=255)
                                                                                                                  ggTTCGCCATCTACCCAGGCATTACCAGTCATCAACATAAGCAGTGCCGTGCCAACAATCCATTGCTTTTTCACCCGAAAAACCTCACGCGGAATAATTTATTCCGGTGATTATCATTAGGGTAATTAAAATTGGCAAt  <  2:507073/139‑1 (MQ=255)
                                                                                                                              |                                                                                                                              
GACCAGGTGGCCTGATTACCCTTAGCATCAAGATTTTTTAAGGTGATAGCATCAAGCTTGAGAATGCCGTGATACTGCTTTTTGAATTGATCTTTTAAGATATTTTCAGTTGGCGGTTCGCCATCTGCCCAGGCATTACCAGTCATCAACATAAGCAGTGCCGTGCCAACAATCCATTGCTTTTTCACCCGAAAAACCTCACGCGGAATAATTTATTCCGGTGATTATCATTAGGGTAATTAAAATTGGCAAT  >  NC_000913/3099378‑3099630

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