Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,125,639 T→G 8.8% N65H (AAT→CAT)  PP_3629 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,125,6390TG 8.8% 121.9 / 2.5 43N65H (AAT→CAT) PP_3629hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (18/20);  new base G (2/2);  total (21/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.00e+00

CCCCATTTGTCCCTGCACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCT  >  NC_002947/4125494‑4125784
                                                                                                                                                 |                                                                                                                                                 
ccccATTTGTCCCTGCACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCCACAAGCGGTTGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAg                                                                                                                                                <  6:121239/149‑1 (MQ=255)
ccccATATGTCCCTGCACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAg                                                                                                                                                <  4:215091/149‑1 (MQ=255)
         tCCCTGCACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGc                                                                                                                                        >  7:83615/1‑148 (MQ=255)
          cccTGCACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGc                                                                                                                                        >  2:184289/1‑147 (MQ=255)
           ccTGCACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAg                                                                                                                                     <  8:38363/149‑1 (MQ=255)
                       gCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTa                                                                                                                         >  6:249354/1‑149 (MQ=255)
                         aaaCGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTAcc                                                                                                                       <  5:249354/149‑1 (MQ=255)
                         aaaCGCTCCCCTCGCCGTCGTTAGCGACAATGGCATCCCACAAGCGGTCGGGCGCTCCCAGGTCCGCGGAAGTCACGGGCAACGAAACAGCCTCACACTGCGTGAAGCCCGGGCCCCCATAGAGGCGCCCGCCAGGGATACCCCTTagc                                                                                                                       >  6:148714/1‑147 (MQ=255)
                                         gtcgtTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTAACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCTCAATGGTGGcccg                                                                                                                                           <  6:188786/113‑2 (MQ=255)
                                         gtcgtTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTAACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCTCAATGGTGGcccg                                                                                                                                           >  5:188786/1‑112 (MQ=255)
                                                   aCAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAACTGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGcc                                                                                             <  6:59421/149‑1 (MQ=255)
                                                   aCAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAg                                                                                                                                     <  1:29134/109‑1 (MQ=255)
                                                   aCAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAg                                                                                                                                     >  2:29134/1‑109 (MQ=255)
                                                           aTCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGc                                                                                     <  5:180595/149‑1 (MQ=255)
                                                                       gtcggtgTCTGCCTGTTCGGCGGTACTTAACTGGAACGAAAATGCCTCACAGTGCTTGATGGGCGGGCCACCATTGAGTCGCACGCAAGGGATACCCATTAGTGTGAAGTCCACGGTGATGACATCGCCCTGCTTGCCTGAAGGGTAAt                                                                         <  6:205514/149‑1 (MQ=255)
                                                                       gtcggtcTCTGCCTGTGCCTCGGTACTGACCTTGAACGAACATGCATAACAGTGCTTGAAGGCCGCTCCAACATTGATGCCCACGCAAGGGATACCCATTACCGTGAATTCCAAGGTTATCACATCGCCCTTCTTGCCTGAAGGGTAAt                                                                         <  6:201875/149‑1 (MQ=255)
                                                                              tctGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATcggcgg                                                                   >  7:205245/1‑148 (MQ=255)
                                                                                tGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATcggcgg                                                                   <  8:83615/146‑1 (MQ=255)
                                                                                      gTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGGAATTCACGGTAAtg                                                                                                     >  7:99903/1‑105 (MQ=255)
                                                                                      gTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGGAATTCACGGTAAtg                                                                                                     <  8:99903/106‑2 (MQ=255)
                                                                                                 tGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGAc                                                >  1:121849/1‑148 (MQ=255)
                                                                                                 tGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGAc                                                >  1:117792/1‑148 (MQ=255)
                                                                                                     cTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAAccc                                           >  3:23591/1‑149 (MQ=255)
                                                                                                      tGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATcggcg                                                                    <  4:190161/123‑1 (MQ=255)
                                                                                                      tGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATcggcg                                                                    >  3:190161/1‑123 (MQ=255)
                                                                                                         aaCGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGc                                        >  1:98988/1‑148 (MQ=255)
                                                                                                           cGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTAt                                     >  3:207603/1‑149 (MQ=255)
                                                                                                                   ccTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGt                                                          <  2:204650/120‑1 (MQ=255)
                                                                                                                       acaGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGaa                                                                                <  2:257582/94‑1 (MQ=255)
                                                                                                                       acaGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGaa                                                                                >  1:257582/1‑94 (MQ=255)
                                                                                                                       acaGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGaa                                                                                >  3:11775/1‑94 (MQ=255)
                                                                                                                       acaGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGaa                                                                                <  4:11775/94‑1 (MQ=255)
                                                                                                                       acaGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGt                           <  2:121849/147‑1 (MQ=255)
                                                                                                                       acaGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGt                           <  2:117792/147‑1 (MQ=255)
                                                                                                                         aGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGTTCgg                       >  4:128434/1‑149 (MQ=255)
                                                                                                                         aGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGTTCgg                       <  4:207603/149‑1 (MQ=255)
                                                                                                                                gAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTa                  <  8:205245/147‑1 (MQ=255)
                                                                                                                                  tgcgtgggCCTCAATGGTGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCgg                                         >  3:109971/6‑122 (MQ=255)
                                                                                                                                  tgcgtgggCCTCAATGGTGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCgg                                         <  4:109971/117‑1 (MQ=255)
                                                                                                                                    gCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGTTCgg                       >  7:242895/1‑138 (MQ=255)
                                                                                                                                    gCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGTTCgg                       <  8:242895/138‑1 (MQ=255)
                                                                                                                                      cGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTcg          >  4:150482/1‑149 (MQ=255)
                                                                                                                                              caTTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCctct  >  1:35763/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
CCCCATTTGTCCCTGCACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCT  >  NC_002947/4125494‑4125784

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