Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,623,786 C→A 7.8% L89L (CTG→CTT)  PP_3194 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,623,7860CA 7.8% 158.7 / 3.5 51L89L (CTG→CTT) PP_3194hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (24/23);  new base A (2/2);  total (26/25)
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.08e-01

CGGCTGCTTGCAGCGATACGGGGGCCTTCGAACAGTCCGGCGTTGACCGTATTGCGCACGGCAATGCTGAGTTGGCCACTGTCGCCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGAACAGCCCCCGCAGGGCGCGAAACTCGTCGCTGGTCCAGATGTCGATGTCTTCCTCGCTCTGGGCATTGCCAAACGCCAGGTGGGTGTGCACGTCCACCAGCCCGGGCATGACGAACTGGCCGTGGGCTTCACGGTCGCCGGGTTGTGGACGGGGTG  >  NC_002947/3623640‑3623931
                                                                                                                                                  |                                                                                                                                                 
cGGCTGCTTGCAGCGATACGGGGGCCTTCGAACAGTCCGGCGTTGACCGTATTGCGCACGGCAATGCTGAGTTGGCCACTGTCGCCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGAACAGCCCCCGCAg                                                                                                                                                 >  5:306247/1‑149 (MQ=255)
 ggCTGCTTGCAGCGATACGGGGGCCTTCGAACAGTCCGGCGTTGACCGTATTGCGCACGGCAATGCTGAGTTGGCCACTGTCGCCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGAACAGCCCCCGCAgg                                                                                                                                                <  2:155551/149‑1 (MQ=255)
          cAGCGATACGGGGTCCTTCGAACAGTCCGGCGTTGACCGTATTGCGCACGGCAATGCTGAGTTGGCCACTGTCGCCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGAACAGCCCCCGCAGGGCGCGAAAc                                                                                                                                       <  5:357964/149‑1 (MQ=255)
           aGCGATACGGGGGCCTTCGAACAGTCCGGCGTTGACCGTATTGCGCACGGCAATGCTGAGTTGGCCACTGTCGCCGGGGCAAACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGAACAGCCCCCGCAGGGCGCGAAACt                                                                                                                                      >  3:95951/1‑149 (MQ=255)
            gCGATACGGGGGCCTTCGAACAGTCCGGCGTTGACCGTATTGCGCACGGCAATGCTGAGTTGGCCACTGTCGCCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGAACAGCCCCCGCAGGGCGCGAAACtc                                                                                                                                     >  8:201593/1‑149 (MQ=255)
                   gggggCCTTCGAACAGTCCGGCGTTGACCGTATTGCGCACGGCAATGCTGAGTTGGCCACTGTCGCCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGAACAGCCCCCGCAGGGCGCGAAACTCGTCGCTg                                                                                                                              <  7:201593/149‑1 (MQ=255)
                        ccTTCGAACAGTCCGGCGTTGACCGTATTGCGCACGGCAATGCTGAGTTGGCCACTGTCGCCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGAACAGCCCCCGCAGGGCGCGAAACTCGTCGCTGGTCCa                                                                                                                         <  7:272106/149‑1 (MQ=255)
                         cTTCGAACAGTCCGGCGTTGACCGTATTGCGCACGGCAATGCTGAGTTGGCCACTGTCGCCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGAACAGCCCCCGCAGGGCGCGAAACTCGTCGCTGGTc                                                                                                                           >  4:134120/1‑146 (MQ=255)
                             gAACAGTCCGGCGTTGACCGTATTGCGCACGGCAATGCTGAGTTGGCCACTGTCGCCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGAACAGCCCCCGCAGGGCGCGAAACTCGTCGCTGGTCCAgatgt                                                                                                                    >  2:185935/1‑149 (MQ=255)
                              aaCAGTCCGGCGTTGACCGTATTGCGCACGGCAATGCTGAGTTGGCCACTGTCGCCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGAACAGCCCCCGCAGGGCGCGAAACTCGTCGCTGGTCCAgatgtc                                                                                                                   >  2:88837/1‑149 (MQ=255)
                                 aGTCCGGCGTTGACCGTATTGCGCACGGCAATGCTGAGTTGGCCACTGTCGCCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGAACAGCCCCCGCAGGGCGCGAAACTCGTCGCTGGTCCAgatgtcgat                                                                                                                <  6:306247/149‑1 (MQ=255)
                                         gTTGACCGTATTGCGCACGGCAATGCTGAGTTGGCCACTGTCGCCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGAACAGCCCCCGAAGGGCGCGAAACTCGTCGCTGGTCCAgatgtcg                                                                                                                  >  4:150724/1‑139 (MQ=255)
                                         gTTGACCGTATTGCGCACGGCAATGCTGAGTTGGCCACTGTCGCCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGAACAGCCCCCGAAGGGCGCGAAACTCGTCGCTGGTCCAgatgtcg                                                                                                                  <  3:150724/139‑1 (MQ=255)
                                                 tATTGCGCACGGCAATGCTGAGTTGGCCACTGTCGCCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGAACAGCCCCCGCAGGGCGCGAAACTCGTCGCTGGTCCAGATGTCGATGTCTTCCTCGCTCTgg                                                                                                >  1:154571/1‑149 (MQ=255)
                                                 tATTGCGCACGGCAATGCTGAGTTGGCCACTGTCGCCGGGGCAAACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGAACAGCCCCCGCAGGGCGCGAAACTCGTCGCTGGTCCAGATGTCGATGTCTTCCTCGCTCTgg                                                                                                <  4:95951/149‑1 (MQ=255)
                                                            gCAATGCTGAGTTGGCCACTGTCGCCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGAACAGCCCCCGCAGGGCGCGAAACTCGTCGCTGGTCCAGATGTCGATGTCTTCCTCGCTCTGGGCATTGCCaaa                                                                                     <  3:134120/149‑1 (MQ=255)
                                                            gCAATGCTGAGTTGGCCACTGTCGCCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGAACAGCCCCCGCAGGGCGCGAAACTCGTCGCTGGTCCAGATGTCGATGTCTTCCTCGCTCTGGGCATTGCCaaa                                                                                     <  1:185935/149‑1 (MQ=255)
                                                             cAATGCTGAGTTGGCCACTGTCGCCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGAACAGCCCCCGCAGGGCGCGAAACTCGTCGCTGGTCCAGATGTCGATGTCTTCCTCGCTCTGGGCATTGCCAAAc                                                                                    >  5:334415/1‑149 (MQ=255)
                                                                    gAGTTGGCCACTGTCGCCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGAACAGCCCCCGCAGGGCGCGAAACTCGTCGCTGGTCCAGATGTCGATGTCTTCCTCGCTCTGGGCATTGCCAAACGCCAggt                                                                             >  1:38904/1‑149 (MQ=255)
                                                                         ggCCACTGTCGCCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGAACAGCCCCCGCAGGGCGCGAAACTCGTCGCTGGTCCAGATGTCGATGTCTTCCTCGCTCTGGGCATTGCCAAACGCCAggtgggt                                                                         <  2:38904/148‑1 (MQ=255)
                                                                             aCTGTCGCCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGAACAGCCCCCGCAGGGCGCGAAACTCGTCGCTGGTCCAGATGTCGATGTGTTCCTCGCTCTGGGCATTGCCAAACGCCAggtgg                                                                           >  5:304844/1‑142 (MQ=255)
                                                                             aCTGTCGCCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGAACAGCCCCCGCAGGGCGCGAAACTCGTCGCTGGTCCAGATGTCGATGTCTTCCTCGCTCTGGGCATTGCCAAACGCCAggtgg                                                                           <  6:304844/142‑1 (MQ=255)
                                                                             aCTGTCGCCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGAACAGCCCCCGAAgcgc                                                                                                                                              >  2:259632/1‑75 (MQ=255)
                                                                             aCTGTCGCCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGAACAGCCCCCGAAGGgc                                                                                                                                              <  1:259632/75‑1 (MQ=255)
                                                                                           aGACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGAACAGCCCCCGCAGGGCGCGAAACTCGTCGCTGGTACAGATGTCGATGTCTTCCTCGCTCTGGGCATTGCCAAACGCCAGGTGGGTGTGAACGTCCACCAGCCCg                                                      >  2:82371/1‑149 (MQ=255)
                                                                                             aCCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGAACAGCCCCCGCAGGGCGCGAAACTCGTCGCTGGTCCAGATGTCGATGTCTTCCTCGCTCTGGGCATTGCCAAACGCCAGGTGGGTGTGCACGTCCACCAGCCCggg                                                    >  6:275865/1‑149 (MQ=255)
                                                                                                   gACGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGAACAGCCCCCGCAGGGCGCGAAACTCGTCGCTGGTCCAGATGTCGATGTCTTCCTCGCTCTGGGCATTGCCAAACGCCAGGTGGGTGTGCACGTCCACCAGCCCGGGCATGa                                               <  5:275865/148‑1 (MQ=255)
                                                                                                                   ccAGAACGTGCTGAGCGAAGAACAGCCCCCGCAGGGCGCGAAACTCGTCGCTGGTCCAGATGTCGATGTCTTCCTCGCTCTGGGCATTGCCAAACGCCAggtgg                                                                           <  4:169456/104‑1 (MQ=255)
                                                                                                                   ccAGAACGTGCTGAGCGAAGAACAGCCCCCGCAGGGCGCGAAACTCGTCGCTGGTCCAGATGTCGATGTCTTCCTCGCTCTGGGCATTGCCAAACGCCAggtgg                                                                           >  3:169456/1‑104 (MQ=255)
                                                                                                                     aGAACGTGCTGAGCGAAGAACAGCCCCCGCAGGGCGCGAAACTCGTCGCTGGTCCAGATGTCGATGTCTTCCTCGCTCTGGGCATTGCCAAACGCCAGGTGGGTGTGCACGTc                                                                >  4:263995/1‑113 (MQ=255)
                                                                                                                     aGAACGTGCTGAGCGAAGAACAGCCCCCGCAGGGCGCGAAACTCGTCGCTGGTCCAGATGTCGATGTCTTCCTCGCTCTGGGCATTGCCAAACGCCAGGTGGGTGTGCACGTc                                                                <  3:263995/113‑1 (MQ=255)
                                                                                                                            gCTGAGCGAAGAACAGCCCCCGCAGGGCGCGAAACTCGTCGCTGGTCCAGATGTCGATGTCTTCCTCGCTCTGGGCATTGCCAAACGCCAGGTGGGTGTGCACGTCCACCAGCCCGGGCATGACGAACTGGcc                                     <  4:185518/133‑1 (MQ=255)
                                                                                                                            gCTGAGCGAAGAACAGCCCCCGCAGGGCGCGAAACTCGTCGCTGGTCCAGATGTCGATGTCTTCCTCGCTCTGGGCATTGCCAAACGCCAGGTGGGTGTGCACGTCCACCAGCCCGGGCATGACGAACTGGcc                                     >  3:185518/1‑133 (MQ=255)
                                                                                                                                aGCGAAGAACAGCCCCCGCAGGGCGCGAAACTCGTCGCTGGTCCAGATGTCGATGTCTTCCTCGCTCTGGGCATTGCCAAACGCCAGGTGGGTGTGc                                                                     <  1:113467/97‑1 (MQ=255)
                                                                                                                                aGCGAAGAACAGCCCCCGCAGGGCGCGAAACTCGTCGCTGGTCCAGATGTCGATGTCTTCCTCGCTCTGGGCATTGCCAAACGCCAGGTGGGTGTGc                                                                     >  2:113467/1‑97 (MQ=255)
                                                                                                                                  cGAAGAACAGCCCCCGCAGGGCGCAAAACTCGTCGCTGGTCCAGATGTCGATGTCTTCCTCGCTCTGGGCATTGCCAAACGCCAGGTGGGTGTGCACGTCCACCAGCCCGGGCAt                                                 >  6:257217/1‑115 (MQ=255)
                                                                                                                                  cGAAGAACAGCCCCCGCAGGGCGCAAAACTCGTCGCTGGTCCAGATGTCGATGTCTTCCTCGCTCTGGGCATTGCCAAACGCCAGGTGGGTGTGCACGTCCACCAGCCCGGGCAt                                                 <  5:257217/115‑1 (MQ=255)
                                                                                                                                   gaagaaCAGCCCCCGCAGGGCGCGAAACTCGTCGCTGGTCCAGATGTCGATGTCTTCCTCGCTCTGGGCATTGCCAAACGCCAGGTGGGTGTGCGCGTCCACCAGCCCgg                                                     <  5:9967/110‑1 (MQ=255)
                                                                                                                                   gaagaaCAGCCCCCGCAGGGCGCGAAACTCGTCGCAGGTCCAGATGTCGATGTCTTCCTCGCTCTGGGCATTGCCAAACGCCAGGTGGGTGTGCGCGTCCACCAGCCCgg                                                     >  6:9967/1‑110 (MQ=255)
                                                                                                                                     agaaCAGCCCCCGCAGGGCGCGAAACTCGTCGCTGGTCCAGATGTCGATGTCTTCCTCGCTCTGGGCATTGCCAAACGCCAGGTGGGTGTGCACGTCCACCAGCCCggg                                                    <  6:25132/109‑1 (MQ=255)
                                                                                                                                     agaaCAGCCCCCGCAGGGCGCGAAACTCGTCGCTGGTCCAGATGTCGATGTCTTCCTCGCTCTGGGCATTGCCAAACGCCAGGTGGGTGTGCACGTCCACCAGCCCggg                                                    >  5:25132/1‑109 (MQ=255)
                                                                                                                                     agaaCAGCCCCCGCAGGGCGCGAAACTCGTCGCTGGTCCAGATGTCGATGTCTTCCTCGCTCTGGGCATTGCCAAACGCCAGGTGGGTGTGCACGTCCACCAGCCCGGGCATGACGAACTGGCCGTGGGCTTCACGGTCGCCGGGTtgt            <  1:88837/149‑1 (MQ=255)
                                                                                                                                              cccGCAGGGCGCGAAACTCGTCGCTGGTCCAGATGTCGATGTCTTCCTCGCTCTGGGCATTGCCAAACGCCAGGTGGGTGTGCACGTCCACCAGCCCGGGc                                                   <  4:143314/101‑1 (MQ=255)
                                                                                                                                              cccGCAGGGCGCGAAACTCGTCGCTGGTCCAGATGTCGATGTCTTCCTCGCTCTGGGCATTGCCAAACGCCAGGTGGGTGTGCACGTCCACCAGCCCGGGc                                                   >  3:143314/1‑101 (MQ=255)
                                                                                                                                              cccGCAGGGCGCGAAACTCGTCGCTGGTCCAGATGTCGATGTCTTCCTCGCTCTGGGCATTGCCAAACGCCAGGTGGGTGTGCACGTCCACCAGCCCGGGCATGACGAACTGGCCGTGGGCTTCACGGTCGCCGGGTTGTGGACGGGGt   <  3:236061/149‑1 (MQ=255)
                                                                                                                                               ccGCAGGGCGCGAAACTCGTCGCTGGTCCAGATGTCGATGTCTTCCTCGCTCTGGGCATTGCCAAACGCCAGGTGGGTGTGCACGTCCACCAGCCCGGGCATGACGAACTg                                        <  7:56011/111‑1 (MQ=255)
                                                                                                                                               ccGCAGGGCGCGAAACTCGTCGCTGGTCCAGATGTCGATGTCTTCCTCGCTCTGGGCATTGCCAAACGCCAGGTGGGTGTGCACGTCCACCAGCCCGGGCATGACGAACTg                                        >  8:56011/1‑111 (MQ=255)
                                                                                                                                               ccGCAGGGCGCGAAACTCGTCGCTGGTCCAGATGTCGATGTCTTCCTCGCTCTGGGCATTGCCAAACGCCAGGTGGGTGTGCACGTCCACCAGCCCGGGCATGACGAACTGGCCGTGGGCTTCACGGTCGCCGGGTTGTGGACGGGGTg  >  3:77181/1‑149 (MQ=255)
                                                                                                                                               ccGCAGGGCGCGAAACTCGTCGCTGGTCCAGACGTCGATGTCTTCCTCGCTCTGGGCATTGCCAAACGCCAGGTGGGTGTGCACGTCCACCAGCCCGGGCATGACGAACTGGCCGTGGGCTTCACGGTCGCCGGGTtg             >  6:134359/1‑138 (MQ=255)
                                                                                                                                                 gCAGGGCTCGAAACTCGTGGCTGGTCCAGATGTCGTTGTCTTCCTCGCTCTGGGCATTGCCAAACGCCAGGTGGGTGTGCACGTCCACCAGCCCGGGCATGACGAACTGGCCGTGGGCTTCACgg                        <  2:259042/125‑1 (MQ=255)
                                                                                                                                                 gCAGGGCGCGAAACTCGTCGCTGGTCCAGATGTCGATGTCATCCTCGCTCTGGGCATTGCCAAACGCCAGGTGGGTGTGCACGTCCACCAGCCCGGGCATGACGAACTTGCCGTGGGCTTCACgc                        >  1:259042/1‑124 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CGGCTGCTTGCAGCGATACGGGGGCCTTCGAACAGTCCGGCGTTGACCGTATTGCGCACGGCAATGCTGAGTTGGCCACTGTCGCCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGAACAGCCCCCGCAGGGCGCGAAACTCGTCGCTGGTCCAGATGTCGATGTCTTCCTCGCTCTGGGCATTGCCAAACGCCAGGTGGGTGTGCACGTCCACCAGCCCGGGCATGACGAACTGGCCGTGGGCTTCACGGTCGCCGGGTTGTGGACGGGGTG  >  NC_002947/3623640‑3623931

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