Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,332,822 G→A 11.2% intergenic (‑87/+40) yhbE ← / ← rpmA inner membrane protein YhbE/50S ribosomal subunit protein L27

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,332,8220GA11.2% 101.2 / 6.5 45intergenic (‑87/+40)yhbE/rpmAinner membrane protein YhbE/50S ribosomal subunit protein L27
Reads supporting (aligned to +/- strand):  ref base G (23/17);  new base A (2/3);  total (25/20)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.01e-01

GCTGTGGTGAGCGCCAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCACTTTACCGTCTGCTTTAGCAAACAGAGTGTGGTCACG  >  NC_000913/3332689‑3332957
                                                                                                                                     |                                                                                                                                       
gCTGTGGTGAGCGCCAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACACGTTGCggg                                                                                                                                <  1:729027/143‑1 (MQ=255)
gCTGGGGTGAGCGCCAAAAGAATTCTAAGGCCTGGCTGCGGCTTCATACCGTATTCCCTGCAAAGCGCTTTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACACGTTGCggg                                                                                                                                <  2:165616/143‑1 (MQ=255)
 cTGTGGTGAGCGCCAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACACGTTGCgggg                                                                                                                               <  2:344247/143‑1 (MQ=255)
       tgAGCGCCAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCa                                                                                                                         >  2:769037/1‑143 (MQ=255)
        gAGCGCCAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCAt                                                                                                                        <  2:845880/143‑1 (MQ=255)
               aaaaGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGGGTTGCGGGGCTTTCATCCGTTAc                                                                                                                 >  1:506706/1‑143 (MQ=255)
                     atgccaatgccTGCCTGCTGCTTCATACCGTATTCCCTGCATAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGAc                                                                                                           <  2:226456/143‑1 (MQ=255)
                      tgccaatgccTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCCTTTTTACCGGTTATCGAACGTAAAAAGCCCCGCAAAGGTTCGCCGGGCGTTCATCCGTTACCGGGAcg                                                                                                          >  2:8723/1‑143 (MQ=255)
                         caatgccTGCCTGCTGTTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGa                                                                                                       >  2:423977/1‑143 (MQ=255)
                               ctgcctgcTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACt                                                                                                 >  2:55863/1‑143 (MQ=255)
                               ctgcctgcTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCAGTTATCGCATGTAAAAACCCCCGCAACGTGTTGCGGGGCTTGCATCCGCTACCTGGAcgcgcagacct                                                                                                 >  2:495100/1‑136 (MQ=255)
                                tgcctgcTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACtt                                                                                                >  2:467975/1‑143 (MQ=255)
                                   ctgctgCTTCATACCGTATTCCCTGCAAAGCGCATTCGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTAtt                                                                                             <  1:41515/143‑1 (MQ=255)
                                       tgctTCATACCGTGTTCCCTGCAAAGCGCAGTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGc                                                                                         >  1:897669/1‑143 (MQ=255)
                                       tgctTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGTTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATACTGTACCTGGACGCGAGAAAATTCTTCAGc                                                                                         >  2:533465/1‑143 (MQ=255)
                                         ctTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTAACCGCGACGCGCACAACTTACACAGCtt                                                                                       >  2:319535/1‑143 (MQ=255)
                                               aCCGTATTCCCTGCAATGCGCATTTGCGGCTTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGc                                                                                 <  2:280297/143‑1 (MQ=255)
                                                     gcctcTGCAAAGCGCATTTGCGGACTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATaa                                                                           <  1:292940/139‑1 (MQ=255)
                                                         cTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAAttt                                                                       >  2:467616/1‑143 (MQ=255)
                                                         cTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGCGCTGCGGGGCGTTCATCCGTTCCCGCGGCGCCACAAACTTACTCAGCTACGCTGCTGATgaaatc                                                                       >  2:73144/1‑140 (MQ=255)
                                                           ttaaaGCTCAATTGCGGAGTAACTTCGCACCATGCGCGTTTTTACTGGTTATCTAATGTAAAAAGCGCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTAc                                                                     <  2:118372/141‑1 (MQ=255)
                                                           gCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTAc                                                                     <  2:208551/143‑1 (MQ=255)
                                                                    aTTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCg                                                            >  1:378149/1‑143 (MQ=255)
                                                                     tttGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCgg                                                           >  1:800371/1‑143 (MQ=255)
                                                                             ttAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTACATGTGGATAAATTTACGGTTTTTCGGGCCTTTaa                                                   <  1:474383/143‑1 (MQ=255)
                                                                             gtgaCGCCGCGCCATGCGTATTGGTACCGGTTAGCGAATGGACAAAGCGCCGCAACGTGTTGCGTGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTaa                                                   <  2:343961/140‑1 (MQ=255)
                                                                                  gTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCg                                              >  1:748370/1‑143 (MQ=255)
                                                                                    cGCACCGTGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGaa                                            <  1:274892/143‑1 (MQ=255)
                                                                                    cGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGaa                                            <  1:498586/143‑1 (MQ=255)
                                                                                     gCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAAt                                           >  1:312639/1‑143 (MQ=255)
                                                                                       aCCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGCTCAATTTCCGGTTTTTCGTGCTTTTAACTTCGAAttt                                         >  1:363395/1‑143 (MQ=255)
                                                                                          aTGCGCATTTTTACCGGTTATCGAATGCAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCAc                                      <  2:202996/143‑1 (MQ=255)
                                                                                          aTGCGCATATCAACCGGATAACCAACGTACTAACCCCCGCAACGTTTGGCGGCGCTATCATCCGATACCGGGACGCGAAAAACGTATTCAGCTTCGATGCTGATAAATGTACGGGATTGCGGGCCTTCAACTTCGAATTacac                                      >  2:199524/1‑143 (MQ=255)
                                                                                              gcATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTGCGGTTCTTCGGGCCTTTATCTTTGACTTTCCCTTTa                                  >  1:30117/1‑143 (MQ=255)
                                                                                                aTTTTGACCGGTGATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCACTTTAcc                                <  1:898172/143‑1 (MQ=255)
                                                                                                      aCCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCACTTTACCGTCTGc                          <  2:638627/143‑1 (MQ=255)
                                                                                                          gTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCACTTTACCGTCTGCTTTa                      >  2:230923/1‑143 (MQ=255)
                                                                                                               cGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCACTTTACCGTCTGCTTTAGCaaa                 <  1:423977/143‑1 (MQ=255)
                                                                                                               cGAACGTAAAAAGCCCCGCCACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCACTTTACCGTCTGCTTTAGCaaa                 <  1:475548/143‑1 (MQ=255)
                                                                                                                            ccccGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGCACAACTTATTCAGCTTCGATGCTGATAAATTTACGGTGTTTCGGGCCTTTAACTTCGAATTTCACTTTCCCGTCTGCTTTGGCAAACAGAGTGTGGTCa    >  1:196859/1‑143 (MQ=255)
                                                                                                                            ccccGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCACTTTACCGTCTGCTTTAGCAAACAGAGTGTGGTCa    >  1:195121/1‑143 (MQ=255)
                                                                                                                            ccccGCAACGTGTTGCGGGGCTCTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGTCCTTTAACTTCGAATTTCACTTTACCGTCTGCTTTAGCAAACAGAGTGTGGTCa    <  1:131988/143‑1 (MQ=255)
                                                                                                                             cccGCAACACGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCACTTTACCGTCTGCTTTAGCAAACAGAGGGTGGTCAc   >  1:514737/1‑143 (MQ=255)
                                                                                                                             cccGCAACACGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCACTTTACCGTCTGCTTTAGCAAACAGAATGTGGTCAc   >  2:471665/1‑143 (MQ=255)
                                                                                                                              ccGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGACAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCGTTTAACTTCGCATTGCACTTTACCGCCTGCTTTAGCAAACAGCGAGTGCTCACg  >  2:112221/1‑143 (MQ=255)
                                                                                                                                     |                                                                                                                                       
GCTGTGGTGAGCGCCAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCACTTTACCGTCTGCTTTAGCAAACAGAGTGTGGTCACG  >  NC_000913/3332689‑3332957

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