Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,165,959 A→G 7.3% G66G (GGT→GGC)  PP_5483 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,165,9590AG 7.3% 155.1 / 3.0 55G66G (GGT→GGC) PP_5483hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (22/29);  new base G (2/2);  total (24/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

CTGCACAAAACCAGATTGACCGGCAGCCAGCTCAATTAATCGCCAAGCTTCAATGATGGCGGCTCGGGGTGAGACTTCAGCGATTCGCTCCAATTGCCTAAGTTCTGGGGGAGTGCTCGGAGCTTTAGGTAGGGTGGGACCGATACTGTCTTCCAGTTTTGCAAGAGTTTCTTCTGTTTTGCTTAACTCTTCGGAAAAGCTTGCTTCAAACTCCTTGAATTTTAGGCTCGATACCCGTGTGAAAACGGCTCGGATTTCCGGGCGGAAAAGGTACGCAAGCCACAAAA  >  NC_002947/2165821‑2166107
                                                                                                                                          |                                                                                                                                                    
cTGCACAAAACCAGATTGACCGGCAGCCAGCTCAATTAATCGCCAAGCTTCAATGATGGCGGCTCGGGGTGAGACTTCAGCGATTCGCTCCAATTGCCTAAGTTCTGGGGGAGTGCTCGGAGCTTTAGGTAGGGTGGGACCGATACTg                                                                                                                                             >  5:284804/1‑148 (MQ=255)
   cacaAAACCAGATTGACCGGCAGCCAGCTCAATTAATCGCCAAGCTTCAATGATGGCGGCTCGGGGTGAGACTTCAGCGATTCGCTCCAATTGCCTAAGTTCTGGGGGAGTGCTCGGAGCTTTAGGTAGGGTGGGACCGATACTGTCtt                                                                                                                                         <  5:308181/149‑1 (MQ=255)
    acaAAACCAGATTGACCGGCAGCCAGCTCAATTAATCGCCAAGCTTCAATGATGGCGGCTCGGGGTGAGACTTCCGCGATTCGCTCCAATTGCCTAAGTTCTGGGGGAGTGCTCGGAGCTTTAGGTAGGGTGGGACCGATACTGTCTTc                                                                                                                                        <  6:260256/149‑1 (MQ=255)
    acaAAACCAGATTGACCGGCAGCCAGCTCAATTAATCGCCAAGCTCCAATGATGGCGGCTCGGGGTGAGACTTCAGCGATTCGCTCCAATTGCCTAAGTTCTGGGGGAGTGCTCGGAGCTTTAGGTAGGGTGGGACCGATACTGTCTTc                                                                                                                                        >  5:265099/1‑149 (MQ=255)
     caAAACCAGATTGACCGGCAGCCAGCTCAATTAATCGCCAAGCTTCAATGATGGCGGCTCGGGGTGAGACTTCAGCGATTCGCTCCAATTGCCTAAGTTCTGGGGGAGTGCTCGGAGCTTTAGGTAGGGTGGGACCGATACTGTCTTc                                                                                                                                        >  6:301578/1‑148 (MQ=255)
      aaaaCCAGATTGACCGGCAGCCAGCTCAATTAATCGCCAAGCTTCAATGATGGCGGCTCGGGGTGAGACTTCAGCGATTCGCTCCAATTGCCTAAGTTCTGGGGGAGTGCTCGGAGCTTTAGGTAGGGTGGGACCGATACTGTCTTc                                                                                                                                        >  4:29953/1‑147 (MQ=255)
      aaaaCCAGATTGACCGGCAGCCAGCTCAATTAATCGCCAAGCTTCAATGATGGCGGCTCGGGGTGAGACTTCAGCGATTCGCTCCAATTGCCTAAGTTCTGGGGGACTGCTCGGAGCTTTAGGTAGTGTGGGACCGATAATGTCTTCCa                                                                                                                                      >  8:303408/1‑149 (MQ=255)
        aaCCAGATTGACCGGCAGCCAGCTCAATTAATCGCCAAGCTTCAATGATGGCGGCTCGGGGTGAGACTTCAGCGATTCGCTCCAATTGCCTAAGTTCTGGGGGAGTGCTCGGAGCTTTAGGTAGGGTGGGACCGATACTGTCTTCCAGt                                                                                                                                    <  8:43209/149‑1 (MQ=255)
         aCCAGATTGACCGGCAGCCAGCTCAATTAATCGCCAAGCTTCAATGATGGCGGCTCGGGGTGAGACTTCAGCGATTCGCTCCAATTGCCTAAGTTCTGGGGGAGTGCTCGGAGCTTTAGGTAGGGTGGGACCGATACTGTCTTCCAGtt                                                                                                                                   >  4:195716/1‑149 (MQ=255)
             gATTGACCGGCAGCCAGCTCAATTAATCGCCAAGCTTCAATGATGGCGGCTCGGGGTGAGACTTCAGCGATTCGCTCCAATTGCCTAAGTTCTGGGGGAGTGCTCGGAGCTTTAGGTAGGGTGGGACCGATACTGTCTTCCAGTTTTGc                                                                                                                               <  5:246611/149‑1 (MQ=255)
              aTTGACCGGCAGCCAGCTCAATTAATCGCCAAGCTTCAATGATGGCGGCTCGGGGTGAGACTTCAGCGATTCGCTCCAATTGCCTAAGTTCTGGGGGAGTGCTCGGAGCTTTAGGTAGGGTGGGACCGATACTGTCTTCCAGTTTTGCa                                                                                                                              <  5:301578/149‑1 (MQ=255)
               ttGACCGGCAGCCAGCTCAATTAATCGCCAAGCTTCAATGATGGCGGCTCGGGGTGAGACTTCAGCGATTCGCTCCAATTGCCTAAGTTCTGGGGGAGTGCTCGGAGCTTTAGGTAGGGTGGGACCGATACTGTCTTCCAGTTTTGCaa                                                                                                                             <  3:48847/149‑1 (MQ=255)
                            agcTCAATTAATCGCCAAGCTTCAATGATGGCGGCTCGGGGTGAGACTTCAGCGATTCGCTCCAATTGCCTAAGTTCTGGGGGAGTGCTCGGAGCTTTAGGTAGGGTGGGACCGATACTGTCTTCCAGTTTTGCAAGAGTTTCTTCTGt                                                                                                                >  8:139929/1‑149 (MQ=255)
                             gcTCAATTAATCGCCAAGCTTCAATGATGGCGGCTCGGGGTGAGACTTCAGCGATTCGCTCCAATTGCCTAAGTTCTGGGGTAGTGCTCGGAGCTTTAGGTAGGGTGGGACCGATACTGTCTTCCAGTTTTGCAAGAGTTTCTTCTGtt                                                                                                               <  3:53330/149‑1 (MQ=255)
                              cTCAATTAATCGCCAAGCTTCAATGATGGCGGCTCGGGGTGAGACTTCAGCGATTCGCTCCAATTGCCTAAGTTCTGGGGGAGTGCTCGGAGCTTTAGGTAGGGTGGGACCGATACTGTCTTCCAGTTTTGCAAGAGTTTCTTCTGttt                                                                                                              <  1:206652/149‑1 (MQ=255)
                                 aaTTAATCGCCACGCTTCAATGATGCCGGCTCGGGGTGAGACTTCAGCGATTCGCTCCAATTGCCTAAGTTCTGGGGGAGTGCTCGGAGCTTTAGGTAGGGGGGGACCGATACTGTCTTCCAGTTTTGCAAGAGTTTCTTCTGTTTTGc                                                                                                           <  4:143216/149‑1 (MQ=255)
                                 aaTTAATCGCCAAGCTTCAATGATGGTGGCTCGGGGTGAGACTTCAGCGATTCGCTCCAATTGCCTAAGTTCTGGGGGAGTGCTCGGAGCTTTAGGTAGGGTGGGGCCGATACTGTCTTCCAGTTTTGCAAGAGTTTCTTCTGTTTTGc                                                                                                           >  6:123591/1‑149 (MQ=255)
                                 aaTTAATCGCCAAGCTTCAATGATGGTGGCTCGGGGTGAGACTTCAGCGATTCGCTCCAATTGCCTAAGTTCTGGGGGAGTGCTCGGAGCTTTAGGTAGGGTGGGGCCGATACTGTCTTCCAGTTTTGCAAGAGTTTCTTCTGTTTTGc                                                                                                           >  6:128813/1‑149 (MQ=255)
                                 aaTTAATCGCCAAGCTTCAATGATGGCGGCTCGGGGTGAGACTTCAGCGATTCGCTCCAATTGCCTAAGTTCTGGGGGAGTGCTCGGAGGTTTAGGTAGGGTGGGACCGATAGTGTCTTCCAGTTTTGCAAGAGTTTCTTCTGTTTTGc                                                                                                           >  8:270476/1‑149 (MQ=255)
                                 aaTTAATCGCCAAGCTCCAATGATGGCGGCTCGGGGTGAGACTTCAGCGATTCGCTCCAATTGCCTAAGTTCTGGGGGAGTGCTCGGAGCTTTAGGTAGGGTGGGACCGATACTGTCTTCCAGTTTTGCAAGAGTTTCTTCTGtttt                                                                                                             <  6:265099/147‑1 (MQ=255)
                                          ccAAGCTTCAATGATGGCGGCTCGGGGTGAGACTTTAGCGATTCGCTCCAATTGCCTAAGTTCTGGGGGAGTGCTCGGAGCTTTAGGTAGGGGGGGACCGATACTGTCTTCCAGTTTTGCAAGAGTTTCTTCTGTTTTGCTTAACTCtt                                                                                                  <  2:97884/149‑1 (MQ=255)
                                           cAAGCTTCAATGATGGTGGCTCGGGGTGAGACTTCAGCGATTCGCTCCAATTGCCTAAGTTCTGGGGGAGTGCTCGGAGCTTTAGGTAGGGTGGGGCCGATACTGTCTTCCAGTTTTGCAAGAGTTTCTTCTGTTTTGCTTAACTCTTc                                                                                                 <  5:123591/149‑1 (MQ=255)
                                           cAAGCTTCAATGATGGTGGCTCGGGGTGAGACTTCAGCGATTCGCTCCAATTGCCTAAGTTCTGGGGGAGTGCTCGGAGCTTTAGGTAGGGTGGGGCCGATACTGTCTTCCAGTTTTGCAAGAGTTTCTTCTGTTTTGCTTAACTCTTc                                                                                                 <  5:128813/149‑1 (MQ=255)
                                            aaGCTTCACTGATGGCGGCTCGGGGTGAGACTTCAGCGATTCGCTCCAATTGCCTAAGTGCTGGGGGAGTGCTCGGAGCTTTAGGTAGGGTGGGACCGATACTGTCTTCCAGTTTTGCAAGAGTTTCTTCTGTTTTGCTTAACTCTTCg                                                                                                <  8:242669/149‑1 (MQ=255)
                                            aaGCTTCAATGATGGCGGCTCGGGGTGAGACTTCAGCGATTCGCTCCAATTGCCTAAGTTCTGGGGGAGTGCTCGGAGCTTTAGGTAGGGTGGGACCGATACTGTCTTCCAGTTTTGCAAGAGTTTCTTCTGTTTTGCTTAACTCTTCg                                                                                                <  3:195716/149‑1 (MQ=255)
                                            aaGCTTCAATGATGGCGGCTCGGGGTGAGACTTCAGCGATTCGCTCCAATTGCCTAAGTTCTGGGGGAGTGCTCGGAGCTTTAGGTAGGGTGGGACCGATACTGTCTTCCAGTTTTGCAAGAGTTTCTTCTGTTTTGCTTAACTCTTCg                                                                                                >  3:161470/1‑149 (MQ=255)
                                            aaGCTTCAATGATGGCGGCTCGGGGTGAGACTTCAGCGATTCGCTCCAATTGCCTAAGTTCTGGGGGAGTGCTCGGAGCTTTAGGTAGGGTGGGACCGATACTGTCTTCCAGTTTTGCAAGAGTTTCTTCTGTTTTGCTTAACTCTTCg                                                                                                <  3:29953/149‑1 (MQ=255)
                                            aaGCTTCAATGATGGCGGCTCGGGGTGAGACTTCAGCGATTCGCTCCAATTGCCTAAGTTCTGGGGGAGTGCTCGGAGCTTTAGGTAGGGTGGGACCGATACTGTCTTCCAGTTTTGCAAGAGTTTCTTCTGTTTTGCTTAACTCTTCg                                                                                                <  5:77556/149‑1 (MQ=255)
                                                ttCAATGATGGCGGCTCGGGGTGAGACTTCAGCGATTCGCTCCAATTGCCTAAGTTCTGGGGGAGTGCTCGGAGCTTTAGGTAGGGTGGGACCGATACTGTCTTCCAGTTTTGCAAGAGTTTCTTCTGTTTTGCTTAACTCTTCGGag                                                                                             >  2:235778/1‑147 (MQ=255)
                                                 tCAATGATGGCGGCTCGGGGTGAGACTTCAGCGATTCGCTCCAATTGCCTAAGTTCTGGGGGAGTGCTCGGAGCTTTAGGTAGGGTGGGACCGATACTGTCTTCCAGTTTTGCAAGAGTTTCTTCTGTTTTGCTTAACTCTTCGGaa                                                                                             >  4:84985/1‑147 (MQ=255)
                                                     tgatgGCGGCTCGGGGTTAGACTTCAGCGATTCGCTCCAATTGCCTAAGTTCTGGGGGAGTGCTCGGAGCTTTAGGTAGGGTGGGACCGATACTGTCTTCCAGTTTTGCAAGAGTTTCTTCTGTTTTGCTTAACTCTTCGGAAAAgctt                                                                                       <  6:284804/149‑1 (MQ=255)
                                                     tgatgGCGGCTCGGGGTGAGACTTCAGCGATTCGCTCCAATTGCCTAAGTTCTAGGGGAGTGCTCGGAGCTTTAGGTAGGGTGGGACCGATACTGTCTTCCAGTTTTGCAAGAGTTTCTTCTGTTTTGCTTAACTCTTCGGAAAAgctt                                                                                       >  8:14684/1‑149 (MQ=255)
                                                       atgGCGGCTCGGGGTGAGACTTCAGCGATTCGCTCCAATTGCCTAAGTTCTAGGGGAGTGCTCGGAGCTTTAGGTAGGGTGGGACCGATACTGTCTTCCAGTTTTGCAAGAGTTTCTTCTGTTTTGCTTAACTCTTCGGAAAAgcttg                                                                                      <  7:14684/148‑1 (MQ=255)
                                                           cggcTCGGGGTGAGACTTCAGCGATTCGCTCCAATTGCCTAAGTTCTGGGGGAGTGCTCGGAGCTTTAGGTAGGGTGGGACCGATACTGTCTTCCAGTTTTGCAAGAGTTTCTTCTGTTTTGCTTAACTCTTCGGAAAAGCTTGCTTc                                                                                  >  4:80299/1‑148 (MQ=255)
                                                              cTCGGGGTGAGACTTCAGCGATTCGCTCCAATTGCCTAAGTTCTGGGGGAGTGCTCGGAGCTTTAGGTAGGGTGGGACCGATACTGTCTTCCAGTTTTGCAaga                                                                                                                           <  8:170473/104‑1 (MQ=255)
                                                                        gaCCTCAGCGATTCGCTCCAATTGCCTAAGTTCTGTGGGCGTGCTCGGAGCTTTAGGTAGGGTGGGACCGATACTGTCTACCAGTTTTGCAAGAGTTTCTTCTGTTTTGCTTAAGTCTTCGGAAAAGCTTGCTTCAAACTCCTTGAAt                                                                     <  8:240157/148‑1 (MQ=255)
                                                                             cAGCGATTCGCTCCAATTGCCTAAGTTCTGGGGGAGTGCTCGGAGCTTTAGGTAGGGTGGGACCGATACTGTCTTCCAGTTTTGCAAGAGTTTCTTCTGTTTTGCTTAACTCTTCGGAAAAGCTTGCTTCAAACTCCTTGAACTTTAg                                                                <  7:303408/148‑1 (MQ=255)
                                                                              aGCTATTCGCTCCAATTGCCTAAGTTCTGGGGGAGTGCTCGGAGCTTTAGGTAGGGTGGGACCGATACTGTCTTCCAGTTTTGCAAGAGTTTCTTCTGTTTTGCTTAACTCTTCGGAAAAGCTTGCTTCAAACTCCTTGAATTTTAGGc                                                              <  3:317616/149‑1 (MQ=255)
                                                                                       cTCCAATTGCCTAAGTTCTGGGGGAGTGCTCGGAGCTTTAGGTAGGGTGGGACCGATACTGTCTTCCAGTTTTGCAAGAGTTTCTTCTGTTTTGCTTAACTCTTCGGAAAAGCTTGCTTCAAACTCCTTGAATTTTAGGCTCGATAccc                                                     <  7:139929/149‑1 (MQ=255)
                                                                                         ccAATTGCCTAAGTTCTGGGGGAGTGCTCGGAGCTTTAGGTAGGGTGGGACCGATACTGTCTTCCAGTTTTGCAAGAGTTTCTTCTGTTTTGCTTAACTCTTCGGAAAAGCTTGCTTCAAACTCCTTGAATTTTAGGCTCGATACCCgt                                                   <  4:161566/149‑1 (MQ=255)
                                                                                          cAATTGCCTAAGTTCTGGGGGAGTGCTCGGAGCTTTAGGTAGGGTGGGACCGATACTGTCTTCCa                                                                                                                                      <  7:176980/65‑1 (MQ=255)
                                                                                          cAATTGCCTAAGTTCTGGGGGAGTGCTCGGAGCTTTAGGTAGGGTGGGACCGATACTGTCTTCCa                                                                                                                                      >  8:176980/1‑65 (MQ=255)
                                                                                          aaaTTGCCTAAGTTCTGGGGGAGTGCTCGGAGCTTTAGGTAGGGTGGGACCGATACTGTCTTCCAGTTTTGCAAGAGTTTCTTCTGTTTTGCTTAACTCTTCGGAAAAGCTTGCTTCAAACTCCTTGAATTTTAGGCTCGATACCCgtg                                                  >  2:69559/2‑149 (MQ=255)
                                                                                                ccTAAGTTCTGGGGGAGTGCTCGGAGCTTTAGGTAGGGTGGGACCGATACTGTCTTCCAGTTTTGCAAGAGTTTCTTCTGTTTTGCTTAACTCTTCGGAAAAGCTTGCTTCAAACTCCTTGAATTTTAGGCTCGATACCCGTGTGaaa                                             >  8:10275/1‑148 (MQ=255)
                                                                                                      ttCTGGGGGAGTGCTCGGAGCTTTAGGTAGGGTGGGACCGATACTGTCTTCCAGTTTTGCAAGAGTTTCTTCTGTTTTGCTTAACTCTTCGGAAAAGCTTGCTTCAAACTCCTTGAATTTTAGGCTCGATACCCGTGTGAAAACGGCTc                                      >  7:17825/1‑149 (MQ=255)
                                                                                                      ttCTGGGGGAGTGCTCGGAGCTTTAGGTAGGGTGGGACCGATACTGTCTTCCAGTTTTGCAAGAGTTTCTTCTGTTTTGCTTAACTCTTCGGAAAAGCTTGCTTCAAACTCCTTGAATTTTAGGCTCGATACCCGTGTGAAAACGGCTc                                      <  3:80299/149‑1 (MQ=255)
                                                                                                         tGGGGGAGTGCTCGGAGCTTTAGGTAGGGTGGGACCGATACTGTCTTCCAGTTTTGCAAGAGTTTCTTCTGTTTTGCTTAACTCTTCGGAAAAGCTTGCTTCAAACTCCTTGAATTTTAGGCTCGATACCCGTGTGAAAACGGCTCGGa                                   >  3:290359/1‑149 (MQ=255)
                                                                                                                   cTCGGAGCTTTAGGTAGGGTGGGACCGATACTGTCTTCCAGTTTTGCAAGAGTTTCTTCTGTTTTGCTTAACTCTTCGGAAAAGCTTGCTTCAAACTCCTTGTATTTTAGGCTCGATACCCGTGTGAAAACGGCTCGGATTTCCGGGc                          <  1:69559/148‑1 (MQ=255)
                                                                                                                    tCGGAGCTTTAGGTAGGGTGGGACCGATACTGTCTTCCAGTTTTGCAAGAGTTTCTTCTGTTTTGCTTAACTCTTCGGAAAAGCTTGCTTCAAACTCCTTGAATTTTAGGCTCGATACCCGTGTGAAAACGGCTCGGATTTCCGGGCg                         >  4:99338/1‑148 (MQ=255)
                                                                                                                              aggtaggGTGGGACCGATACTGTCTTCCAGTTTTGCAAGAGTTTCTTCTGTTTTGCTTAACTCTTCGGAAAAGCTTGCTTCAAACTCCTTGAATTTTAGGCTCGATACCCGTGTGAAAACGGCTCGGATTTCCGGGCGGAAAAGGTAc               <  6:8491/148‑1 (MQ=255)
                                                                                                                               ggtaggGTGGGACCGATACTGTCTTCCAGTTTTGCAAGAGTTTCTTCTGTTTTGCTTAACTCTTCGGAAAAGCTTGCTTCAAACTCCTTGAATTTTAGGCTCGATACCCGTGTGAAAACGGCTCGGATTTCCGGGCGGAAAAGGTACg              <  7:10275/148‑1 (MQ=255)
                                                                                                                               ggtaggGTGGGACCGATACTGTCTTCCAGTTTTGCAAGAGTTTCTTCTGTTTTGCTTAACTCTTCGGAAAAGCTTGCTTCAAACTCCTTGAATTTTAGGCTCGATACCCGAGTGAAAACGGCTCGGATTTCCGGGCGGAAAAGGTACg              >  5:186395/1‑148 (MQ=255)
                                                                                                                                  aggGTGGGACCGATACTGTCTTCCAGTTTTGCAAGAGTTTCTTCTGTTTTGCTTAACTCTTCGGAAAAGCTTGCTTCAAACTCCTTGAATTTTAGGCTCGATACCCGTGTGAAAACGGCTCGGATTTCCGGGCGGAAAAGGTACGCAAg          >  6:107353/1‑149 (MQ=255)
                                                                                                                                     gTGGGACCGATACTGTCTTCCAGTTTTGCAAGAGTTTCTTCTGTTTTGCTTAACTCTTCGGAAAAGCTTGCTTCAAACTCCTTGAATTTTAGGCTCGATACCCGTGTGAAAACGGCTCgg                                    <  5:331707/120‑1 (MQ=255)
                                                                                                                                     gTGGGACCGATACTGTCTTCCAGTTTTGCAAGAGTTTCTTCTGTTTTGCTTAACTCTTCGGAAAAGCTTGCTTCAAACTCCTTGAATTTTAGGCTCGATACCCGTGTGAAAACGGCTCgg                                    >  6:331707/1‑120 (MQ=255)
                                                                                                                                          aCCGATACTGTCTTCCAGTTTTGCAAGAGTTTCTTCTGTTTTGCTTAACTCTTCGGAAAAGCTTGCTTCAAACTCCTTGAATTTTAGGCTCGATACCCGTGTGAAAACGGCTCGGATTTCCGGGCGGAAAAGGTACGCAAGCCACaaaa  <  3:84985/149‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                    
CTGCACAAAACCAGATTGACCGGCAGCCAGCTCAATTAATCGCCAAGCTTCAATGATGGCGGCTCGGGGTGAGACTTCAGCGATTCGCTCCAATTGCCTAAGTTCTGGGGGAGTGCTCGGAGCTTTAGGTAGGGTGGGACCGATACTGTCTTCCAGTTTTGCAAGAGTTTCTTCTGTTTTGCTTAACTCTTCGGAAAAGCTTGCTTCAAACTCCTTGAATTTTAGGCTCGATACCCGTGTGAAAACGGCTCGGATTTCCGGGCGGAAAAGGTACGCAAGCCACAAAA  >  NC_002947/2165821‑2166107

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