Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,741,445 C→T 7.9% V156V (GTG→GTA)  PP_1538 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,741,4450CT 7.9% 141.9 / 2.0 51V156V (GTG→GTA) PP_1538hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (27/19);  new base T (2/2);  total (30/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

AGTAGTGTTCGCCATCGCCGTGGTGAATGCGGTTTATGCCCTCGGGGACGTAGCGGCTTTCGAAATCCGTCTCGCCCTCGCCATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCGGGGCATGTGGATGCTGGGCTTCCATGCTGACCGTGGATCGCCAACGCCATCCTCGTCACCTGCCCATTCGCAGGGTCCGTCTGCTGCGTAGACAATGTGGCCTGAGTAATAGCCGCGCCCGAACCGCATTTCGCGGATAGCCT  >  NC_002947/1741309‑1741593
                                                                                                                                        |                                                                                                                                                    
agtagtGTTCGCCATCGCCGTGGTGAATGCGGTTTATGCCCTCGGGGACGTAGCGGCTTTCGAAATCCGTCTCGCCCTCGCCATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCGGGGCAtg                                                                                                                                          >  6:188833/1‑149 (MQ=255)
  tagtGTTCGCCATCGCCGTGGTGAATGCGGTTTATGCCCTCGGGGACGTAGCGGCTTTTGAAATCCGTCTCGCCCTCGCCATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCGGGGCAtgtg                                                                                                                                        >  1:323412/1‑149 (MQ=255)
     tgtTCGCCATCGCCGTGGTGAATGCGGTTTATGCCCTCGGGGACGTAGCGGCTTTCGAAATCCGTCTCGCCCTCGCCATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCGGGGCATGTGGAt                                                                                                                                     >  8:311351/1‑149 (MQ=255)
       ttCGCCATCGCCGTGGTGAATGCGGTTTATGCCCTCGGGGACGTAGCGGCTTTCGAAATCCGTCTCGCCCTCGCCATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCGGGGCAt                                                                                                                                           <  7:126269/141‑1 (MQ=255)
       ttCGCCATCGCCGTGGTGAATGCGGTTTATGCCCTAGGGGACGTAGCGCCTTTCGAAATCCGTCTCGCCATCGCCATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCCGCCACTCGCGGCAt                                                                                                                                           >  8:126269/1‑141 (MQ=255)
          gCCATCGCCGTGGTGAATGCGGTTTATGCCCTCGGGGACGTAGCGGCTTTCGAAATCCGTCTCGCCCTCGCCATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCGGGGCATGTGGATGCTg                                                                                                                                 <  8:197363/148‑1 (MQ=255)
            cATCGCCGTGGTGAATGCGGTTTATGCCCTCGGGGACGTAGCGGCTTTCGAAATCCGTCTCGCCCTCGCCATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCGGGGCATGTGGATGCTGGGc                                                                                                                              >  3:111202/1‑149 (MQ=255)
             aTCGCCGTGGTGAATGCGGTTTATGCCCTCGGGGACGTAGCGGCTTTCGAAATCCGTCTCGCCCTCGCCATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCGGGGCATGTGGATGCTGGGCt                                                                                                                             <  4:111202/149‑1 (MQ=255)
                    tggtgAATGCGGTTTATGCCCTCGGGGACGTAGCGGCTTTCGAAATCCGTCTCGCCCTCGCCATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCGGGGCATGTGGATGCTGGGCTTCCATGc                                                                                                                      >  1:314565/1‑149 (MQ=255)
                    tggtgAATGCGGTTTATGCCCTCGGGGACGTAGCGGCTTTCGAAATCCGTCTCGCCCTCGCCATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCGGGGCATGTGGATGCTGGGCTTCCATGc                                                                                                                      <  7:257311/149‑1 (MQ=255)
                        gAATGCGGTTTATGCCCTCGGGGACGTAGCGGCTTTCGAAATCCGTCTCGCCCTCGCCATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCGGGGCATGTGGATGCTGGGCTTCCATGCTGAc                                                                                                                  >  5:184710/1‑149 (MQ=255)
                           tGCGGTTTATGCCCTCGGGGACGTAGCGGCTTTCGAAATCCGTCTCGCCCTCGCCATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCGGGGCATGTGGATGCTGGGCTTCCATGCTGACCGt                                                                                                               >  8:139447/1‑149 (MQ=255)
                             cGGTTTATGCCCTCGGGGACGTAGCGGCTTTCGAAATCCGTCTCGCCCTCGCCATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCGGGGCATGTGGATGCTGGGCTTCCATGCTGACCGTg                                                                                                              >  7:110072/1‑148 (MQ=255)
                                tttATGCCCTCGGGGACGTAGCGGCTTTCGAAATCCGTCTCGCCCTCGCCATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCGGGGCATGTGGATGCTGGGCTTCCATGCTGACCGTGGATc                                                                                                          <  5:188833/149‑1 (MQ=255)
                                tttATGCCCTCGGGGACGTAGCGGCTTTCGAAATCCGTCTCGCCCTCGCCATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCGGGGCATGTGGATGCTGGGCTTCCATGCTGACCGTGGATc                                                                                                          >  4:108069/1‑149 (MQ=255)
                                   aTGCCCTCGGGGACGTAGCGGCTTTCGAAATCCGTCTCGCCCTCGCCATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCGGGGCATGTGGATGCTGGGCTTCCAt                                                                                                                        >  5:161045/1‑132 (MQ=255)
                                   aTGCCCTCGGGGACGTAGCGGCTTTCGAAATCCGTCTCGCCCTCGCCATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCGGGGCATGTGGATGCTGGGCTTCCAt                                                                                                                        <  6:161045/132‑1 (MQ=255)
                                         tCGGGGACGTAGCGGCTTTCGAAATCCGTCTCGCCCTCGCCATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCGGGGCATGTGGATGCTGGGCTTCCATGCTGACCGTGGATCGCCAACGc                                                                                                  >  7:41546/1‑148 (MQ=255)
                                         tCGGGGACGTAGCGGCTTTCGAAATCCGTCTCGCCCTCGCCATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCGGGGCATGTGGATGCTGGGCTTCCATGCTGACCGTGGATCGCCAACGc                                                                                                  >  7:47230/1‑148 (MQ=255)
                                               aCGTAGCGGCTTTCGAAATCCGTCTCGCCCTCGCCATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCGGGGCATGTGGATGCTGGGCTTCCATGCTGACCGTGGATCGCCAACGCCATCCTc                                                                                           <  2:314565/149‑1 (MQ=255)
                                                    gCGGCTTTCGAAATCCGTCTCGCCCTCGCCATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCggg                                                                                                                                               <  2:352073/92‑1 (MQ=255)
                                                    gCGGCTTTCGAAATCCGTCTCGCCCTCGCCATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCggg                                                                                                                                               >  1:352073/1‑92 (MQ=255)
                                                    gCGGCTTTCGAAATCCGTCTCGCCCTCGCCATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCGGGGCATGTGGATGCTGGGCTTCCATGCTGACCGTGGATCGCCAACGCCATCCTCGTCAc                                                                                      <  7:139447/149‑1 (MQ=255)
                                                      ggCTTTCGAAATCCGTCTCGCCCTCGCCATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCGGGGCATGTGGATGCTGGGCTTCCATGCTGACCGTGGATCGCCAACGCCATCCTCGTCACCt                                                                                    <  6:184710/149‑1 (MQ=255)
                                                          ttCGAAATCCGTCTCGCCCTCGCCATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCAACGCGGGGCATGTGGATGCTGGGCTTCCATGCTGACCGTGGATCGCCAACGCCATCCTCGTCACCTggc                                                                                 >  5:164128/1‑146 (MQ=255)
                                                            cGAAATCCGTCTCGCCCTCGCCATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCCGCAGGATTCGGCCGGCCCCGCGGCGCATGTGGATGCTGGGCTTCCATGCTGACCGTGGATCGCCAACGCCCTCCTCGTCACCTGCCCAt                                                                              >  2:245783/1‑149 (MQ=255)
                                                              aaaTCCGTCTCGCCCTCGCCATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCGGGGCATGTGGATGCTGGGTTTCCATGCTGACCGTGGATCGCCAAc                                                                                                    <  6:17386/125‑1 (MQ=255)
                                                              aaaTCCGTCTCGCCCTCGCCATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCGGGGCATGTGGATGCTGGGTTTCCATGCTGACCGTGGATCGCCAAc                                                                                                    >  5:17386/1‑125 (MQ=255)
                                                                      ctcgccctcgccATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCGGGGCATGTGGATGCTGGGCTTCCATGCTGACCGTGGATCGCCAACGCCATCCTCGTCACCTGCCCATTCGCAGGGTc                                                                    >  5:133367/1‑149 (MQ=255)
                                                                       tcgccctcgccATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCGGGGCATGTGGATGCTGGGCTTCCATGCTGACCGTGGATCGCCAACGCCATCCTCGTCACCTGCCCATTCGCAGGGTcc                                                                   >  8:274461/1‑149 (MQ=255)
                                                                        cgccctcgccATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCTACGCGGGGCATGTGGATGCTGGGCTTCCATGCTGACCGTGGATCGCCAACGCCATCCTCGTCACCTGCCCATTCGCAGGGTCCg                                                                  >  2:324926/1‑149 (MQ=255)
                                                                        cgccctcgccATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCTACGCGGGGCATGTGGATGCTGGGCTTCCATGCTGACCGTGGATCGCCAACGCCATCCTCGTCACCTGCCCATTCGCAGGGTCCg                                                                  >  2:320757/1‑149 (MQ=255)
                                                                          ccctcgccATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCGGGGCATGTGGATGCTGGGCTTCCATGCTGACCGTGGATCGCCAACGCCATCCTCGTCACCTGc                                                                                  <  4:235909/131‑1 (MQ=255)
                                                                          ccctcgccATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCGGGGCATGTGGATGCTGGGCTTCCATGCTGACCGTGGATCGCCAACGCCATCCTCGTCACCTGc                                                                                  >  3:235909/1‑131 (MQ=255)
                                                                               gccATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCGGGGCATGTGGATGCTGGGCTTCCATGCTGACCGTGGATCGCCAACGCCATCCTCGTCACCTGCCCATTCGCAGGGTCCGTctgctg                                                           >  2:279833/1‑149 (MQ=255)
                                                                                 cATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCGGGGCATGTGGATGCTGGGCTTCCATGCTGACCGTGGATCGCCAACGCCATCCTCGTCACCTGCCCATTCGCAGGGTCCGTCTGCTGCg                                                         >  2:188174/1‑149 (MQ=255)
                                                                                  aTCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCTACGCGGGGCATGTGGATGCTGGGCTTCCATGCTGACCGTGGATCGCCAACGCCATCCTCGTCACCTGCCCATTCGCAGGGTCCGTCTGCTGCGt                                                        <  1:324926/149‑1 (MQ=255)
                                                                                  aTCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCTACGCGGGGCATGTGGATGCTGGGCTTCCATGCTGACCGTGGATCGCCAACGCCATCCTCGTCACCTGCCCATTCGCAGGGTCCGTCTGCTGCGt                                                        <  1:320757/149‑1 (MQ=255)
                                                                                  aTCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCGGGGCATGTGGATGCTGGGCTTCCATGCTGACCGTGGATCGCCAACGCCATCCTCGTCACCTGCCCATTCGCAGGGTCCGTCTGCTGCGt                                                        <  4:175552/149‑1 (MQ=255)
                                                                                        cAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCGGGGCATGTGGATGCTGGGCTTCCATGCTGACCGTGGATCGCCAACGCCATCCTCGTCACCTGCCCATTCGCAGGGTCCGTCTGCTGCGTAGAc                                                    <  1:279833/147‑1 (MQ=255)
                                                                                             gcTCGACTCGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCGGGGCATGTGGATGCTGGGCTTCCATGCTGACCGTGGATCGCCAACGCCATCCTCGTCACCTGCCCATTCGCAGGGTCCGTCTGCTGCGTAGACAATGTgg                                             <  8:41546/149‑1 (MQ=255)
                                                                                             gcTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCGGGGCATGTGGATGCTGGGCTTCCATGCTGACCGTGGATCGCCAACGCCATCCTCGTCACCTGCCCATTCGCAGGGTCCGTCTGCTGCGTAGACAATGTgg                                             <  8:47230/149‑1 (MQ=255)
                                                                                                   cacgcacgTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCGGGGCATGTGGATGCTGGGCTTCCATGCTGACCGTGGATCGCCAACGCCATCCTCGTCACCTGCCCATTCGCAGGGTCCGTCTGCTGCGTAGACAATGTGGCCTGAg                                       <  7:311351/149‑1 (MQ=255)
                                                                                                      gcacgTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCGGGGCATGTGGATGCTGGGCTTCCATGCTGACCGTGGATCGCCAACGGCAGCCTGGTCCCCTGCCGATTCGCAGGGTCCGTCTGCTGCGTTGACATGGGGGCGTGGGTaa                                    >  4:29753/1‑149 (MQ=255)
                                                                                                      gcacgTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCGGGGCATGTGGATGCTGGGCTTCCATGCTGACCGTGGATCGCCAACGCCATcgtcgt                                                                                         >  6:45563/1‑96 (MQ=255)
                                                                                                      gcacgTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCGGGGCATGTGGATGCTGGGCTTCCATGCTGACCGTGGATCGCCAACGCCATCCTCGTCACCTGCCCATTCGCAGGGTCCGTCTGCTGCGTAGACAAt                                                 <  2:251006/136‑1 (MQ=255)
                                                                                                      gcacgTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCGGGGCATGTGGATGCTGGGCTTCCATGCTGACCGTGGATCGCCAACGCCATCCTCGTCACCTGCCCATTCGCAGGGTCCGTCTGCTGCGTAGACAAt                                                 >  1:251006/1‑136 (MQ=255)
                                                                                                                    tCCAGCAGGATTCGGCCGGCCACGCGGGGCATGTGGATGCTGGGCTTCCATGCTGACCGTGGATCGCCAACGCCATCCTCGTCACCAGCCCATTCGCAGGGTCCGTCTGCTGCGTAGACAATGTGGCCTGAGTAATAGCCGCGCCCGaa                      >  2:7037/1‑149 (MQ=255)
                                                                                                                       agcagGATTCGGCCGGCCACGCGGGGCATGTGGATGCTGGGCTTCCATGCTGACCGTGGATCGCCAACGCCATCCTCGTCACCTGCCCATTCGCAGGGTCCGTCTGCTGCGTAGACAATGTGGCCTGAGTAATAGCCGCGCCCGAACCg                   >  5:267965/1‑149 (MQ=255)
                                                                                                                              ttCGGCCGGCCACGCGGGGCATGTGGATGCTGGGCTTCCATGCTGACCGTGGATCGCCAACGCCATCCTCGTCACCTGCCCATTCGCAGGGTCCGTCTGCTGCGTAGACAATGTGGCCTGAGTAATAGCCGCGCCCGAACCGCATTTcg            >  3:268558/1‑149 (MQ=255)
                                                                                                                              ttCGGCCGGCCACGCGGGGCATGTGGATGCTGGGCTTCCATGCTGACCGTGGATCGCCAACGCCATCCTCGTCACCAGCCCATTCGCAGGGTCCGTCTGCTGCGTAGACAATGTGGCCTGAGTAATAGCCGCGCCCGAACCGCATTTcg            <  1:7037/149‑1 (MQ=255)
                                                                                                                                    cggccACGCGGGGCATGTGGATGCTGGGCTTCCATGCTGACCGTGGATCGCCAACGCCATCCTCGTCACCTGCCCATTCGCAGGGTCCGTCTGCTGCGTAGACAATGTGGCCTGAGTAATAGCCGCGCCCGAACCGCATTTCGCGGa        >  5:243597/1‑147 (MQ=255)
                                                                                                                                      gccACGCGGGGCATGTGGATGCTGGGCTTCCATGCTGACCGTGGATCGCCAACGCCATCCTCGTCACCTGCCCATTCGCAGGGTCCGTCTGCTGCGTAGACAATGTGGCCTGAGTAATAGCCGCGCCCGAACCGCATTTCGCGGATAGc    >  5:220850/1‑149 (MQ=255)
                                                                                                                                       ccACGCGGGGCATGTGGATGCTGGGCTTCCATGCTGACCGTGGATCGCCAACGCCATCCTCGTCACCTGCCCATTCGCAGGGTCCGTCTGCTGCGTAGACAATGTGGCCTGAGTAATAGCCGCGCCCGAACCGCATTTCGCGGATAGcc   >  3:246372/1‑149 (MQ=255)
                                                                                                                                        cACGCGGGGCATGTGGATGCTGGGCTTCCATGCTGACCGTGGATCGCCAACGCCATCCTCGTCACCTGCCCATTCGCAGGGTCCGTCTGCTGCGTAGACAATGTGGCCTGAGTAATAGCCGCGCCCGAACCGCATTTCGCGGATAGCCt  <  7:274461/149‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                                    
AGTAGTGTTCGCCATCGCCGTGGTGAATGCGGTTTATGCCCTCGGGGACGTAGCGGCTTTCGAAATCCGTCTCGCCCTCGCCATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCGGGGCATGTGGATGCTGGGCTTCCATGCTGACCGTGGATCGCCAACGCCATCCTCGTCACCTGCCCATTCGCAGGGTCCGTCTGCTGCGTAGACAATGTGGCCTGAGTAATAGCCGCGCCCGAACCGCATTTCGCGGATAGCCT  >  NC_002947/1741309‑1741593

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: ‑