Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,366,601 A→T 8.5% M45K (ATG→AAG)  PP_3843 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,366,6010AT 8.5% 132.6 / 3.1 47M45K (ATG→AAG) PP_3843hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (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 = 7.76e-01

AAAGCCTGCCGCAGGCGCACCTTCAGCGGCGACACGTCGGGCCCGGCCAACGCCCCACGCCAGCGCCCTGGTCATCGACTCGCCAGGGCGAGAATCGAAGGCCTCCTCGAGGA‑TAGCCATGCCACCAGGGGCGTACACCCCGATGAACATCTGCGTTTTGCCTGTCCGCGACAGCCGCGTCTGAACATCAATGTATGTCCCATCGTCGAGTGTCTCGTCGTGAGTTCGGTGATGAAGCGTCGGGTCAGCCCAGGCCCAGAAAACATCACCGCGTATTCTCATGAGACGCTCCTACGA  >  NC_002947/4366453‑4366749
                                                                                                                                                     |                                                                                                                                                    
aaaGCCTGACGCAGGCGCACCTTCAGCGGCGACACGTCGGGCCCGGCCAACGCCCCACGCCAGCGCCCTGGTCATCGACTCGCCAGGGCGAGAATCGAAGGCCTCCTCGAGGA‑TAGCCATGCCACCAGGGGCGGACACCCCGATGAACa                                                                                                                                                      >  6:212697/1‑149 (MQ=255)
          gCAGGCGCACCTTCAGCGGCGACACGTCGGGCCCGGCCAACGCCCCACGCCAGCGCCCTGGTCATCGACTCGCCAGGGCGAGAATCGAAGGCCTCCTCGAGGA‑TAGCCATGCCACCAGGGGCGTACACCCCGATGAACATCTGCGtttt                                                                                                                                            <  3:264269/149‑1 (MQ=255)
          gCAGGCGCACCTTCAGCGGCGACACGTCGGGCCCGGCCAACGCCCCACGCCAGCGCCCTGGTCATCGACTCGCCAGGGCGAGAATCGAAGGCCTCCTCGAGGA‑TAGCCATGCCACCAGGGGCGTACACCCCGATGAACATCTGCGtttt                                                                                                                                            <  3:259469/149‑1 (MQ=255)
           cAGGCGCACCTTCAGCGGCGACACGTCGGGCCCGGCCAACGCCCCACGCCAGCGCCCTGGTCATCGACTCGCCAGGGCGAGAATCGAAGGCCTCCTCGAGGA‑TAGCCATGCCACCAGGGGCGTACACCCCGATGAACATCTGCGTTTTg                                                                                                                                           <  1:264099/149‑1 (MQ=255)
                               acacGTCGGGCCCGGCCAACGCCCCACGCCAGCGCCCTGGTCATCGACTCGCCAGGTCGAGAATCGAAGGCCTCCTCGAGGA‑TAGCCATGCCACCAGGGGCGTACACCCCGATGAACATCTGCGTTTTGCCTGTCCGCGACAGCCGCGt                                                                                                                       <  3:256917/149‑1 (MQ=255)
                                 acGGCGGGCCCGGCCAACGCCCCACGCCAGCGACCTGGTCATCGACTCGCCAGGTCGAGAATCGAAGGCCTCCTCGAGTA‑TAGCCATGCCACCAGGGGCGTACACCACGATGAACATCTGCGTTTTGCCTGTCCGCGACAGCCGCGTCt                                                                                                                     <  2:360229/149‑1 (MQ=255)
                                         cccGGCCAACGCCCCACGCCAGCGCCCTGGTCATCGACTCGCCAGGGCGAGAATCGAAGGCCTCCTCGAGGA‑TAGCCATGCCACCAGGGGCGTACACCCCGATGAACATCTGCGTTTTGCCTGTCCGCGACAGCCGCGTCTGAACATCa                                                                                                             <  3:96810/149‑1 (MQ=255)
                                               cAACGCCCCACGCCAGCGCCCTGGTCATCGACTCGCCAGGGCGAGAATCGAAGGCCTCCTCGAGGA‑TAGCCATGCCACCAGGGGCGTACACCCCGATGAACATCTGCGTTTTGCCTGTCCGCGACAGCCGCGTCTGAACATCAatgtat                                                                                                       >  3:292572/1‑149 (MQ=255)
                                                           ccAGCGCCCTGGTCATCGACTCGCCAGGGCGAGAATCGAAGGCCTCCTCGAGGA‑TAGCCATGCCACCAGGGGCGTACACCCCGATGAACATCTGCGTTTTGCCTGTCCGCGACAGCCGCGTCTGAACATCAATGTATGTCCCAtcgtc                                                                                            >  1:118015/1‑148 (MQ=255)
                                                            cAGCGCCCTGGTCATCGACTCGCCAGGGCGAGAATCGAAGGCCTCCTCGAGGA‑TAGCCATGCCACCAGGGGCGTACACCCCGATGAACATCTGCGTTTTGCCTGTCCGCGACAGCCGCGTCTGAACATCAATGTATGTCCCAtcgtc                                                                                            >  7:76254/1‑147 (MQ=255)
                                                                  ccTGGTCATCGACTCGCCAGGGCGAGAATCGAAGGCCTCCTCGAGGA‑TAGCCATGCCACCAGGGGCGTACAACCCGATGAACATCTGCGTTTTGCCTGTCCGCGACACCCGCGTCTGAACATCAATGTATGTCCCATAGTCGAGTGtc                                                                                     >  2:208501/1‑148 (MQ=255)
                                                                      gTCATCGACTCGCCAGGGCGAGAATCGAAGGCCTCCTCGAGGA‑TAGCCATGCCACCAGGGGCGTACACCCCGATGAACATCTGCGTTTTGCCTGTCCGCGACAGCCGCGTCTGAACATCAATGTATGTCCCATCGTCGAGTGTCtcgtc                                                                                <  8:76254/149‑1 (MQ=255)
                                                                      gTCATCGACTCGCCAGGGCGAGAATCGAAGGCCTCCTCGAGGA‑TAGCCATGCCACCAGGGGCGTACACCCCGATGAACATCTGCGTTTTGCCTGTCCGCGACAGCCGCGTCTGAACATCAATGTATGTCCCATCGTCGAGTGTCtcgtc                                                                                >  5:197515/1‑149 (MQ=255)
                                                                      gTCATCGACTCGCCAGGGCGAGAATCGAAGGCCTCCTCGAGGA‑TAGCCATGCCACCAGGGGCGTACACCCCGATGAACATCTGCGTTTTGCCTGTCCGCGACAGCCGCGTCTGAACATCAATGTATGTCCCATCGTCGAGTGTCtcgtc                                                                                <  6:329269/149‑1 (MQ=255)
                                                                         aTCGACTCGCCAGGGCGAGAATCGAAGGCCTCCTCGAGGA‑TAGCCATGCCACCAGGGGCGTACACCCCGATGAACATCTGCCTTTTGCCTGTCCGCGACAGCCGCGTCTg                                                                                                                    <  4:189532/110‑1 (MQ=255)
                                                                               tCGCCAGGGCGAGAATCGAAGGCCTCCTCGAGGA‑TAGCCATGCCACCAGGGGCGTACACCCCGATGAACATCTGCGTTTTGCCTGTCCGCGACAGCCGCGTCTGAACATCAATGTATGTCCCATCGTCGAGTGTCTCGTCGTGAGTTCg                                                                       >  5:178638/1‑149 (MQ=255)
                                                                                 gCCAGGGCGAGAATCGAAGGCCTCCTCGAGGA‑TAGCCATGCCACCAGGGGCGTACACCCCGATGAACATCTGCGTTTTGCCTGTCCGCGAc                                                                                                                               <  6:93833/91‑1 (MQ=255)
                                                                                 gCCAGGGCGAGAATCGAAGGCCTCCTCGAGGA‑TAGCCATGCCACCAGGGGCGTACACCCCGATGAACATCTGCGTTTTGCCTGTCCGCGAc                                                                                                                               >  5:93833/1‑91 (MQ=255)
                                                                                    aGGGCGAGAATCGAAGGCCTCCTCGAGGATTAGCCATGCCACCAGGGGCGTACACCCCGATGAACATCTGCGTTTTGCCTGTCCGCGACAGCCGCGTCTGAACATCAATGTATGTCCCATCGTCGAGTGTCTCGTCGTGAGTTCGGtga                                                                   >  1:26951/1‑149 (MQ=255)
                                                                                    aGGGCGAGAATCGAAGGCCTCCTCGAGGA‑TAGCCATGCCACCAGGGGCGTACACCCCGATGAACATCTGCGTTTTGCCTGTCCGCGACAGCCGCGTCTGAACATCAATGTATGTCCCATCGTCGAGTGTCTCGTCGTGAGTTCGGtgat                                                                  >  5:288144/1‑149 (MQ=255)
                                                                                       tcGAGAATCGAATGCCTCCTCGAGGA‑TATCCATGCCACCAGGGGCTTACACCCCGATGAAAATCTGCGTTTTGCCTGTCCGCTACAGCCGCGTCTGAACATCAATGTATGTCCCATCGTCGAGTTTCTCGTCGTGAGTTCGGTGATGaa                                                               <  1:208501/148‑1 (MQ=255)
                                                                                       gCGAGAATCGAAGGCCTCCTCGAGGA‑TAGCCATGCCACCAGGGGCGTACACCCCGATGAACATCTGCGTTTTGCCTGTCCGCGACAGCCGCGTCTGAACATCAATGTATGTCCCATCGTCGAGTGTCTCGTCGCGAGTTCGGTGATGaa                                                               <  3:293867/149‑1 (MQ=255)
                                                                                         gagaATCGAAGGCCTCCTCGAGGA‑TAGCCATGCCACCAGGGGCGTACACCCCGATGAACATCTGCGTTTTGCCTGTCCGCGACAGCCGCGTCTGAACATCAATGTATGTCCCATCGTCGAGTGtctc                                                                                   <  1:14395/127‑1 (MQ=255)
                                                                                         gagaATCGAAGGCCTCCTCGAGGA‑TAGCCATGCCACCAGGGGCGTACACCCCGATGAACATCTGCGTTTTGCCTGTCCGCGACAGCCGCGTCTGAACATCAATGTATGTCCCATCGTCGAGTGtctc                                                                                   >  2:14395/1‑127 (MQ=255)
                                                                                         gagaATCGAAGGCCTCCTCGAGGA‑TAGCCATGCCACCAGGGGCGTACACCCCGATGAACATCTGCGTTTTGCCTGTCCGCGACAGCCGCGTCTGAACATCAATGTATGTCCCATCGTCGAGTGTCTCGTCGTGAGTTCGGTGATGAAGc                                                             >  7:131863/1‑149 (MQ=255)
                                                                                              tCGAAGGCCTCCTCGAGGA‑TAGCCATGCCACCAGGGGCGTACACCCCGATGAACATCTGCGTTTTGCCTGTCCGCGACAGCCGCGTCTGAACATCAATGTATGTCCCATCGTCGAGTGTCTCGTCGTGAGTTCGGTGATGAAGCGTCgg                                                        >  1:77884/1‑149 (MQ=255)
                                                                                                      ctcctcGAGGA‑TAGCCATGCCACCAGGGGCGTACACCCCGATGAACATCTGCGTTTTGCCTGTCCGCGACAGCCGCGTCTGAACATCAATGTATGTCCCATCGTCGAGTGTCTCGTCGTGAGTTCGGTGATGAAGCGTCGGGTCAgcc                                                 >  6:104958/1‑148 (MQ=255)
                                                                                                         ctcGAGGA‑TAGCCATGCCACCAGGGGCGTACACCCCGATGAACATCTGCGTTTTGCCTGTCCGCGACAGCCGCGTCTGAACATCAATGTATGTCCCATCGTCGAGTGTCTCGTCGTGAGTTCGGTGATGAAGCGTCGGGTCAgcccag                                              <  6:197515/148‑1 (MQ=255)
                                                                                                          tcGAGGA‑TAGCCATGCCACCAGGGGCGTACACCCCGATGAACATCTGCGTTTTGCCTGTCCGCGACAGCCGCGTCTGAACATCAATGTATGTCCCATCGTCGAGTGTCTCGTCGTGAGTTCGGTGATGAAGCGTCGGGTCAgcccaggc                                            >  6:246538/1‑149 (MQ=255)
                                                                                                              ggA‑TAGCCATGCCACCAGGGGCGTACACCCCGATGAACTTCTGCGTTTTGCCTGTCCGCGACAGCCGCGTCTGAACATCAATGTATGTCCCATCGTCGAGTGTCtcgtc                                                                                >  1:130761/1‑109 (MQ=255)
                                                                                                              ggA‑TAGCCATGCCACCAGGGGCGTACACCCCGATGAACTTCTGCGTTTTGCCTGTCCGCGACAGCCGCGTCTGAACATCAATGTATGTCCCATCGTCGAGTGTCtcgtc                                                                                <  2:130761/109‑1 (MQ=255)
                                                                                                                       aTGCCACCAGGGGCGTACACCCCGATGAACATCTGCGTTTTGCCTGTCCGCGACAGCCGCGTCTGAACATCAATGTATGTCCCATCGTCGAGTGTCTCGTCGTGAGTTCGGTGATGAAGCGTCGGGTCAGCCCAGGCCCAGAAAACAt                                 >  1:357753/1‑148 (MQ=255)
                                                                                                                       aTGCCACCAGGGGCGTACACCCCGATGAACATCTGCGTTTTGCCTGTCCGCGACAGCCGCGTCTGAACATCAATGTATGTCCCATCGTCGAGTGTCTCGTCGTGAGTTCGGTGATGAAGCGTCGGGTCAGCCCAGGCCCAGAAAACATc                                <  5:104958/149‑1 (MQ=255)
                                                                                                                          ccaccaGGGGCGTACACCCCGATGAACATCTGCGTTTTGCCTGTCCGCGACAGCCGCGTCTGAACATCAATGTATGTCCCATCGTCGAGTGTCTCGTCGTGAGTTCGGTGATGAAGCGTCGGGTCAGCCCAGGCCCAGAAAACATCAc                              <  6:178638/148‑1 (MQ=255)
                                                                                                                           caccaGGGGCGTACACCCCGATGAACATCTGCGTTTTGCCTGTCCGCGACAGCCGCGTCTGAACATCAATGTATGTCCCATCGTCGAGTGTCTCGTCGTGAGTTCCGTGATGAAGCGTCGGGTCAGCCCATGCCCACAAAACATCACcg                            >  8:346618/1‑149 (MQ=255)
                                                                                                                            accaGGGGCGTACACCCCGATGAACATCTGCGTTTTGCCTGTCCGCGACAGCCGCGTCTGAACATCAATGTATGTCCCATCGTCGAGTGTCTCGTCGTGAGTTCGGTGATGAAGCGTCGGGTCAGCCCAGGCCCAGAAAACATCACcgc                           >  6:306121/1‑149 (MQ=255)
                                                                                                                            accaGGGGCGTACACCCCGATGAACATCTGCGTTTTGCCTGTCCGCGACAGCCGCGTCTGAACATCAATGTATGTCCCATCGTCGAGAGTCTCGTCGTGAGTTCGGTGATGAAGCGTCGGGTCAGCCCAGGCCCAGAAAACATCACcgc                           >  7:79018/1‑149 (MQ=255)
                                                                                                                                     gTACACCCCGATGAACATCTGCGTTTTGCCTGTCCGCGACAGCCGCGTCTGAACATCAATGTATGTCCCATCGTCGAGTGTCTCGTCGTGAGTTCGGTGATGAAGCGTCGGGTCAGCCCAGGCTCAGAAAACATCACCGCGTATTCTCa                  >  5:116667/1‑149 (MQ=255)
                                                                                                                                      tACACCCCGATGAACATCTGCGTTTTGCCTGTCCGCGACAGCCGCGTCTGAACATCAATGTATGTCCCATCGTCGAGTGTCTCGTCGTGAGTTCGGTGATGAAGCGTCGGGTCAGCCCAGGCCCAGAAAACATCACCGCGTATTCTCAt                 <  2:357753/149‑1 (MQ=255)
                                                                                                                                      tACACCCCGATGAACATCTGCGTTTTGCCTGTCCGCGACAGCCGCGTCTGAACATCAATGTATGTCCCATCGTCGAGTGTCTCGTCGTGAGTTCGGTGATGAAGCGTCCGGTCAGCCCAGGCCCAGAAAACATCACCGCGTATTCTCAt                 >  6:89703/1‑149 (MQ=255)
                                                                                                                                        cacCCCGATGAACATCTGCGTTTTGCCTGTCCGCGACAGCCGCGTCTGAACATCAATGTATGTCCCATCGTCGAGTGTCTCGTCGTGAGTTCGGTGATGAAGCGTCGGGTCAGCCCAGGCCCAGAAAACATCACCGCGTATTCTCATg                <  4:292572/148‑1 (MQ=255)
                                                                                                                                        cacCCCGATGAACATCTGCGTTTTGCCTGTCCGCGACAGCCGCGTCTGAACATCAATGTATGTCCCATCGTCGAGTGTCTCGTCGTGAGTTCGGTGATGAAGCGTCGGGTCAGCCCAGGCCCAGAAAACATCACCGCGTATTCTCATg                <  6:288144/148‑1 (MQ=255)
                                                                                                                                        cacCCCGATGAACATCTGCGTTTTGCCTGTCCGCGACAGCCGCGTCTGAACATCAATGTATGTCCCATCGTCGAGTGTCTCGTCGTGAGTTCGGTGATGAAGCGTCGGGTCAGCCCAGGCCCAGAAAACATCACCGCGTATTCTCATg                >  6:10033/1‑148 (MQ=255)
                                                                                                                                          ccccGATGAACATCTGCGTTTTGCCTGTCCGCGACAGCCGCGTCTGAACATCAATGTATGTCCCATCGTCGAGTGTCTCGTCGTGAGTTAGGTGATGAAGCGTAGGGTCAGACAAAGCCCACAAAACATCACCGCGAATTCTCATgag              >  4:289354/1‑148 (MQ=255)
                                                                                                                                           cccGATGAACTTCTGCGTTTTGCCTGTCCGCGACAGCCGCGTCTGAACATCAATGTATGTCCCATCGTCGAGTGTCTCGTCGTGAGTTCGGTGATGAAGCGTCGGGTCAGCCCAGGCCCAGAAAACATCACCGCGt                         <  5:357689/136‑1 (MQ=255)
                                                                                                                                           cccGATGAACTTCTGCGTTTTGCCTGTCCGCGACAGCCGCGTCTGAACATCAATGTATGTCCCATCGTCGAGTGTCTCGTCGTGAGTTCGGTGATGAAGCGTCGGGTCAGCCCAGGCCCAGAAAACATCACCGCGt                         >  6:357689/1‑136 (MQ=255)
                                                                                                                                                   aCATCTGCGTTTTGCCTGTCCGCGACAGCCGCGTCTGAACATCAATGTATGTCCCATCGTCGAGTGTCTCGTCGTGAGTTCGGTGATGAAGCGTCGGGTCAGCCCAGGCCCAGAAAACATCACCGCGTATTCTCATGAGACGCTCCTAc    >  4:352619/1‑149 (MQ=255)
                                                                                                                                                   aCATCTGCGTTTTGCCTGTCCGCGACAGCCGCGTCTGAACATCAATGTATGTCCCATCGTCGAGTGTCTCGTCGTGAGTTCGGTGATGAAGCGTCGGGTCAGCCCAGGCCCAGAAAACATCACCGCGTATTCTCATGAGACGCTCCTAc    <  5:212697/149‑1 (MQ=255)
                                                                                                                                                     aTCTGCGTTTTGCCTGTCCGCGACAGCCGCGTCTGAACATCAATGTATGTCCCATCGTCGAGTGTCTCGTCGTGAGTTCGGTGATGAAGCGTCGGGTCAGCCCAGGCCCACAAAAAATGACCGCGTATTCTCATGAGACGCTACTACGa  >  2:205353/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AAAGCCTGCCGCAGGCGCACCTTCAGCGGCGACACGTCGGGCCCGGCCAACGCCCCACGCCAGCGCCCTGGTCATCGACTCGCCAGGGCGAGAATCGAAGGCCTCCTCGAGGA‑TAGCCATGCCACCAGGGGCGTACACCCCGATGAACATCTGCGTTTTGCCTGTCCGCGACAGCCGCGTCTGAACATCAATGTATGTCCCATCGTCGAGTGTCTCGTCGTGAGTTCGGTGATGAAGCGTCGGGTCAGCCCAGGCCCAGAAAACATCACCGCGTATTCTCATGAGACGCTCCTACGA  >  NC_002947/4366453‑4366749

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