Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 905,720 C→A L270M (CTG→ATG)  amiD → 1,6‑anhydro‑N‑acetylmuramyl‑L‑alanine amidase, Zn‑dependent; OM lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913905,7200CA94.6% 183.2 / 3.1 56L270M (CTG→ATG) amiD1,6‑anhydro‑N‑acetylmuramyl‑L‑alanine amidase, Zn‑dependent; OM lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base A (30/23);  total (32/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.96e-01

CAGCGGCGCGTGATTATGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGCTGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAGGGTGATGACCGG  >  NC_000913/905624‑905814
                                                                                                |                                                                                              
cagcGGCGCGTGATTATGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATg                                                                                              >  2:368378/1‑99 (MQ=255)
  gcggcgCGTGATTATGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGgag                                                                                           <  2:525919/100‑1 (MQ=255)
  gcggcgCGTGATTATGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGgag                                                                                           >  1:607961/1‑100 (MQ=255)
  gcggcgCGTGATTATGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGgag                                                                                           <  2:1260249/100‑1 (MQ=255)
     gcgcGTGATTATGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGaaa                                                                                        <  1:1058989/100‑1 (MQ=255)
      cgcgTGATTATGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAAt                                                                                       <  2:1448355/100‑1 (MQ=255)
          tGATTATGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACgg                                                                                   <  2:352083/100‑1 (MQ=255)
               aTGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAgg                                                                              >  1:1515351/1‑100 (MQ=255)
               aTGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAgg                                                                              >  2:1309876/1‑100 (MQ=255)
                tGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAgg                                                                              >  1:1391386/1‑99 (MQ=255)
                tGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGa                                                                             >  1:1412022/1‑100 (MQ=255)
                tGGCATGCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGa                                                                             >  2:313199/1‑100 (MQ=255)
                   cATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTa                                                                          <  2:571939/100‑1 (MQ=255)
                    aTTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAg                                                                         >  2:549203/1‑100 (MQ=255)
                      tCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAgcg                                                                       <  1:238886/100‑1 (MQ=255)
                      tCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAgcg                                                                       >  1:1382161/1‑100 (MQ=255)
                       ccAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAgcgc                                                                      <  1:590845/100‑1 (MQ=255)
                         aGATGCATTTCCGCCCGCCGTTATATAACGGCGACGCGGATGAAGACACTCCGGCGATTGCCGAAGCATTGAAGGAGAAATACGGGCAGGATTCGCGCgg                                                                    >  2:294158/1‑100 (MQ=255)
                          gATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGAGGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGc                                                                   <  2:384372/100‑1 (MQ=255)
                            tGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAg                                                                 <  1:594624/100‑1 (MQ=255)
                                tttCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGtttt                                                             >  1:954856/1‑100 (MQ=255)
                                 ttCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTc                                                            <  2:1035319/100‑1 (MQ=255)
                                  tCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTcc                                                           >  2:576345/1‑100 (MQ=255)
                                   ccGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCg                                                          <  1:632099/100‑1 (MQ=255)
                                    cGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCg                                                          <  2:88168/99‑1 (MQ=255)
                                      cccGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTgg                                                       <  1:1417499/100‑1 (MQ=255)
                                                 ataACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCa                                            >  2:1163700/1‑100 (MQ=255)
                                                 ataACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCa                                            >  1:701276/1‑100 (MQ=255)
                                                  taACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCTCGTAGCCAg                                           >  2:1317748/1‑100 (MQ=255)
                                                  taACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAg                                           <  1:31347/100‑1 (MQ=255)
                                                   aaCGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAgg                                          >  2:1221502/1‑100 (MQ=255)
                                                   aaCGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGAGCGGGAGTTTTCCGTGGTCGCGTAGCCAgg                                          >  2:148390/1‑100 (MQ=255)
                                                     cGGCGACGCGGATGCAGAAACCAAGACGATTGCCGAAGCATTGATGGAGACATACGGGCAGGAATAGCGCGGCCGTTTTCCGTGGTAGCGCAGCCAg                                           >  2:831097/1‑97 (MQ=255)
                                                     cGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAggcg                                        <  2:333498/100‑1 (MQ=255)
                                                            gaggATGCGAAAACTCAGGCGATGGCCGAGGAATTGCTGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGTAGttt                                 <  2:569773/98‑1 (MQ=255)
                                                              ggATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTtc                               <  1:375820/100‑1 (MQ=255)
                                                               gATGCAGAAATTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTtct                              <  2:713062/100‑1 (MQ=255)
                                                               gATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTtct                              <  1:1221507/100‑1 (MQ=255)
                                                                 tGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCg                            >  1:299027/1‑100 (MQ=255)
                                                                  gCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGa                           <  1:1392942/100‑1 (MQ=255)
                                                                   cAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGAt                          <  2:1412014/100‑1 (MQ=255)
                                                                   cAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGAt                          >  2:1381443/1‑100 (MQ=255)
                                                                   cAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGAt                          >  2:1546334/1‑100 (MQ=255)
                                                                   cAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGAt                          >  1:1207574/1‑100 (MQ=255)
                                                                    aGAAACTCAGGCGATTGCCGAAGCATTGCTGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATa                         >  2:1577393/1‑100 (MQ=255)
                                                                        aCTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCtt                     >  1:399636/1‑100 (MQ=255)
                                                                          tCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCa                   <  2:1449555/100‑1 (MQ=255)
                                                                           cAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCAt                  >  2:941522/1‑100 (MQ=255)
                                                                           cAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCAt                  >  2:272089/1‑100 (MQ=255)
                                                                                   tGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAGGGtg          >  1:1286699/1‑100 (MQ=255)
                                                                                    gCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAGGGtga         >  2:297544/1‑100 (MQ=255)
                                                                                     ccGAAGCATTGCTGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAGGGtgat        >  1:453328/1‑100 (MQ=255)
                                                                                     ccGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAGGGtgat        >  2:1292549/1‑100 (MQ=255)
                                                                                     ccGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAGGGtgat        <  1:1494655/100‑1 (MQ=255)
                                                                                      cGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAGGGtgatg       >  1:574457/1‑100 (MQ=255)
                                                                                           cATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGTTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATTCAGGGTTCTGACCgg  >  2:482010/1‑100 (MQ=255)
                                                                                                |                                                                                              
CAGCGGCGCGTGATTATGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGCTGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAGGGTGATGACCGG  >  NC_000913/905624‑905814

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