Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,884,434 T→A 9.8% R33R (CGT→CGA)  rpmH → 50S ribosomal subunit protein L34

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,884,4340TA 9.8% 105.8 / 5.2 48R33R (CGT→CGA) rpmH50S ribosomal subunit protein L34
Reads supporting (aligned to +/- strand):  ref base T (15/28);  new base A (3/2);  total (18/30)
Fisher's exact test for biased strand distribution p-value = 3.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.33e-01

AGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCACTTTTCAACCGTCTGTACTGAAGCGCAACCGTTCTCACGGCTTCCGTGCTCGTATGGCTACTAAAAATGGTCGTCAGGTTCTGGCACGTCGTCGTGCTAAAGGCCGCGCTCGTCTGACCGTTTCTAAGTAATAAAGCTAACCCCTGAGTGGTTAAGCTCGCATTTCCCAGGGAGTTACGCTTGTTAACTCCCAGTCAATTCACATTCGTCTTCCAGCAGCCACAACGGGCT  >  NC_000913/3884299‑3884576
                                                                                                                                       |                                                                                                                                              
aGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCACTTTTCAACCGTCTGTACTGAAGCGCAACCGTTCTCACGGCTTCCGTGCTCGTATGGCTACTAAAAATGGTCGTCAGGTTCTGGCACGTCGTCGTg                                                                                                                                         <  1:669945/143‑1 (MQ=255)
   cATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCACTTTTCAACCGTCTGTACTGAAGCGCAACCGTTCTCACGGCTTCCGTGCTCGTATGGCTACTAAAAATGGTCGTCAGGTTCTGGCACGTCGTCGTGCTa                                                                                                                                      >  1:70499/1‑143 (MQ=255)
     ttGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCACTTTTCATCCGTCTGTACTGAAGCGCAACCGTTCTCACGGCTTCCTTGCTCGTATGGCCACTAAAAATGGTCGTCAGGTTCTGGCACGTCGTCGTGCTaaa                                                                                                                                    <  2:676070/143‑1 (MQ=255)
     ttGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCACTTTTCAACCGTCTGTACTGAAGCGCAACCGTTCTCACGGCTTCCGTGCTCGTATGGCTACTAAAAATGGGCGTCAGGTTCTGGCACGTCGTCGTgatcca                                                                                                                                    >  2:759542/1‑138 (MQ=255)
       gAATTTCTTCAAGTTTAGGTAGAAATCGCCATGAAACGCACTTTTCAACCGTCTGTACTGAAGCGCAACCGTTCTCACGGCTTCCGTGCTCGTATGGCTACTAAAAATGGTCGTCAGGTTCTGGCACGACGACGTGccagaac                                                                                                                                  <  2:154780/143‑7 (MQ=255)
       gAACTTAGTCAAGGTTGGGGAGACATCGCCATGAAGCGCACTTTTCAACCGTCTGTACTGAAGCGCAACCGTTCTCACGGCGTCCGTGCTCGTATGGCTACTAAAAATGGTCGTCAGGTTCTGGCACGACGACGTGccagaac                                                                                                                                  <  2:841792/143‑7 (MQ=255)
           ttaagcAAGTTTAAGTAGAAATCGCCATGAAACGCACTTTTCAACCGTCTGTACTGAAGCGCAACCGTTCTCACAGCTTCCGTGCTCGTATGGCTACTAAAAATGGTCGTCAGGTTCTGGCACGTCGTCGTGCTAAAGGCcgc                                                                                                                              <  2:102243/138‑1 (MQ=255)
           ttATTCAAGTTTAGGTAGAAATCGCCATGAAACGCACTTCTCAACCGTCTGTACTGAAGCGCAACCGTTCTCACGGCTTCCGTGCTCGTATGGCTACTAAAAATGGTCGTCAGGTTCTGGCACGTCGTCGTGCTAAAGGCcgc                                                                                                                              <  1:147416/143‑1 (MQ=255)
              ttCAAGTTTAGGTAGAAATCGCCATGACACACACTTTTCAACCGTCTGTACTGACGCGCAACCGTTCTCACGGCTTCCGTGCTCGTATGTCTCCTATAAATCGTCGTCAAGATCAGACACATCGTAGTGCTAAAGACAGCGCt                                                                                                                           >  1:38093/1‑143 (MQ=255)
               tCAAGTTTAGGTAGAAATCGCCATGAAACGCACTTTTCAACCGTCTGTACTGAAGCGCAACCGTTCTCACGGCTTCCGTGCTCGTATGGCTACTAAAAATGGTCGTCAGGTTCTGGCACGTCGTAGTGCTAAAGGCCGCGCtg                                                                                                                          >  2:606286/1‑142 (MQ=255)
                          tAGAAATCGCCATGAAACGCACTTTTCAACCGTCTGTACTGAAGCACAACCGTTCTCACGGCTTCCGTGCTCGTATGGCTACTAAAAATGGTCGTCAGGTTCTGGCACGTCGTCGTGCTAAAGGCCGCGCTCGTCTGACCGtt                                                                                                               >  1:300109/1‑143 (MQ=255)
                            gAAATCGCCATGAACCGCCCTTTTCAACCGTCCGTACTGAAGCGCAACCGTTCTCACGGCTTCCGTGCTCGTATGGCTACTAAAAATGGTCGTCAGGTTCTGGCACGTCGTCGTGCTAAAGGCCGCGCTCGTCTGACCGTTTc                                                                                                             <  2:486900/143‑1 (MQ=255)
                               atggCCATGGAACGCACTTTCCAACCGTCTGTACTGAAGCGCAGCCGTTCTCACGGCTTCCGGGCTCGTATGGCTACTAAAAGTGGTCGTCAGGTTCTGGCACGTCGTCGTGCTAAAGGCCGCGCTCGTCTGACCGTTTCTaa                                                                                                          <  1:27788/140‑1 (MQ=255)
                                    cATGAAACGCACTTTGCAACCGTCTGTACTGAAGCGCAACCGTTCGCAAGGCTTCCGTGCTCGTATGGCTACTAAAAATGGTCGTCAGGTTCTGGCACGTGGTCGTGCTAAAGACCGCGCTCGTGTGAACGTTTCTAAGgaat                                                                                                     >  2:240257/1‑143 (MQ=255)
                                      tGACACGCACTTGTGAACCGTCTTTACTGAAGCGCAACGGTTGTCACGTCTTCCGTGCTCGTATGGCTACTAAAAATGGTGGTCAGGTTCTGGCACGTCGTCGTGCTAAAGGCCGCGCTCGTCTGACCGTTTCTAAGtaataa                                                                                                   <  2:32835/143‑1 (MQ=255)
                                        acaCGCACTTTTCAGCCGGCTGGAGTGATGCGCCTCCGTTCGCACGGCTTCCGTGGTCGTATGGCTACTCAAAATGGTCGTCAGGTTCTGGCACGTCGTCGTGCTAAAGGCCGCGCTCGTCTGACCGTTTCTAAGTAATAAAg                                                                                                 <  2:734813/141‑1 (MQ=255)
                                        aaaCGCACTTTTCAACCGTCTGTACTGAAGCGCAACCGTTCTCACGGCTTCCGTGCTCGTATGGCTACTAAAAATGGTCGTCAGGTTCTGGCACGTCGTCGTGCTAAAGGCCGCGCTCGTCTGACCGTTTCTAAGTAATAAAg                                                                                                 <  1:680244/143‑1 (MQ=255)
                                        aaaCGCACGTTTCAACCGTCTGTACTGAAGCGCAACCGTTCTCACGGCTTCCGTGCTCGTATGGCTACTAAAAATGGTCGTCAGGTTCTGGCACGTCGTCGTGCTAAAGGCCGCGCTCGTCTGACCGTTTCTAAGTAATAAAg                                                                                                 <  1:582318/143‑1 (MQ=255)
                                            gCACTTTTCAACCGTCTGTACTGAAGCGCAACCGTTCTCACGGCTTCCGTCGTCGTATCCCTACCAAAAATGACCCTCCGGTCCCGGCACGAGGCCCAGCTAAAGGCCGCACCCGTCCGACCGTTGTCAAGTAAAAAAGCaaa                                                                                             >  2:379018/1‑143 (MQ=255)
                                               cTTTTCACACTTCTGTAGTGAGGCGCATCCGTTCGCTTGGCTGCCGTGCTCGTATGGCTACTAAGAATGGTCGTCCGGTTCTGGCACGTCGTCGTGCTAAAGGCCGCGCTCGTCTGACCGTTTCTAAGTAATAAAGCTAAccc                                                                                          <  2:301657/143‑1 (MQ=255)
                                                 tttCAACCGTCTGTACTGAAGCGCAACCGTTCTCACGGCTTCCGTGCTCGTATGGCTACTAAAAATGGTCGTCAGGTTCTGGCACGTCGTCGTGCTAAAGGCCGCGCTCGTCTGACCGTTTCTAAGTAATAAAGCTAACCCCt                                                                                        <  1:509860/143‑1 (MQ=255)
                                                              tACTGAAGCGCAACCGTTCTCACGGCTTCCGTGCTCGTATGGCTACTAAGAATGCTCGTCAGGTTTTGGCCTCTTGACGTGGTAAAGAAGGTATACGTCTCACCGTGGCGACGTACTAAAGCTAACCCCGGAGTggcgagtct                                                                           >  2:888996/1‑136 (MQ=255)
                                                              tACTGAAGCGCAACCGTTCTCACGGCTTCCGTGCTCGTATGGCTACTAAAAATGGTCGTCAGGTTCTGGCACGTCGTCGTGCTAAAGGCCGCGCTCGTCTGACCGATACTAAGTAATAAAGCTAACCCCTGAGTGGTTAAGCt                                                                           >  1:103097/1‑143 (MQ=255)
                                                                 tGAAGCGCAACCGTTCTCACGGCTTCCGTGCTCGTATGGCTACTAAAAATGGTCGTCAGGTTCTGGCACGTCCTCGTGCTAAAGGCCGCGCTCGTCAGACCGTGTGTAAGTACTAACGCTAACGCCTGACTTGTTAAGCACGc                                                                        >  2:120805/1‑143 (MQ=255)
                                                                        cAACCGTTCTCACGGCTTCCGTGCTCGTATGGCTACTAAAAATGGTCGTCAGGTTCTGGCACGTCGTCGTGCTAAAGGCCGCGCTCGTCTGACCGTTTCTAAGTAATAAAGCTAACCCCTGAGTGGTTAAGCTCGCATTTccc                                                                 >  2:251145/1‑143 (MQ=255)
                                                                             gTTCTCACGGCTTCCGTGCTCGTATGGCTACTAAAAATGGTCGTCCGGTTCTGGCACGTCGTCGTGCTAACGGCCGCGCTCGTCTGACCGTTACTAAGTCGGCAAGCTAGCCCCGCAGCGGCTAAGCCCGCCTCTCCCaccgc                                                            >  1:166201/1‑139 (MQ=255)
                                                                                 tcACGGCTTCCGTGCTCGTATGGCTACTAAAAATGGTCGTCAGGTTCTGGCACGTCGTCGTGCTAAAGGCCGCGCTCGTCTGACCGTTTCTAAGTAATAAAGCTAACCCCTGAGTGGTTAAGCTCGCATTTCCCAGGGAGTTa                                                        <  2:253919/143‑1 (MQ=255)
                                                                                   aCGGCTTCCTTGCTCGTATGGCTACTAAAAATGGTCGTCAGGTTCTGGCACGTCGTCGTGCTAAAGCCCGCGCTCGTCTGACCGTTTCTAAGTAATAAAGCTAACCCCTGAGTGGTTACGCTCGCATTTCCCAGGGAGATAcc                                                      >  2:449198/1‑142 (MQ=255)
                                                                                   aCGGCTTCCGTGCTCGTATGGCTACTAAAAATGGTCGTCAGGTTCTGGCACGTCGTCGTGCTAAAGGCCGCGCTCGTCTGACCGTTTCTAAGTAATAAAGCTAACCCCTGAGTGGTTAAGCTCGCATTTCCCAGGGAGTTACg                                                      >  1:774271/1‑143 (MQ=255)
                                                                                    cGGCTTCCGTGCTCGTATGGCTACTAAAAATGGTCGTCAGGTTCTGGCACGTCGTCGTGCTAAAGGCCGCGCTCGTCTGACCGTTTCTAAGTAATAAAGCTAACCCCTGAGTGGTTAAGCTCGCATTTCCCAGGGAGTTACGc                                                     <  2:461330/143‑1 (MQ=255)
                                                                                     ggCTTCCGTGCTCGCATGGCTACTAAAAATGGTCGTCAGGTTCGGGCACGTCGTCGTGCTAAAGGCCGCGCTCGTCTGACCGTTTCTAAGTAATAAAGCTAACCCCTGAGTGGTTAAGCTCGCATTTCCCAGGGAGTTACGCt                                                    <  1:759542/143‑1 (MQ=255)
                                                                                      gCTTCCGTGCTCTTATGGCTACTAAAAATGGTCGTCAGGTTCTGGCACGTCGTCGTGCTAAAGGCCGCGCTAATACTACCGTTTCTAAGTAATAAAGCTAACCCCTGAGTGGTTAAGCTCGCATTTCCCAGGGAGTTACGCtt                                                   <  1:818406/143‑1 (MQ=255)
                                                                                        ttCCGGTCTCTTATGGCTACTAAAAATGGTCGTCAGGTTCTTGCACGTCGTCGTGCTAAAGGCCGCGCTCGTCTGACCGTTTCTAAGTAATAAAGCTAATCCCTGAGTGGTTTAGCTCGCATTTCCCAGGGAGTTACGCTTGt                                                 <  2:110848/143‑1 (MQ=255)
                                                                                        tgccagccTCGTATGCCTACTAATAATGGTCGTCAGGTTCTGGCACGTCGGCATGCTAAAGGCCGCGCTCGTCTGACCGTTTCTAAGTAATAAAGCTCACCCCTTAGTGGTTAAGCTCGCATTTCCCAGGGAGTTACGCTTGt                                                 <  2:588698/136‑1 (MQ=255)
                                                                                             tGCTCGTATGGGTACTAAAAATGGTGGTCAGGTTCTGGCCCGTCGTCGTGCTAAAGGCCGCGCTCGTCTGACCGTTTCTAAGTAATAAAGCTAACCCCTGAGTGGTTAAGCTCGCATTTCCCAGGGAGTTACGCTTGTTAACt                                            <  2:824619/143‑1 (MQ=255)
                                                                                                 cggaTGGCTACTAAAAATGGTCGTCAGGTTCTGGCACGTCGTCGTGCTAAAGGCCGCGCTCGTCTGACCGTTTCTAAGTAATAAAGCTAACCCCTGAGTGGTTAAGCTCGCATTTCCCAGGGAGTTACGCTTGTTAACTCCCa                                        <  2:85634/140‑1 (MQ=255)
                                                                                                 cGTATGGCTACTAAAAATGGTCGTCAGGTTCTGGCACGTCGTCGTGCTAAAGGCCGCGCTCGTCTGACCGTTTCTAAGTAATAAAGCTAACCCCTGAGTGGTTAAGCTCGCATTTCCCAGGGAGTTACGCTTGTTAACTCCCa                                        <  1:739503/143‑1 (MQ=255)
                                                                                                 cGTATGGCTACTAAAAATGGTCGTCAGGTTCTGGCACGTCGTCGTGCTAAAGGCCGCGCTCGTCTGACCGTTTCTAAGTAATAAAGCTAACCCCTGAGTGGTTAAGCTCGCATTTCCCAGGGAGTTACGCTTGTTAACTCCCa                                        <  1:160247/143‑1 (MQ=255)
                                                                                                        ctgctaAAAATGGTCGTCAGGTTCTGGCACGTCGTCGTGCTAAAGGCCGCGCTCGTCTGACCGTTTCTAAGTAATAAAGCTAACCCCTGAGTGGTTAAGCTCGCATTTCCCAGGGAGTTACGCTTGTTAACTCCCAGTCAAtt                                 <  1:532210/140‑1 (MQ=255)
                                                                                                         tactaAAAATGGTCGTCAGGTTCTGGCACGTCGTCACGCTAAAGGCCGCGCTCGTCTGACCGTTTCTAAGTAATAAAGCTAACCCCTGAGTGGTTAAGCTCGCATTTCCCAGGGAGTTACGCTTGTTAACTCCCAGTCAATTc                                <  1:812559/143‑1 (MQ=255)
                                                                                                         gcctaCAAAGGGTGGGCAGGTGTTGGCACGTCGTATTATTCAACGCTTCTCTCGTCTGACCGTTTCTAAGTAATAACGCTAACCCCTGAGTGGTTAAGCTCGCATTTCCCAGGGAGTTACGCTTGTTAACTCCCAGTCAATTc                                <  2:399460/141‑1 (MQ=255)
                                                                                                            taACAATGTTCGTCAGGTTCTGGCACGTCGTCGTGCTAAAGGCCGCACTCGTCTGACCGTTCCTAAGTAATAAAGCTAACCCCTCAGTGGTTAAGCTCGCATTTCCCAGAGACTTACGCTCGTTGACCACCAGGCAATTcacc                             >  2:852991/1‑142 (MQ=255)
                                                                                                            taAAAATGGTCGTCAGGTTGTGGCACGTCGTCGTGCTAAAGGCCGCGCTCGTCTGACCGTCTCTGAGTAATAAAGCTAACCCCTGAGTGGTTAAGCTCGCATTTCCCAGGGCGTTACGCTTGTTAACTCCCAGTCAATTcacc                             >  2:196039/1‑142 (MQ=255)
                                                                                                              acaagggTCGTCAGGTTCTGGCACGTCGTCGTGCGAAAGGCCGCGCTCGTCTGACCGTTTCTAAGTAATAAAGCTAACCCCTGAGTGGTTAAGCTCGCATTTCCCAGGGAGTTACGCTTGTTAACTCCCAGTCAATTCACAtt                           <  2:573178/138‑1 (MQ=255)
                                                                                                                       gtcAGGTTCTGGCACGTCGTCGTGCTAAAGGCCGCGCTCGTCTGACCGTTTCTAAGTAATAAAGCTAACCCCTGAGTGGTTAAGCTCGCATTTCCCAGGGAGTTACGCTTGTTAACTCCCAGTCAATTCACATTCGTCTTCca                  <  1:649222/143‑1 (MQ=255)
                                                                                                                             ctttagCACGACGACGTGCTAAAGGCCGCGCCCGTCTGTCCGCTACTAAGTAATAAAGCTAACCCCTGAGGGGTTAAGCGCGCATTTCCCAGGGAGGTACGCTTGTTACCTCCGAGTCAATGCACATTCGTCTTCCAGCAGcc            >  2:650337/6‑143 (MQ=255)
                                                                                                                              tttagCACGACGACGTGCTAAAGGCCGCGCTCGTCTGACCGTTTCTACGTAATAAAGCTAACCCCTGAGTGGTTAAGCTCGCATTTCCCAGGGAGTTACGCTTGTTAACTCCCAGTCAATTCACATTCGTCTTCCAGCAGCca           >  1:596878/5‑143 (MQ=255)
                                                                                                                                    aCGTCGTCGTGCTAGAGGCCGCGCTCGTCTGCCCGTTTCTAAGTAATAAAGCTAACCCCTGAGTGGTTAAGCTCGCATTTCCCAGGGAGTTACGCTTGTTAACTCCCAGTCAATTCACATTCGTCTTCCAGCAGCCACAACgg     <  2:807473/143‑1 (MQ=255)
                                                                                                                                       tcgtcgtGCTAAAGGCCGCGCTCGTCTGACCGTTTCTAAGTAATAAAGCTAACCCCTGAGTGGTTAAGCTCGCATTTCCCAGGGAGTTACGCTTGTTAACTCCCAGTCAATTCACATTCGTCTTCCAGCAGCCACAACGGGCt  <  2:300109/143‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                              
AGCCATTGAATTTATTCAAGTTTAGGTAGAAATCGCCATGAAACGCACTTTTCAACCGTCTGTACTGAAGCGCAACCGTTCTCACGGCTTCCGTGCTCGTATGGCTACTAAAAATGGTCGTCAGGTTCTGGCACGTCGTCGTGCTAAAGGCCGCGCTCGTCTGACCGTTTCTAAGTAATAAAGCTAACCCCTGAGTGGTTAAGCTCGCATTTCCCAGGGAGTTACGCTTGTTAACTCCCAGTCAATTCACATTCGTCTTCCAGCAGCCACAACGGGCT  >  NC_000913/3884299‑3884576

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