Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,946,631 G→C 9.3% P312A (CCG→GCG) ‡ yijP ← phosphoethanolamine transferase CptA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,946,6310GC 9.3% 165.3 / 8.5 54P312A (CCG→GCG) ‡yijPphosphoethanolamine transferase CptA
Reads supporting (aligned to +/- strand):  ref base G (28/21);  new base C (3/2);  total (31/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.33e-01

GGTTGCGCTGGTTGTTCAGGTACACCGGCGCGTCTGTTTGCTGCGAGAAGGTGGTGAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGCTCTTGTCGCTGGCGGCCAGCGCGTCGAGGTTGGGCGTG  >  NC_002947/2946487‑2946779
                                                                                                                                                |                                                                                                                                                    
ggTTGCGCTGGTTGTTCAGGTACACCGGCGCGTCTGTTTGCTGCGAGAAGGTGGTGAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGc                                                                                                                                                 >  1:517331/1‑150 (MQ=255)
 gTTGCGCTGGTTGTTCAGGTACACCGGCGCGTCTGTTTGCTGCGAGAAGGTGGTGAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCt                                                                                                                                                <  1:1002371/150‑1 (MQ=255)
         ggTTGTTCAGGTACACCGGCGCGTCTGTTTGCTGCGAGAAGGTGGTGAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAAc                                                                                                                                        <  1:944898/150‑1 (MQ=255)
         ggTTGTTCAGGTACACCGGCGCGTCTGTTTGCTGCGAGAAGGTGGTGAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAAc                                                                                                                                        <  2:517331/150‑1 (MQ=255)
              ttCAGGTACACCGGCGCGTCTGTTTGCTGCGAGAAGGTGGTGAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAgg                                                                                                                                            >  2:382557/1‑141 (MQ=255)
              ttCAGGTACACCGGCGCGTCTGTTTGCTGCGAGAAGGTGGTGAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAgg                                                                                                                                            <  1:382557/141‑1 (MQ=255)
                         cGGCGCGTCTGTTTGCTGCGAGAAGGTGGTGAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCtc                                                                                                                        >  1:511395/1‑150 (MQ=255)
                                 cTGTTTGCTGCGAGAAGGTGGTGAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGa                                                                                                                >  1:138866/1‑150 (MQ=255)
                                         tgcGAGAAGGTGGTGAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAgg                                                                                                        >  1:1199412/1‑150 (MQ=255)
                                            gagaAGGTGGTGAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATc                                                                                                     <  2:138866/150‑1 (MQ=255)
                                               aaGGTGGTGAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCtcgtc                                                                                                                     <  2:77356/131‑1 (MQ=255)
                                               aaGGTGGTGAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCtcgtc                                                                                                                     >  1:77356/1‑131 (MQ=255)
                                                 ggtggtGAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATctgctg                                                                                                <  1:53883/150‑1 (MQ=255)
                                                  gtggtgAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCGAGTAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGCTCAGCGCCGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCtgctgc                                                                                               >  2:165863/1‑150 (MQ=255)
                                                   tggtgAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCa                                                                                              >  2:634828/1‑150 (MQ=255)
                                                          cATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAATGTGAGGATCTGCTGCATCACTTc                                                                                       >  2:883115/1‑150 (MQ=255)
                                                              gtgtTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATg                                                                                   >  1:498117/1‑150 (MQ=255)
                                                                   gcgcTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTa                                                                              >  2:1378/1‑150 (MQ=255)
                                                                    cgcTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTAc                                                                             <  1:672457/150‑1 (MQ=255)
                                                                      cTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACgg                                                                           >  2:176640/1‑150 (MQ=255)
                                                                           tCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGa                                                                                     <  1:730879/135‑1 (MQ=255)
                                                                           tCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGa                                                                                     >  2:730879/1‑135 (MQ=255)
                                                                           tCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGcgcg                                                                      >  2:78225/1‑150 (MQ=255)
                                                                                gTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGc                                                                                                                          >  1:778836/1‑93 (MQ=255)
                                                                                gTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGc                                                                                                                          <  2:778836/93‑1 (MQ=255)
                                                                                gTCTGATGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGcgcgacgtg                                                                 >  2:791495/1‑145 (MQ=255)
                                                                                 tCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACCCGTCGCTGa                                                                                                                                              <  1:924542/72‑1 (MQ=255)
                                                                                 tCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACCCGTCGCTGa                                                                                                                                              >  2:924542/1‑72 (MQ=255)
                                                                                   tgctgGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATAcc                                                              >  2:958665/1‑150 (MQ=255)
                                                                                        gttggtGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTc                                                                                       <  1:41421/120‑1 (MQ=255)
                                                                                        gttggtGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTc                                                                                       >  2:41421/1‑120 (MQ=255)
                                                                                                tCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAaca                                                 >  2:157468/1‑150 (MQ=255)
                                                                                                 ccAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAacac                                                >  1:98461/1‑150 (MQ=255)
                                                                                                   aGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACgg                                              >  1:498500/1‑150 (MQ=255)
                                                                                                        gTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAgg                                         >  1:95054/1‑150 (MQ=255)
                                                                                                         tCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAATGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGAGCGGCGATACCACGTCGCGGAACACGGTCAGGc                                        <  1:883115/150‑1 (MQ=255)
                                                                                                         tCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGc                                        >  1:1052722/1‑150 (MQ=255)
                                                                                                           ttGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAgg                                         >  2:1064578/1‑147 (MQ=255)
                                                                                                           ttGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAgg                                         <  1:1064578/147‑1 (MQ=255)
                                                                                                             gTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGTTCtt                                    <  1:958665/150‑1 (MQ=255)
                                                                                                                gcctgcctGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGCTCTTGTc                                 >  1:54743/1‑150 (MQ=255)
                                                                                                                  ctgcctgTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGCTCTTGTCGc                               >  1:296600/1‑150 (MQ=255)
                                                                                                                       tgTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGCTCTTGTCGCTggcg                          <  1:176640/150‑1 (MQ=255)
                                                                                                                       ggTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGGCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCCCATCACACGGTCAGGCTCTTGTCGCTggcg                          <  1:997978/149‑1 (MQ=255)
                                                                                                                               cagGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGCTCTTGt                                  >  2:237386/1‑134 (MQ=255)
                                                                                                                               cagGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGCTCTTGt                                  <  1:237386/134‑1 (MQ=255)
                                                                                                                                 ggTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGTGGAACACGGTCAGGCTCTTGTCGCTGGCGGCCAGCGCGt                <  1:78225/150‑1 (MQ=255)
                                                                                                                                  gTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCAGGTCGCGGAACACGGTCAGGCTCTTTTCGCTGGCGGCCAGCGCGTc               >  1:77676/1‑150 (MQ=255)
                                                                                                                                  gTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGCTCTTGTCGCTGGCGGCCAGCGCGTc               <  2:498117/150‑1 (MQ=255)
                                                                                                                                  gTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGCTCTTGTCGCTGGCGGCCAGCGCGTc               <  2:1052722/150‑1 (MQ=255)
                                                                                                                                      ctCAGCGACCCGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACgg                                                                           >  1:1062116/2‑86 (MQ=255)
                                                                                                                                      ctCAGCGACCCGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACgg                                                                           <  2:1062116/85‑1 (MQ=255)
                                                                                                                                       tCAGCGACCCGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGCTCTTGTCGCTGGCGGCCAGCGCGTCGAGGt          >  2:805995/1‑150 (MQ=255)
                                                                                                                                         aGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGGGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGCTCTTGTCGCTGGCGGCCAGCGCGTCGAGGTTg        <  1:165863/150‑1 (MQ=255)
                                                                                                                                               gggTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGCTCTTGTCGCTGGCGGCCAGCGCGTCGAGGTTGGGCgtg  <  1:634828/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GGTTGCGCTGGTTGTTCAGGTACACCGGCGCGTCTGTTTGCTGCGAGAAGGTGGTGAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGCTCTTGTCGCTGGCGGCCAGCGCGTCGAGGTTGGGCGTG  >  NC_002947/2946487‑2946779

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