Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,707,539 C→T 8.5% D12N (GAC→AAC)  PP_2370 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,707,5390CT 8.5% 129.2 / 2.5 47D12N (GAC→AAC) PP_2370hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (23/20);  new base T (2/2);  total (25/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 = 1.23e-01

GCTGTCCAGCCATACACCCAGGCCGCCATCGCGGGCCTTGAGCACGTGGGCATGGTCGCGCCACATTTGCAGCGGCACGATGTAGTCGTCGCAGTTGGTCAGCTCGTCGAACGAGGTTTCCTTGGGCAGCAGGTGCCAGGTTTCGTCGACGATCTGGTTGTTCTTAATGATTCGCTGCATAGACGCGTTCCTTGAAGGGGTCGATGCCGATACGCTGGTAGGTGTCGATGAAACGCTCTTCCTCGGTACGTTGCTCGACGTACACGGCGATCAGCTTCTCGATCACATCGG  >  NC_002947/2707393‑2707683
                                                                                                                                                  |                                                                                                                                                
gCTGTCCAGCCATACACCCAGGCCGCCATCGCGGGCCTTGAGCACGTGGGCATGGTCGCGCCACATTTGCAGCGGCACGATGTAGTCGTCGCAGTTGGTCAGCTCGTCGAACGAGGTTTCCTTGGGCAGCAGGTGCCAGGTTTCGTcga                                                                                                                                                <  3:245747/149‑1 (MQ=255)
cctctCCAGAAATACAACCAGGCCGCCATCGCGGGCCTTGAGAACGAGGGCATGGTCGCGCCACATTTCCAGCGGCACGATGTAGTCGTCTCAGTTGCACAGCTCGTCGAACGAGGTTTCCTTGGGCAGCAGGAGCAAGGATTCGTcga                                                                                                                                                >  2:177323/5‑149 (MQ=255)
      cAGCCATACACCCAGGCCGCCATCGCGGGCCTTGAGCACGTGGGCATGGTCGCGCCACATTTGCAGCGGCACGATGTAGTCGTCGCAGTTGGTCAGCTCGTCGAACGAGGTTTCCTTGGGCAGCAGGTGCCAGGTTTCGTCGACGATCt                                                                                                                                          >  2:295158/1‑149 (MQ=255)
         ccATACACCCAGGCCGCCATCGCGGGCCTTGAGCACGTGGGCATGGTCGCGCCACATTTGCAGCGGCACGATGTAGTCGTCGCAGTTGGTCAGCTCGTCGAACGAGGTTTCCTTGGGCAGCAGGTGCCAGGTTTCGTTGACGATCTGgt                                                                                                                                       >  1:13288/1‑149 (MQ=255)
         ccATACACCCAGGCCGCCATCGCGGGCCTTGAGCACGTGGGCATGGTCGCGCCACATTTGCAGCGGCACGATGTAGTCGTCGCAGTTGGTCAGCTCGTCGAACGAGGTTTCCTTGGGCAGCAGGTGCCAGGTTTCGTTGACGATCTGgt                                                                                                                                       >  1:8083/1‑149 (MQ=255)
           aTACACCCAGGCCGCCATCGCGGGCCTTGAGCACGTGGGCATGGTCGCGCCACATTTGCAGCGGCACGATGTAGTCGTCGCAGTTGGTCAGCTCGTCGAACGAGGTTTCCTTGGGCAGCAGGTGCCAGGTTTCGTCGACGATCTGgttg                                                                                                                                     >  6:83116/1‑149 (MQ=255)
             acacCCAGGCCGCCATCGCGGGCCTTGAGCACGTGGGCATGGTCGCGCCACATTTGCAGCGGCACGATGTAGTCGTCGCAGTTGGTCAGCTCGTCGAACGAGGTTTCCTTGGGCAGCAGGTGCCAGGTTTCGTCGACGATCTGgttgtt                                                                                                                                   >  1:68717/1‑149 (MQ=255)
              cacCCAGGCCGCCATCGCGGGCCTTGAGCACGTGGGCTTGGTCGCGCCACATTTGCAGCGGCACGATGTAGTCGTCGCAGTTGGTCAGCTCGTCGAACGAGGTTTCCTTGGGCAGCAGGTGCCAGGTTTCGTCGACGATCTGGTTGTTc                                                                                                                                  <  1:177323/149‑1 (MQ=255)
              cacCCAGGCCGCCATCGCGGGCCTTGAGCACGTGGGCATGGTCGCGCCACATTTGCAGCGGCACGATGTAGTCGTCGCAGTTGGTCAGCTCGTCGAACGAGGTTTCCTTGGGCAGCAGGTGCCAGGTTTCGTCGACGATCTGGTTGTTc                                                                                                                                  <  2:68717/149‑1 (MQ=255)
               acCCAGGCCGCCATCGCGGGCCTTGAGCACGTGGGCATGGTCGCGCCACATTTGCAGCGGCACGATGTAGTCGTCGCAGTTGGTCAGCTCGTCGAACGAGGTTTCCTTGGGCAGCAGGTGCCAGGTTTCGTCGACGATc                                                                                                                                           >  1:211357/1‑139 (MQ=255)
               acCCAGGCCGCCATCGCGGGCCTTGAGCACGTGGGCATGGTCGCGCCACATTTGCAGCGGCACGATGTAGTCGTCGCAGTTGGTCAGCTCGTCGAACGAGGTTTCCTTGGGCAGCAGGTGCCAGGTTTCGTCGACGATc                                                                                                                                           <  2:211357/139‑1 (MQ=255)
                 ccAGGCCGCCATCGCGGGCCTTGAGCACGTGGGCATGGTCGCGCCACATTTGCAGCGGCACGATGTAGTGGTCGCAGTTGGTCAGCTCGTCGAACGAGGTTTCCTTGGGCAGCAGGTGCCAGGTTTCGTCGACGATc                                                                                                                                           <  4:166089/137‑1 (MQ=255)
                 ccAGGCCGCCATCGCGGGCCTTGAGCACGTGGGCATGGTCGCGCCACATTTGCAGCGGCACGATGTAGTGGTCGCAGTTGGTCAGCTCGTCGAACGAGGTTTCCTTGGGCAGCAGGTGCCAGGTTTCGTCGACGATc                                                                                                                                           >  3:166089/1‑137 (MQ=255)
                 ccAGGCCGCCATCGCGGGCCTTGAGCACGTGGGCATGGTCGCGCCACATTTGCAGCGGCACGATGTAGTCGTCGCAGTTGGTCAGCTCGTCGAACGAGGTTTCCTTGGGCAGCAGGTGCCAGGTTTCGTcga                                                                                                                                                >  4:13287/1‑132 (MQ=255)
                 ccAGGCCGCCATCGCGGGCCTTGAGCACGTGGGCATGGTCGCGCCACATTTGCAGCGGCACGATGTAGTCGTCGCAGTTGGTCAGCTCGTCGAACGAGGTTTCCTTGGGCAGCAGGTGCCAGGTTTCGTcga                                                                                                                                                <  3:13287/132‑1 (MQ=255)
                 ccAGGCCGCCATCGCGGGCCTTGAGCACGTGGGCATGGTCGCGCCACATTTGCAGCGGCACGATGTAGTCGTCGCAGTTGGTCAGCTCGTCGAACGAGGTTTCCTTGGGCAGCAGGTGCCAGGTTTCGTCGACGATCTGGTTGTTCtt                                                                                                                                >  4:185841/1‑148 (MQ=255)
                 ccAGGCCGCCATCGCGGGCCTTGAGCACGTGGGCATGGTCGCGCCACATTTGCAGCGGCACGATGTAGTCGTCGCAGTTGGTCAGCTCGTCGAACGAGGTTTCCTTGGGCAGCAGGTGCCAGGTTTCGTCGACGATCTGGTTGTTCtt                                                                                                                                >  4:36439/1‑148 (MQ=255)
                     gccgccATCGCGGGCCTTGAGCACGTGGGCATGGTCGCGCCACATTTGCAGCGGCACGATGTAGTCGTCGCAGTTGGTCAGCTCGTCGAACGAGGTTTCCTTGGGCAGCAGGTGCCAGGTTTCGTCGACGATCTGGTTGTTCTTaa                                                                                                                              >  6:79142/1‑146 (MQ=255)
                     gccgccATCGCGGGCCTTGAGCACGTGGGCATGGTCGCGCCACATTTGCAGCGGCACGATGTAGTCGTCGCAGTTGGTCAGCTCGTCGAACGAGGTTTCCTTGGGCAGCAGGTGCCAGGTTTCGTCGACGATCTGGTTGTTCTTaa                                                                                                                              <  5:79142/146‑1 (MQ=255)
                        gccATCGCGGGCCTTGAGCACGTGGGCATGGTCGCGCCACATTTGCAGCGGCACGATGTAGTCGTCGCAGTTGGTCAGCTCGTCGAACGAGGTTTCCTTGGGCAGCAGGTGCCAGGTTTCGTCGACGATCTGGTTGTTCTTAATGATTc                                                                                                                        >  7:87913/1‑149 (MQ=255)
                           aTCGCGGGCCTTGAGCACGTGGGCATGGTCGCGCCACATTTGCAGCGGCACGATGTAGTCGTCGCAGTTGGTCAGCTCGTCGAACGAGGTTTCCTTGGGCAGCAGGTGCCAGGTTTCGTCGACGATCTGGTTGTTCTTAATGATTCGCt                                                                                                                     <  3:185841/149‑1 (MQ=255)
                               cgGGCCTTGAGCACGTGGGCATGGTCGCGCCACATTTGCAGCGGCACGATGTAGTCGTCGCAGTTGGTCAGCTCGTCGAACGAGGTTTCCTTGGGCAGCAGGTGCCAGGTTTCGTCGACGATCTGGTTGTTCTTAATGATTCGCTGCAt                                                                                                                 <  5:83116/149‑1 (MQ=255)
                                   ccTTGAGCACGTGGGCATGGTCGCGCCACATTTGCAGCGGCACGATGTAGTCGTCGCAGTTGGTCAGCTCGTCGAACGAGGTTTCCTTGGGCAGCATGTGCCAGGTTTCGTCGACGATAGTGTTGTTCTTAATGATTCGCTGCATAgcc                                                                                                             >  8:201092/1‑147 (MQ=255)
                                         ggaCGTGGGCATGGTCGCGCCAAATTTGCAGCGGCACGATGTAGTCGTCGCAGTTGGTCAGCTCGTCGAACGAGGTTTCCTTGGGCAGCAGGTGCCAGGTTTCGTCGACGATCTGGTTGTTCTTAATGATTCGCTGCATAGACGCGtt                                                                                                        <  5:263672/146‑1 (MQ=255)
                                         gCACGTGGGCATGGTCGCGCCACATTTGCAGCGGCACGATGTAGTCGTCGCAGTTGGTCAGCTCGTCGAACGAGGTTTCCTTGGGCAGCAGGTGCCAGGTTTCGTTGACGATCTGGTTGTTCTTAATGATTCGCTGCATAGACGCGTTc                                                                                                       <  2:8083/149‑1 (MQ=255)
                                         gCACGTGGGCATGGTCGCGCCACATTTGCAGCGGCACGATGTAGTCGTCGCAGTTGGTCAGCTCGTCGAACGAGGTTTCCTTGGGCAGCAGGTGCCAGGTTTCGTTGACGATCTGGTTGTTCTTAATGATTCGCTGCATAGACGCGTTc                                                                                                       <  2:13288/149‑1 (MQ=255)
                                         gCACGTGGGCATGGTCGCGCCACATTTGCAGCGGCACGATGTAGTCGTCGCAGTTGGTCAGCTCGTCGAACGAGGTTTCCTTGGGCAGCAGGTGCCAGGTTTCGTCGACGATCTGGTTGTTCTTAATGATTCGCTGCATAGACGCGTTc                                                                                                       <  3:74743/149‑1 (MQ=255)
                                                 gCATGGTCGCGCCACATTTGCAGCGGCACGATGTAGTCGTCGCAGTTGGTCAGCTCGTCGAACGAGGTTTCCTTGGGCAGCAGGTGCCAGGTTTCGTCGACGATCTGGTTGTTCTTAATGATTCGCTGCATAGACGCGTTCCTTGAAgg                                                                                               >  1:132894/1‑149 (MQ=255)
                                                  cATGGTCGCGCCTCTTTTTCAGCTGGACGATGTTGTCGTCGCAGTTGGTCAGCGCGTCGATCGAGTTTTCGTTGGGCAGCTGGTGCCAGGTTTCGTCGACGATGTGGTTGTTGTTAATGATTCGCTGCATAGACGCGTTCCTTGAAggg                                                                                              <  8:95674/149‑1 (MQ=255)
                                                               caTTTGCAGCGGCACGATGTAGTCGTCGCAGTTGGTCAGCTCGTCGAACGAGGTTTCCTTGGGCAGCAGGTGCCAGGTTTCGTCGACGATCTGGTTGTTCTTAATGATTCGCTGCATAGACGCGTTCCTTGAAGGGGTCGATGCCGATa                                                                                 <  3:36439/149‑1 (MQ=255)
                                                                      aGCGGCACGATGTAGTCGTCGCAGTTGGTCAGCTCGTCGAACGAGGTTTCCTTGGGCAGCAGGTGCCAGGTTTCGTCGACGATCTGGTTGTTCTTAATGATTCGCTGCATAGACGCGTTCCTTGAAGGGGTCGATGCCGATACGCTgg                                                                           >  8:37690/1‑148 (MQ=255)
                                                                          gCACGATGTAGTCGTCGCAGTTGGTCCGCTCGTCGAACGAGGTTTCCTTGGGCAGCAGGTGCCAGGTTTCGTCGACGATCTGGTTGTTCTTAATGATTCGCTGCATAGACGCGTTCCTTGAAGGGGTCGATGCCGATACGCTGGTAGgt                                                                      <  7:201092/149‑1 (MQ=255)
                                                                           cACGATGTAGTCGTCGCAGTTGGTCAGCTCGTCGAACGAGGTTTCCTTGGGCAGCAGGTGCCAGGTTTCGTCGACGATCTGGTTGTTCTTAATGATTCGCTGCATAGACGCGTTCCTTGAAGGGGTCGATGCCGATACGCTGGTAGgtg                                                                     >  5:281074/1‑149 (MQ=255)
                                                                            aCGATGTAGTCGTCGCAGTTGGTCAGCTCGTCGAACGAGGTTTCCTTGGGCAGCAGGTGCCAGGTTTCGTCGACGATCTGGTTGTTCTTAATGATTCGCTGCATAGACGCGTTCCTTGAAggg                                                                                              <  3:208752/123‑1 (MQ=255)
                                                                            aCGATGTAGTCGTCGCAGTTGGTCAGCTCGTCGAACGAGGTTTCCTTGGGCAGCAGGTGCCAGGTTTCGTCGACGATCTGGTTGTTCTTAATGATTCGCTGCATAGACGCGTTCCTTGAAggg                                                                                              >  4:208752/1‑123 (MQ=255)
                                                                               aTGTAGTCGTCGCAGTTGGTCAGCTCGTCGAACGAGGTTTCCTTGGGCAGCAGGTGCCAGGTTTCGTCGACGATCTGGTTGTTCTTAATGATTCGCTGCATAGACGCGTTCCTTGAAGGGGTCGATGCCGATACGCTGGTAGGTGTCGa                                                                 >  3:206717/1‑149 (MQ=255)
                                                                                tGTAGTCGTCGCAGTTGGTCAGCTCGTCGAACGAGGTTTCCTTGGGCAGCAGGTGCCAGGTTTCGTCGACGATCTGGTTGTTCTTAATGATTCGCTGCATAGACGCGTTCCTTGAAggg                                                                                              >  7:288897/1‑119 (MQ=255)
                                                                                    gtcgtcGCAGTTGGTCAGCTCGTCGAACGAGGTTTCCTTGGGCAGCAGGTGCCAGGTTTCGTCGACGATCTGGTTGTTCTTAATGATTCGCTGCATAGACGCGTTCCTTGAAGGGGTCGATGCCGATACGCTGGTAGGTGTCGATGaaa                                                            >  7:113748/1‑149 (MQ=255)
                                                                                          gCAGTTGGTCAGCTCGTCGAACGAGGTTTCCTTGGGCAGCAGGTGCCAGGTTTCGTCGACGATCTGGTTGTTCTTAATGATTCGCTGCATAGACGCGTTCCTTGAAGGGGTCGATGCCGATACGCTGGTAGGTGTCGATGAAACGctct                                                      <  3:7640/149‑1 (MQ=255)
                                                                                                 gTCAGCTCGTCGAACGAGGTTTCCTTGGGCAGCAGGTGCCAGGTTTCGTCGACGATCTGGTTGTTCTTAATGATTCGCTGCATAGACGCGTTCCTTGAAGGGGTCGATGCCGATACGCTGGTAGGTGTCGAt                                                                >  7:176218/1‑132 (MQ=255)
                                                                                                 gTCAGCTCGTCGAACGAGGTTTCCTTGGGCAGCAGGTGCCAGGTTTCGTCGACGATCTGGTTGTTCTTAATGATTCGCTGCATAGACGCGTTCCTTGAAGGGGTCGATGCCGATACGCTGGTAGGTGTCGAt                                                                <  8:176218/132‑1 (MQ=255)
                                                                                                 gTCAGCTCGTCGAACGAGGTTTCCTTGGGCAGCAGGTGCCAGGTTTCGTCGACGATCTGGTTGTTCTTAATGATTCGCTGCATAGACGCGTTCCTTGAAGGGGTCGATGCCGATACGCTGGTAGGTGTCGATGAAACGCTCTTCCTCgg                                               >  6:50546/1‑149 (MQ=255)
                                                                                                             aaCGAGGTTTCCTTGGGCAGCAGGTGCCAGGTTTCGTCGACGATCTGGTTGTTCTTAATGATTCGCTGCATAGACGCGTTCCTTGAAGGGGTCGATGCCGATACGCTGGTAGGTGTCGATGAAACGCTCTTCCTCGGTACGTTGCTCGa                                   <  8:87913/149‑1 (MQ=255)
                                                                                                                    tttCCTTGGGCAGCAGGTGCCAGGTTTCGTCGACGATCTGGTTGTTCTTAATGATTCGCTGCATAGACGCGTTCCTTGAAGGGGTCGATGcc                                                                                     >  7:187382/1‑92 (MQ=255)
                                                                                                                    tttCCTTGGGCAGCAGGTGCCAGGTTTCGTCGACGATCTGGTTGTTCTTAATGATTCGCTGCATAGACGCGTTCCTTGAAGGGGTCGATGcc                                                                                     <  8:187382/92‑1 (MQ=255)
                                                                                                                        cTTGGGCAGCAGGTGCCAGGTTTCGTCGACGATCTGGTTGTTCTTAATGATTCGCTGCATAGACGCGTTCCTTGAAGGGGTCGATGCCGATACGCTGGTAGGTGTCGATGAAACGCTCTTCCTCGGTACGTTGCTCGACGTACACGGCg                        <  2:132894/149‑1 (MQ=255)
                                                                                                                           gggCAGCAGGTGCCAGGTTTCGTCGACGATCTGGTTGTTCTTAATGATTCGCTGCATAGACGCGTTCCTTGAAGGGGTCGATGCCGATACGCTGGTAGGTGTCGATGAAACGCTCTTCCTCGGTACGTTGCTCGACGTACACGGCGAt                      >  2:56540/1‑148 (MQ=255)
                                                                                                                              cagcagGTGCCAGGTTTCGTCGACGATCTGGTTGTTCTTAATGATTCGCTGCATAGACGCGTTCCTTGAAGGGGTCGATGCCGAAACGCTGGTAGGTGTCGATGAAACGCTCTTCCTCGGTACGTTGCTCGACGTACACGGCGATCAGc                  >  2:2328/1‑149 (MQ=255)
                                                                                                                               agcagGTGCCAGGTTTCGTCGACGATCTGGTTGTTCTTAATGATTCGCTGCATAGACGCGTTCCTTGAAGGGGTCGATGCCGATACGCTGGTAGGTGTCGATGAAACGCTCTTCCTCGGTACGTTGCTCGAAGTACACGGCGATCAGCt                 >  8:245022/1‑149 (MQ=255)
                                                                                                                                    gTGCCAGGTTTCGTCGACGATCTGGTTGTTCTTAATGATTCGCTGCATAGACGCGTTCCTTGAAGGGGTCGATGCCGATACGCTGGTAGGTGTCGATGAAACGCTCTTCCTCGGTACGTTGCTCGACGTACACGGCGATCAGCTTCTCg            <  4:206717/149‑1 (MQ=255)
                                                                                                                                         aGGTTTCGTCGACGATCTGGTTGTTCTTAATGATTCGCTGCATAGACGCGTTCCTTGAAGGGGTCGATGCCGATACGCTGGTAGGTGTCGATGAAACGCTCTTCCTCGGTACGTTGCTCGACGTACACGGCGATCAGCTTCTCGATcac       <  8:113748/149‑1 (MQ=255)
                                                                                                                                               cgtcgACGATCTGGTTGTTCTTAATGATTCGCTGCATAGACGCGTTCCTTGAAGGGGTCGATGCCGATACGCTGGTAGGTGTCGATGAAACGCTCTTCCTCGGTACGTTGCTCGACGTACACGGCGATCAGCTTCTCGATCACATCgg  <  6:281074/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GCTGTCCAGCCATACACCCAGGCCGCCATCGCGGGCCTTGAGCACGTGGGCATGGTCGCGCCACATTTGCAGCGGCACGATGTAGTCGTCGCAGTTGGTCAGCTCGTCGAACGAGGTTTCCTTGGGCAGCAGGTGCCAGGTTTCGTCGACGATCTGGTTGTTCTTAATGATTCGCTGCATAGACGCGTTCCTTGAAGGGGTCGATGCCGATACGCTGGTAGGTGTCGATGAAACGCTCTTCCTCGGTACGTTGCTCGACGTACACGGCGATCAGCTTCTCGATCACATCGG  >  NC_002947/2707393‑2707683

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