Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,452,057 A→T I77N (ATC→AAC)  rplD ← 50S ribosomal subunit protein L4

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,452,0570AT93.8% 213.2 / 3.5 64I77N (ATC→AAC) rplD50S ribosomal subunit protein L4
Reads supporting (aligned to +/- strand):  ref base A (3/1);  new base T (27/33);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 3.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

GGATGCTTTTCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGATCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGTCTTCTG  >  NC_000913/3451962‑3452151
                                                                                               |                                                                                              
ggATGCTTTTCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCgg                                                                                            >  1:240760/1‑100 (MQ=255)
 gATGCTTTTCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCggg                                                                                           <  1:1121272/100‑1 (MQ=255)
 gATGCTTTTCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCggg                                                                                           >  1:507850/1‑100 (MQ=255)
 gATGCTTTTCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCggg                                                                                           <  1:662608/100‑1 (MQ=255)
 gATGCTTTTCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCggg                                                                                           >  1:1366123/1‑100 (MQ=255)
   tGCTTTTCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGct                                                                                         >  1:93976/1‑100 (MQ=255)
   tGCTTTTCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGc                                                                                          >  1:248290/1‑99 (MQ=255)
     cTTTTCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGctct                                                                                       >  2:1389063/1‑100 (MQ=255)
     cTTTTCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGctct                                                                                       >  1:384066/1‑100 (MQ=255)
      ttttCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCtt                                                                                      >  2:693688/1‑100 (MQ=255)
      ttttCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCtt                                                                                      <  1:1557786/100‑1 (MQ=255)
          cAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATa                                                                                  <  2:567000/100‑1 (MQ=255)
          cAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATa                                                                                  <  1:241429/100‑1 (MQ=255)
            gcgcgcCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGa                                                                                <  1:46663/100‑1 (MQ=255)
                 ccgcGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAg                                                                           <  1:275116/100‑1 (MQ=255)
                  cgcgGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGa                                                                          <  1:316413/100‑1 (MQ=255)
                           tcttctTGTTAACTTTTTGACTGTGGTTCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGATCGGGCTCTTGATAGAACCAGAAcgcg                                                                     >  2:749415/1‑96 (MQ=255)
                            cttcttGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACgg                                                                >  1:751856/1‑100 (MQ=255)
                                 tGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGt                                                           >  1:177573/1‑100 (MQ=255)
                                 tGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGt                                                           >  1:378018/1‑100 (MQ=255)
                                  gTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTg                                                          <  1:1210339/100‑1 (MQ=255)
                                  gTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTg                                                          >  1:652994/1‑100 (MQ=255)
                                    tAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGcc                                                        >  1:302504/1‑100 (MQ=255)
                                     aaCTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGATCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCt                                                       >  2:1248520/1‑100 (MQ=255)
                                       cTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCttt                                                     <  1:1257960/100‑1 (MQ=255)
                                          tttGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTg                                                  <  1:1510556/100‑1 (MQ=255)
                                           tccgatctgGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTgg                                                 <  2:852447/93‑1 (MQ=255)
                                            tGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGgc                                                <  1:1467893/100‑1 (MQ=255)
                                             gACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGgcg                                               <  2:1483700/100‑1 (MQ=255)
                                             gACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGgcg                                               <  2:1392976/100‑1 (MQ=255)
                                               cTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGcc                                             >  2:977117/1‑100 (MQ=255)
                                               cTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGcc                                             >  2:1439106/1‑100 (MQ=255)
                                               cTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGcc                                             <  2:379881/100‑1 (MQ=255)
                                                   ggTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTggagatcgga                                         >  1:852453/1‑92 (MQ=255)
                                                     tCCTGCGAACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGt                                        <  2:174728/99‑1 (MQ=255)
                                                      ccTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGttt                                      <  1:334714/100‑1 (MQ=255)
                                                       cTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGtttt                                     >  1:980622/1‑100 (MQ=255)
                                                         acGGACGAGCAACAAAGATTACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCCCGGCCGGATCCTTTCTGGCGCCACGGtttttt                                   >  2:1501712/2‑100 (MQ=255)
                                                           ggACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTAc                                 >  1:698998/1‑100 (MQ=255)
                                                           ggACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTAc                                 >  1:504738/1‑100 (MQ=255)
                                                           ggACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTAc                                 <  1:252382/100‑1 (MQ=255)
                                                           ggACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTAc                                 <  1:28848/100‑1 (MQ=255)
                                                             aCGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCg                               <  1:1194116/100‑1 (MQ=255)
                                                              cGAGCAGCAAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCgg                              <  1:348713/100‑1 (MQ=255)
                                                                  cagcaAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAAcc                          <  1:884822/100‑1 (MQ=255)
                                                                    gcaAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAg                        >  2:1357917/1‑100 (MQ=255)
                                                                    gcaAAGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAg                        >  1:1238718/1‑100 (MQ=255)
                                                                      aaaGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGtt                      <  1:334155/100‑1 (MQ=255)
                                                                        aGGTCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTAc                    <  2:1367662/100‑1 (MQ=255)
                                                                        aGGTCACGCCACCAGAACGCCAGATCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTAc                    <  1:1248528/100‑1 (MQ=255)
                                                                           tCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTc                 >  1:1304162/1‑100 (MQ=255)
                                                                           tCACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTc                 <  2:507847/100‑1 (MQ=255)
                                                                            cACGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTc                 >  2:1515646/1‑99 (MQ=255)
                                                                             aCGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCa                >  2:1026365/1‑99 (MQ=255)
                                                                              cGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGc              <  1:62875/100‑1 (MQ=255)
                                                                              cGCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGc              <  1:693691/100‑1 (MQ=255)
                                                                              cGCCACCAGAACGCCAGATCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGc              >  2:791746/1‑100 (MQ=255)
                                                                               gCCACCAGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCa             >  1:1467285/1‑100 (MQ=255)
                                                                                ccaccaGAACGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCAc            <  2:352360/100‑1 (MQ=255)
                                                                                        aCGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGtcttc    <  1:1480375/100‑1 (MQ=255)
                                                                                        aCGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGtcttc    <  2:525697/100‑1 (MQ=255)
                                                                                        aCGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGtcttc    >  1:1320931/1‑100 (MQ=255)
                                                                                        aCGCCAGTTCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGtcttc    <  1:1501722/100‑1 (MQ=255)
                                                                                          gCCAGTTAGGGCTCCTGATACAACCAGACCGCGCAAGGCCGGTGCCTTTCTGCCGCCACGGTTTTTTACCGGAACCAGTTCCTTCAGCCCGAGTCTTCTg  <  2:285675/100‑1 (MQ=255)
                                                                                               |                                                                                              
GGATGCTTTTCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGATCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGTCTTCTG  >  NC_000913/3451962‑3452151

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