Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,492,111 A→G 11.3% L103P (CTG→CCG)  PP_3096 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,492,1110AG11.3% 85.9 / 2.9 37L103P (CTG→CCG) PP_3096hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (21/10);  new base G (2/2);  total (23/14)
Fisher's exact test for biased strand distribution p-value = 5.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

TGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGA‑‑TACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGCCTTCGCGCTCGAAGTGCAGGTCGCTGACCTCGCTGGCCGGGAACGCCAGCGACAGCTGGGTACGAAAA  >  NC_002947/3491974‑3492257
                                                                                                                                           |                                                                                                                                                  
tgtgGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGA‑‑TACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGAc                                                                                                                                         >  3:277277/1‑149 (MQ=255)
  tgGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGA‑‑TACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGCTGCCACGCCCgtaca                                                                                                                                       <  5:218096/149‑5 (MQ=255)
          gTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGA‑‑TACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGAc                                                                                                                                 >  6:315095/1‑147 (MQ=255)
             tttttGCCAGGCCTCTTATAACAAGGTGGTAATGGGGTGCTATAAGATGTCGAGAAACGCTTGGGCGGCATCGTCGCGA‑‑TACTGCACGGGGCGCTCGAGCAGCAGCTCGGTTTAGGGGCGTGGCAGCACGCCCGACGGCCCGACCAgg                                                                                                                             <  6:271583/148‑1 (MQ=255)
                      ggcttggtagtaCATCGTGGTCATGCGGTGCTTGAATATGTCTTTAAACGCTTGGGCGGCATCGTCGCGA‑‑TACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGGGTGGCTGCACGCCCGACGGCCCGACCAGGCCCATGAAg                                                                                                                    <  5:365571/142‑1 (MQ=255)
                          tCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGA‑‑TACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGAc                                                                                                               >  1:313849/1‑149 (MQ=255)
                            gTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGAGGTACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGAc                                                                                                               >  1:302264/1‑149 (MQ=255)
                              aGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGA‑‑TACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGc                                                                                                           >  6:109612/1‑149 (MQ=255)
                                      gTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGA‑‑TACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGa                                                                                                   >  1:256628/1‑149 (MQ=255)
                                         gTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGA‑‑TACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAg                                                                                                >  7:289268/1‑149 (MQ=255)
                                                 gTGCGAGAAGATGGCGAGAAACGCATGGGCGGCATGGTCGCGA‑‑TACTGCACGTGGCGCTCGAGCAGCAGCGCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAgc                                                                                        <  1:87523/149‑1 (MQ=255)
                                                  tGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGA‑‑TACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAgcc                                                                                       >  6:1672/1‑149 (MQ=255)
                                                           aTGTCGAGAAACGCATGGGCGGCATCGTCGCGA‑‑TACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGa                                                                              >  8:153604/1‑149 (MQ=255)
                                                             gTCGAGAAACGCATGGGCGGCATCGTCGCGA‑‑TACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGCTGCCACGCCCg                                                                                                                                           >  6:177967/1‑86 (MQ=255)
                                                             gTCGAGAAACGCATGGGCGGCATCGTCGCGA‑‑TACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGCTGCCACGCCCg                                                                                                                                           <  5:177967/86‑1 (MQ=255)
                                                                         aTGGGCGGCATCGTCGCGA‑‑TACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGCCTTcgc                                                                >  2:1862/1‑149 (MQ=255)
                                                                             gcggcATCGTCGCGA‑‑TACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGCCTTCGCGCTc                                                            >  3:38312/1‑149 (MQ=255)
                                                                                     gtcgCGA‑‑TACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGCCTTCGCGCTCGAAGTg                                                      >  3:149385/1‑147 (MQ=255)
                                                                                     ggcgctA‑‑TACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGCCTTCGCGCTCGAAGTGCa                                                    <  5:118697/147‑1 (MQ=255)
                                                                                        gcgA‑‑TACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGCCTTCGCGCTCGAAGTGCAGGt                                                 <  2:256628/149‑1 (MQ=255)
                                                                                        gcgA‑‑TACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGCCTTCGCGCTCGAAGTGCAGGt                                                 >  8:160150/1‑149 (MQ=255)
                                                                                                 tGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGCCTTCGCGCTCGAAGTGCAGGTCGCTGAc                                          >  3:104924/1‑149 (MQ=255)
                                                                                                    aCGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGCCTTCGCGCTCGAAGTGCAGGTCGCTGACCTc                                       >  2:314390/1‑149 (MQ=255)
                                                                                                      gTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGCCTTCGCGCTCGAAGTGCAGGTCGCTGACCTCGc                                     <  1:1862/149‑1 (MQ=255)
                                                                                                        gtacaccGAGCTGCTGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGCCTTCGCGCTCGAAGTGCAGGTCGCTGACCTCGCTg                                   >  7:308399/7‑149 (MQ=255)
                                                                                                            accGAGCTGCTGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGCCTTCGCGCTCGAAGTGCAGGTCGCTGACCTCGCTGGCCg                               >  4:133864/3‑149 (MQ=255)
                                                                                                            accGAGCTGCTGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGCCTTCGCGCTCGAAGTGCAGGTCGCTGACCTCGCTGGCCg                               >  1:242926/3‑149 (MQ=255)
                                                                                                               gAGCCGCAGCGCGGTGACCGGGCGGGGCCGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGCCTTCGCGCTCGAAGTGCAGGTCGCTGACCTCGCTGGCCGGGa                            <  8:67520/149‑1 (MQ=255)
                                                                                                               gAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGCCTTCGCGCTCGAAGTGCAGGTCGCTGACCTCGCTGGCCggg                             >  6:363341/1‑148 (MQ=255)
                                                                                                                       tgggtgtgtACGTGGGTGGAAGCCTGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACTCCTTGCCTTCGCGCTCGAAGTGCAGGTCGCTGACCTCGCTGGCCGGGAACGCCAGc                    <  6:262005/144‑1 (MQ=255)
                                                                                                                        cTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGCCTTCGCGCTCGAAGTGCAGGTCGCTGACCTCGCTGGCCGGGAACGCCAGCg                   >  1:360902/1‑149 (MQ=255)
                                                                                                                              gtACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGCCTTCGCGCTCGAAGTGCAGGTCGCTGACCTCGCTGGCCGGGAACGCCAGCGACAGCt             <  8:289268/149‑1 (MQ=255)
                                                                                                                              gtACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGCCTTCGCGCTCGAAGTGCAGGTCGCTGACCTCGCTGGCCGGGAACGCCAGCGACAGCt             <  2:302264/149‑1 (MQ=255)
                                                                                                                                    gCGTGGCGGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGCCTTCGCGCTCGAAGTGCAGGTCGCTGACCTCGCTGGCCGGGAACGCCAGCGACAGCTGGGTAc       >  5:150623/1‑149 (MQ=255)
                                                                                                                                    gCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGCCTTCGCGCTCGAAGTGCAGGTCGCTGACCTCGCTGGCCGGGAACGCCAGCGACAGCTGGGTAc       >  1:203917/1‑149 (MQ=255)
                                                                                                                                    gCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGCCTTCGCGCTCGAAGTGCAGGTCGCTGACCTCGCTGGCCGGGAACGCCAGCGACAGCTGGGTAc       <  2:203917/149‑1 (MQ=255)
                                                                                                                                       gggCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGCCTTCGCGCTCGAAGTGCAGGTCGCTGACCTCGCTGGCCGGGAACGCCAGCGACAGCTGGGTACGaa    <  4:104924/148‑1 (MQ=255)
                                                                                                                                         gcagcaCGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGCCTTCGCGCTCGAAGTGCAGGTCGCTGACCTCGCTGGCCGGGAACGCCAGCGACAGCTGGGTACGaaaa  <  5:109612/149‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                  
TGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGA‑‑TACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGCCTTCGCGCTCGAAGTGCAGGTCGCTGACCTCGCTGGCCGGGAACGCCAGCGACAGCTGGGTACGAAAA  >  NC_002947/3491974‑3492257

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