Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,125,641 A→T 9.1% L64H (CTC→CAC)  PP_3629 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,125,6410AT 9.1% 118.4 / 3.7 44L64H (CTC→CAC) PP_3629hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (19/21);  new base T (2/2);  total (21/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.60e-01

CCCCATTTGTCCCTGCACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGG  >  NC_002947/4125494‑4125788
                                                                                                                                                   |                                                                                                                                                   
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)
                                                                                                                                                  gAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCg   <  1:184289/148‑1 (MQ=255)
                                                                                                                                                   aGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCgg  >  2:265919/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CCCCATTTGTCCCTGCACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGG  >  NC_002947/4125494‑4125788

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