Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 665,094 G→C 9.3% P33A (CCT→GCT)  rlpA ← septal ring protein, suppressor of prc, minor lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913665,0940GC 9.3% 141.8 / 6.9 43P33A (CCT→GCT) rlpAseptal ring protein, suppressor of prc, minor lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (19/20);  new base C (2/2);  total (21/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.93e-01

CGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTG  >  NC_000913/664951‑665241
                                                                                                                                               |                                                                                                                                                   
cGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGCCTGCg                                                                                                                                                <  1:325223/149‑1 (MQ=255)
 gCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCg                                                                                                                                                >  3:53367/1‑148 (MQ=255)
 gCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGCCTGCgg                                                                                                                                               <  2:338520/149‑1 (MQ=255)
    tGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTa                                                                                                                                             <  4:175989/148‑1 (MQ=255)
     ggCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTaca                                                                                                                                           <  2:335809/149‑1 (MQ=255)
       cTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACt                                                                                                                                         >  1:275242/1‑149 (MQ=255)
        tAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACt                                                                                                                                         <  4:101435/148‑1 (MQ=255)
               agaCGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCt                                                                                                                                 >  2:370422/1‑149 (MQ=255)
                  cGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTacac                                                                                                                                          >  3:166153/1‑137 (MQ=255)
                  cGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTacac                                                                                                                                          <  4:166153/137‑1 (MQ=255)
                     aTCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCt                                                                                                                           <  1:79239/149‑1 (MQ=255)
                       ccTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTg                                                                                                                          <  4:102098/148‑1 (MQ=255)
                                   tGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCt                                                                                                              <  4:284951/148‑1 (MQ=255)
                                         ttttACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACAt                                                                                                       >  1:36073/1‑149 (MQ=255)
                                                       tGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTcc                                                                                         >  1:229061/1‑149 (MQ=255)
                                                       tGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTc                                                                                          >  4:68906/1‑148 (MQ=255)
                                                        ggTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTcc                                                                                         <  4:238626/148‑1 (MQ=255)
                                                          tAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGc                                                                                      <  1:41599/149‑1 (MQ=255)
                                                                tGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCt                                                                                        <  2:272220/141‑1 (MQ=255)
                                                                tGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCt                                                                                        >  1:272220/1‑141 (MQ=255)
                                                                 gATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGAt                                                                                >  3:85053/1‑148 (MQ=255)
                                                                        cGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATcc                                                                        <  2:129805/149‑1 (MQ=255)
                                                                                tCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACcat                                                                                                                     >  4:265689/1‑96 (MQ=255)
                                                                                tCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACcat                                                                                                                     <  3:265689/96‑1 (MQ=255)
                                                                                       ttCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGa                                                                     >  1:152825/1‑137 (MQ=255)
                                                                                       ttCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGa                                                                     <  2:152825/137‑1 (MQ=255)
                                                                                       ttCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTa                                                          >  4:241225/1‑148 (MQ=255)
                                                                                          gAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCa                                                      >  2:264609/1‑149 (MQ=255)
                                                                                              cgcgGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTg                                                   >  4:81405/1‑148 (MQ=255)
                                                                                                       gCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCtt                                         >  1:110544/1‑149 (MQ=255)
                                                                                                                           aGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACAt                      >  3:141946/1‑148 (MQ=255)
                                                                                                                              gCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACAtttt                   <  3:319557/148‑1 (MQ=255)
                                                                                                                               ccGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTc                  >  3:211261/1‑148 (MQ=255)
                                                                                                                                  ttACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCtg               <  3:128321/148‑1 (MQ=255)
                                                                                                                                  aagtgtACCGCAGCCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATACCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCtg               >  4:330864/6‑148 (MQ=255)
                                                                                                                                   tACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCtgt              <  3:340383/148‑1 (MQ=255)
                                                                                                                                    aCATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCtgtg             >  3:310647/1‑148 (MQ=255)
                                                                                                                                    aCATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTgg            <  2:35298/149‑1 (MQ=255)
                                                                                                                                      aTACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTcc                 >  2:404695/1‑142 (MQ=255)
                                                                                                                                      aTACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTcc                 <  1:404695/142‑1 (MQ=255)
                                                                                                                                       tACCGCAGCCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGgtg         >  1:240928/1‑149 (MQ=255)
                                                                                                                                        aCCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGgtg         <  4:211261/148‑1 (MQ=255)
                                                                                                                                        aCCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGgtg         <  4:53367/148‑1 (MQ=255)
                                                                                                                                               gCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTg  <  3:241225/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
CGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTG  >  NC_000913/664951‑665241

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