Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,154,396 C→A 9.6% G126W (GGG→TGG)  PP_1909 ← septum formation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,154,3960CA 9.6% 144.9 / 6.7 52G126W (GGG→TGG) PP_1909septum formation protein
Reads supporting (aligned to +/- strand):  ref base C (26/21);  new base A (2/3);  total (28/24)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.55e-01

CTGGAACAAGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGCCCAATACCGCAACCTGGTCCGA  >  NC_002947/2154250‑2154534
                                                                                                                                                  |                                                                                                                                          
ctGGAACAAGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCgg                                                                                                                                          >  2:259418/1‑149 (MQ=255)
     aCAAGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCGGGCAATGCCAGGTggat                                                                                                                                     <  3:37455/149‑3 (MQ=255)
       aaGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCt                                                                                                                                   <  7:251243/149‑1 (MQ=255)
         gCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGt                                                                                                                                 <  1:190364/149‑1 (MQ=255)
            cACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCGGGCAATGCCAGGTgga                                                                                                                                      <  6:82829/141‑2 (MQ=255)
            cACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCGGGCAATGCCAGGTgga                                                                                                                                      >  5:82829/1‑140 (MQ=255)
             aCCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTcag                                                                                                                             <  4:121906/149‑1 (MQ=255)
                  cAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTg                                                                                                                        >  7:38360/1‑149 (MQ=255)
                  cAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTACTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTg                                                                                                                        >  8:246110/1‑149 (MQ=255)
                   aGCCCCTCCGAATTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACCTAACGCTCGACCAGCTCCCTGCTCAGTACACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGACATTCCCCGCTGGCGCAGTTCAGCAGTGc                                                                                                                       >  4:41602/1‑149 (MQ=255)
                        cTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGcc                                                                                                                  >  6:92384/1‑149 (MQ=255)
                        cTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGcc                                                                                                                  <  8:38360/149‑1 (MQ=255)
                         tCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGt                                                                                                                         >  6:91483/1‑141 (MQ=255)
                         tCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGt                                                                                                                         <  5:91483/141‑1 (MQ=255)
                         tCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGccc                                                                                                                 >  3:73663/1‑149 (MQ=255)
                            gCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCgg                                                                                                               >  1:104354/1‑148 (MQ=255)
                            gCCTTGAAACTCCCTGCGCAACCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCgg                                                                                                               >  6:230450/1‑148 (MQ=255)
                                 gAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGa                                                                                                          <  5:240237/148‑1 (MQ=255)
                                     cTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCt                                                                                                     >  6:309031/1‑149 (MQ=255)
                                     cTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCt                                                                                                     >  5:59300/1‑149 (MQ=255)
                                         cTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTaa                                                                                                            <  6:65059/138‑1 (MQ=255)
                                         cTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTaa                                                                                                            >  5:65059/1‑138 (MQ=255)
                                          tGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACct                                                                                                >  4:13858/1‑149 (MQ=255)
                                                   cAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTcagca                                                                                                                           <  7:45907/113‑1 (MQ=255)
                                                   cAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTcagca                                                                                                                           >  8:45907/1‑113 (MQ=255)
                                                   cAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGc                                                                                       >  5:119600/1‑149 (MQ=255)
                                                    aGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCt                                                                                      <  4:73663/149‑1 (MQ=255)
                                                     gCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTg                                                                                     <  3:13858/149‑1 (MQ=255)
                                                        gCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCAc                                                                                  >  5:303977/1‑149 (MQ=255)
                                                            ttCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCa                                                                              <  6:303977/149‑1 (MQ=255)
                                                             tCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCa                                                                              <  2:104354/148‑1 (MQ=255)
                                                                    ccACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGc                                                                                                                   <  5:221205/104‑1 (MQ=255)
                                                                    ccACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGc                                                                                                                   >  6:221205/1‑104 (MQ=255)
                                                                    ccACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCagtagtc                                                                        <  7:246110/147‑2 (MQ=255)
                                                                    ccACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGc                                                                      >  3:151358/1‑149 (MQ=255)
                                                                        gTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCg                                                                   >  2:194159/1‑148 (MQ=255)
                                                                        gTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGc                                                                  <  5:92384/149‑1 (MQ=255)
                                                                             gCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAgcg                                                             >  1:234375/1‑149 (MQ=255)
                                                                              cTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCCCTCCAGTAGTTGCTCGCAAgcgc                                                            >  8:237717/1‑149 (MQ=255)
                                                                                       ctccctgctcAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAAc                                                    <  5:230450/148‑1 (MQ=255)
                                                                                              ctcAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGc                                              <  1:194159/147‑1 (MQ=255)
                                                                                               tcAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACt                                                                                 <  5:172733/111‑1 (MQ=255)
                                                                                               tcAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACt                                                                                 >  6:172733/1‑111 (MQ=255)
                                                                                                  gTTCACGCAGTGTCACGGTAAAAGGCACGAAATCCACCTGGCATTGCCCGGTGGCGGTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTg                                        >  4:187808/1‑149 (MQ=255)
                                                                                                        gCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGc                                  <  3:41602/149‑1 (MQ=255)
                                                                                                                        aGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGCCCAATa                  >  3:65091/1‑149 (MQ=255)
                                                                                                                            aCGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGCCCAATACCg               >  3:160853/1‑148 (MQ=255)
                                                                                                                                   ccACCGGGCAATGCCAGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTAt                                                <  1:131863/108‑1 (MQ=255)
                                                                                                                                   ccACCGGGCAATGCCAGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTAt                                                >  2:131863/1‑108 (MQ=255)
                                                                                                                                    cACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGCCCAATACCGCAACCTGGt      >  5:52493/1‑149 (MQ=255)
                                                                                                                                       cTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGCCCAATACCGCAACCTGGTCCg   <  6:113262/149‑1 (MQ=255)
                                                                                                                                          gCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGCCCAATACCGCAACCTGGTCCGa  <  2:234375/147‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                          
CTGGAACAAGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGCCCAATACCGCAACCTGGTCCGA  >  NC_002947/2154250‑2154534

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