Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,181,401 A→G 9.7% E422E (GAA→GAG)  PP_3678 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,181,4010AG 9.7% 105.1 / 2.8 41E422E (GAA→GAG) PP_3678hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (21/16);  new base G (2/2);  total (23/18)
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.76e-01

AGAATTCGCCTCGAGCTTAGGTTCGCTAAGCATCCGGAGTTAAAGAATGGATGGACTGCCGGGCTTCGAGATGAGGTGTGTCGAGAAGTTATAGAACTGACTGGCCACGCAGATCCAGCTTCACTGCATCATTACTTTCACAGTGAATATACAGCCGTCACTTCAGATGCTACCTACAGCTCTATGCTTAAACTGAAGGATGAAATCGATGCTGCAAATGAGACTATCACAGCCTTAGAGCACAAAGCGAAGATTGGCAAAATAGACTGTTCAGTAGAGATTACGTCACTTAAGT  >  NC_002947/4181255‑4181549
                                                                                                                                                  |                                                                                                                                                    
aGAATTCGCCTCGAGCTTAGGTTCGCTAAGCATCCGGAGTTAAAGAATGGATGGACTGCCGGGCTTCGAGATGAGGTGTGTCGAGAAGTTATAGAACTGACTGGCCACGCAGATCCAGCTTCACTGCATCATTACTTTCACAGTGAata                                                                                                                                                    >  3:292102/1‑149 (MQ=255)
       gCCTCGAGCTTAGGTTCGCTAAGCATCCGGAGTTAAAGAATGGATGGACTGCCGGGCTTCGAGATGAGGGGTGTCGAGAAGTTATAGAACTGACTGGCCACGCAGATCCAGCTTCACTGCATCATTACTTTCACAGTGAATATACAGc                                                                                                                                              <  2:111033/148‑1 (MQ=255)
                        gCTAAGCATCCGGAGTTAAAGAATGGATGGACTGCCGGGCTTCGAGATGAGGTGTGTCGAGAAGTTATAGAACTGACTGGCCACGCAGATCCAGCTTTACTGCATCATTACTTTCACAGCGAATATACAGCCGTCACTTCAGATGcta                                                                                                                             <  6:136897/148‑1 (MQ=255)
                        gCTAAGCATCCGGAGTTAAAGAATGGATGGACTGCCGGGCTTCGAGATGACGTGTGTCGAGAAGTTATAGAACTGACTGGCCACGCAGATCCAGCTTCACTGCATCATTACTTTCACAGTGAATATACAGCCGTCACTTCAGATGctac                                                                                                                            <  2:41271/149‑1 (MQ=255)
                             acATTCGGAGTTAAAGAATGGATGGACTGCCGGGCTTCGAGATGAGGTGTGTCGAGAAGTTATAGAACTGACTGGCCACGCAGATCCAGCTTCACTGCATCATTACTTTCACAGTGAATATACAGCCGTCACTTCAGATGctaccta                                                                                                                         >  2:142259/2‑147 (MQ=255)
                                         aaaGAATGGATGGACTGCCGGGCTTCGAGATGAGGTGTGTCGAGAAGTTATAGAACTGACTGGCCACGCAGATCCAGCTTCACTGCATCATTACTTTCACAGTGAATATACAGCCGTCACTTCa                                                                                                                                    >  1:66547/1‑124 (MQ=255)
                                         aaaGAATGGATGGACTGCCGGGCTTCGAGATGAGGTGTGTCGAGAAGTTATAGAACTGACTGGCCACGCAGATCCAGCTTCACTGCATCATTACTTTCACAGTGAATATACAGCCGTCACTTCa                                                                                                                                    <  2:66547/124‑1 (MQ=255)
                                         aaaGAATGGATGGACTGCCGGGCTTCGAGATGAGGTGTGTCGAGAAGTTATAGAACTGACTGGCCACGCAGATCCAGCTTCACTGCATCATTACTTTCACAGTGAATATACAGCCGTCACTTCAGATGCTACCTACAGCTCTATGCtt                                                                                                            >  3:75533/1‑148 (MQ=255)
                                           aGAATGGATGGACTGCCGGGCTTCGAGATGAGGTGTGTCGAGAAGTTATAGAACTGACTGGCCACGCAGATCCAGCTTCACTGCATCATTACTTTCACATTGAATATACAGCCGTCACTTCAGATGCTACCTACAGCTCTATGCTTaa                                                                                                          >  5:37807/1‑148 (MQ=255)
                                                        tGCCGGGCTTCGAGATGAGGTGTGTCGAGAAGTTATAGAACTGACTGGCCACGCAGATCCAGCTTCACTGCATCATTACTTTCACAGTGAatat                                                                                                                                                   <  3:49896/94‑1 (MQ=255)
                                                        tGCCGGGCTTCGAGATGAGGTGTGTCGAGAAGTTATAGAACTGACTGGCCACGCAGATCCAGCTTCACTGCATCATTACTTTCACAGTGAatat                                                                                                                                                   >  4:49896/1‑94 (MQ=255)
                                                        tGCCGGGCTTCGAGATGAGGTGTGTCGAGAAGTTATAGAACTGACTGGCCACGCAGATCCAGCTTCACTGCATCATTACTTTCACAGTGAATATACAGCCGTCACTTCAGATGCTACCTACAGCTCTATGCTTAAACTGAAGGATGaaa                                                                                            >  5:181765/1‑149 (MQ=255)
                                                          ccGGGCTTCGAGATGAGGTGTGTCGAGAAGTTATAGAACTGACTGGCCACGCAGATCCAGCTTCACTGCATCATTACTTTCACAGTGAATATACAGCCGTCACTTCAGATGCTACCTACAGCTCTATGCTTAAACTGAAGGa                                                                                                 >  2:303080/1‑142 (MQ=255)
                                                          ccGGGCTTCGAGATGAGGTGTGTCGAGAAGTTATAGAACTGACTGGCCACGCAGATCCAGCTTCACTGCATCATTACTTTCACAGTGAATATACAGCCGTCACTTCAGATGCTACCTACAGCTCTATGCTTAAACTGAAGGa                                                                                                 <  1:303080/142‑1 (MQ=255)
                                                              ncTTCGAGATGAGGTGTGTCGAGAAGTTATAGAACTGACTGGCCACGCAGATCCAGCTTCACTGCATCATTACTTTCACAGTGAGTATACAGCCGTCACTTCAGATGCTACCTACAGCTCTATGCTTAAACTGAAGGATGAAATCGAtg                                                                                      >  4:179886/2‑149 (MQ=255)
                                                              gCTTCGAGATGAGGTGTGTCGAGAAGTTATAGAACTGACTGGCCACGCAGATCCAGCTTCACTGCATCATTACTTTCACAGTGAGTATACAGCCGTCACTTCAGATGCTACCTACAGCTCTATGCTTAAACTGAAGGATGAAATCGAtg                                                                                      >  4:175596/1‑149 (MQ=255)
                                                                ttCGAGATGAGGTGTGTCGAGAAGTTATAGAACTGACTGGCCACGCAGATCCAGCTTCACTGCATCATTACTTTCACAGTGAATATACAGCCGTCACTTCAGATGCTACCTACAGctct                                                                                                                  >  2:208481/1‑119 (MQ=255)
                                                                ttCGAGATGAGGTGTGTCGAGAAGTTATAGAACTGACTGGCCACGCAGATCCAGCTTCACTGCATCATTACTTTCACAGTGAATATACAGCCGTCACTTCAGATGCTACCTACAGCTCTATGCTTAAACTGAAGGATGAAATCGAtgct                                                                                    >  1:283244/1‑149 (MQ=255)
                                                                ttCGAGATGAGGTGTGTCGAGAAGTTATAGAACTGACTGGCCACGCAGATCCAGCTTCACTGCATCATTACTTTCACAGTGAATATACAGCCGGCACTTCAGATGCTACCTACAGCGCTATGCTTAAACTGATGGCTGAAATCGAtgct                                                                                    >  8:231863/1‑149 (MQ=255)
                                                                  cGAGATGAGGTGTGTCGAGAAGTTATAGAACTGACTGGCCACGCAGATCCAGCTTCACTGCATCATTACTTTCACAGTGAATATACAGCCGTCACTTCAGATGCTACCTACAGCTCTATGCTTAAACTGAAGGATGAAATCGAtgctg                                                                                   >  7:99503/1‑148 (MQ=255)
                                                                   gagaTGAGGTGTGTCGAGAAGTTATAGAACTGACTGGCCACGCAGATCCAGCTTCACTGCATCATTACTTTCACAGTGAATATACAGCCGTCACTTCAGATGCTACCTACAGCTCTATGCTTAAACTGAAGGATGAAATCGAtgctgc                                                                                  <  8:99503/148‑1 (MQ=255)
                                                                    agaTGAGGTGTGTCGAGAAGTTATAGAACTGACTGGCCACGCAGATCCAGCTTCACTGCATCATTACTTTCACAGTGAATATACAGCCGTCACTTCAGATGCTACCTACAGCTCTATGCTTAAACTGAAGGATGAAATCGATGCTGCa                                                                                 >  6:315850/1‑148 (MQ=255)
                                                                     gaTGAGGTGTGTCGAGAAGTTATAGAACTGACTGGCCACGCAGATCCAGCTTCACTGCATCATTACTTTCACAGTGAGTATACAGCCGTCACTTCAGATGCTACCTACAGCTCTATGCTTAAACTGAAGGATGAAATCGATGCTGCaaa                                                                               <  3:179886/149‑1 (MQ=255)
                                                                     gaTGAGGTGTGTCGAGAAGTTATAGAACTGACTGGCCACGCAGATCCAGCTTCACTGCATCATTACTTTCACAGTGAGTATACAGCCGTCACTTCAGATGCTACCTACAGCTCTATGCTTAAACTGAAGGATGAAATCGATGCTGCaaa                                                                               <  3:175596/149‑1 (MQ=255)
                                                                     gaTGAGGTGTGTCGAGAAGTTATAGAACTGACTGGCCACGCAGATCCAGCTTCACTGCATCATTACTTTCACAGTGAATATACAGCCGTCACTTCAGATGCTACCTACAGATCTATGCTTAAACTGAAGGATTAAATCGATGCTGCaaa                                                                               >  1:153455/1‑149 (MQ=255)
                                                                       tGAGGTGTGTCGAGAAGTTATAGAACTGACTGGCCACGCAGATCCAGCTTCACTGCATCATTACTTTCACAGTGAATATACAGCCGTCACTTCAGATGCTACCTACAGCTCTATGCTTAAACTGAAGGATGAAATCGATGCTGCAAAt                                                                              <  4:75533/148‑1 (MQ=255)
                                                                         aGGTGTGTCGAGAAGTTATAGAACTGACTGGCCACGCAGATCCAGCTTCACTGCATCATTACTTTCACAGTGAATATACAGCCGTCACTTCAGATGCTACCTACAGCTCTATGCTTAAACTGAAGGATGAAATCGATGCTGCAAATga                                                                            <  2:283244/148‑1 (MQ=255)
                                                                          ggTGTGTCGAGAAGTTATAGAACTGACTGTCCACGCAGATCCAGCTTCAATGCATCATTACTTTAACAGTGAATATACATCCGTCACTTCAGATGCTACCTACAGCTCTATGCTTAAACTGAAGGATGAAATCGATGCTGCAAATgaga                                                                          <  1:288659/149‑1 (MQ=255)
                                                                           gtgtgtCGAGAAGTTATAGAACTGACTGGCCACGCAGATCCAGCTTCACTGCATCATTACTTTCACAGTGAATATACAGCCGTCACTTCAGATGCTACCTACAGCTCTATGCTTAAACTGAAGGATGAAATCGATGCTGCAAATGAGAc                                                                         >  3:136085/1‑149 (MQ=255)
                                                                               gtCGAGAAGTTATAGAACTGACTGGCCACGCAGATCCAGCTTCACTGCATCATTACTTTCACAGTGAATATACAGCCGTCACTTCAGATGCTACCTACAGCTCTATGCTTAAACTGAAGGATGAAATCGATGCTGCAAATGAGACTAt                                                                      <  5:315850/148‑1 (MQ=255)
                                                                                           taGAACTCACTGGCCACGCAGATCCAGCTTCACTGCATCATTACTTTCACAGTGAATATCCAGCCTTCACTTCAGATGCTACCTAGAGCTCTATGCTTAAACTGAAGGATGAAATCGATGCTGCAAATGTGACTATCACAGCCTTagag                                                         <  6:90871/149‑1 (MQ=255)
                                                                                                        ccACGCAGATCCAGCTTCACTGCATCATTACTTTCACAGTGAATATACAGCCGTCACTTCAGATGCTACCTACAGCTCTATGCTTAAACTGAAGGATGAAATCGATGCTGCAAATGAGACTATCACAGCCTTAGAGCACAAAGCg                                                >  6:127829/1‑145 (MQ=255)
                                                                                                         cACGCAGATCCAGCTTCACTGCATCATTACTTTCTCAGTGAATATACAGCCGTCACTTCAGATGCTACCTACAGCTCTATGCTTAAACTGAAGGATGAAATCGATGCTGCAAATGAGACTATCACAGCCTTAGAGCACAAAGCGAAGAt                                           >  1:131881/1‑149 (MQ=255)
                                                                                                          aCGCAGATCCAGCTTCACTGCATCATTACTTTCACAGTGAATATACAGCCGTCACTTCAGATGCTACCTACAGCTCTATGCTTAAACTGAAGGATGAAATCGATGCTGCAAATGAGACTATCACAGCCTTAGAGCACAAAGCGAAGAtt                                          >  5:285152/1‑149 (MQ=255)
                                                                                                          aCGCAGATCCAGCTTCACTGCATCATTACTTTCACAGTGAATATACAGCCGTCACTTCAGATGCTACCTACAGCTCTATGCTTAAACTGAAGGATGAAATCGATGCTGCAAATGAGACTATCACAGCCTTAGAGCACAAAGCGAAGAtt                                          <  4:136085/149‑1 (MQ=255)
                                                                                                             cAGATCCAGCTTCACTGCATCATTACTTTCACAGTGAATATACAGCCGTCACTTCAGATGCTACCTACAGCTCTATGCTTAAACTGAAGGATGAAATCGATGCTGCAAATGAGACTATCACAGCCTTAGAGCACAAAGCGAAGATTGGc                                       <  6:37807/149‑1 (MQ=255)
                                                                                                                     gCTTCACTGCATCATTACTTTCACAGTGAATATACAGCCGTCACTTCAGATGCTACCTACAGCTCTATGCTTAAACTGAAGGATGAAATCGATGCTGCAAATGAGACTATCACAGCCTTAGAGCACAAAGCGAAGATTGGCAAAATAg                                >  8:221497/1‑148 (MQ=255)
                                                                                                                     gCTTCACTGCATCATTACTTTCACAGTGAATATACAGCCGTCACTTCAGATGCTACCTACAGCTCTATGCTTAAACTGAAGGATGAAATCGATGCTGCAAATGAGACTATCACAGCCTTAGAGCACAAAGCGAAGATTGGCAAAATAGa                               <  6:181765/149‑1 (MQ=255)
                                                                                                                              catcatTACTTTCACAGTGAATATACAGCCGTCACTTCAGATGCTACCTACAGCTCTATGCTTAAACTGAAGGATGAAATCGATGCTGCAAATGAGACTATCACAGCCTTAGAGCACAAAGCGAAGATTGGCAAAATAGACTGTTCAg                       >  8:33221/1‑148 (MQ=255)
                                                                                                                                   ttACTTTCACAGTGAATATACAGCCGTCACTTCAGATGCTACCTACAGCTCTATGCTTAAACTGAAGGATGAAATCGATGCTGCAAATGAGACTATCACAGCCTTAGAGCACAAAGCGAAGATTGGCAAAATAGACTGTTCAGTa                     >  3:193496/1‑145 (MQ=255)
                                                                                                                                      cTTTCACAGTGAATATACAGCCGTCACTTCAGATGCTACCTACAGCTCTATGCTTAAACTGAAGGATGAAATCGATGCTGCAATTGAGACTATCACAGACTTAGAGCTCAAAGCGAAGTTTGGCAAAATAGACTGTGCAGTAGAGGTTa              >  8:275409/1‑149 (MQ=255)
                                                                                                                                           acaGTGAATATACAGCCGTCACTTCAGATGCTACCTACAGCTCTATGCTTAAACTGAAGGATGAAATCGATGCTGCAAATGAGACTATCACAGCCTTAGAGCACAAAGCGAAGATTGGCAAAATAGACTGTTCAGTAGAGATTACGTCa         <  6:285152/149‑1 (MQ=255)
                                                                                                                                           acaGTGAATATACAGACGTCACTTCAGATGCTACCTACAGCTCTAAGCTTAATCTTAAGGATCAAATCGATTCTTCAAATGAGACTATCACAGCCTTAGAGCACAAAGCGAAGATTGGCAAAATAGACTGTTCAGTAGAGATTACGTc          <  4:101307/148‑1 (MQ=255)
                                                                                                                                               tGAATATACAGCCGTCACTTCAGATGCTACCTACAGCTCTATGCTTAAACTGAAGGATGAAATCGATGCTGCAAATGAGACTATCACAGCCTTAGAGCACAAAGCGAAGATTGGCAAAATAGACTGTTCAGTAGAGATTACGTCACtt      >  3:279915/1‑148 (MQ=255)
                                                                                                                                                 aaTATACAGCCGTCACTTCAGATGCTACCTACAGCTCTATGCTTAAACTGAAGGATGAAATCGATGCTGCAAATGAGACTATCACAGCCTTAGAGCACAAAGCGAAGATTGGCAAAATAGACTGTTCAGTAGAGATTACGTCACTTaa    <  4:193496/148‑1 (MQ=255)
                                                                                                                                                 aaTATACAGCCGTCACTTCAGATGCTACCTACAGCTCTATGCTTAAACTGAAGGATGAAATCGATGCTGCAAATGAGACTATCACAGCCTTAGAGCACAAAGCGAAGATTGGCAAAATAGACTGTTCAGTAGAGATTACGTCACTTaa    <  7:221497/148‑1 (MQ=255)
                                                                                                                                                  atatACAGCCGTCACTTCAGATGCTACCTACATCTCTATGCTTAAACTGAAGGATGAAATCGATGCTGCAAATGAGACTATCACAGCCTTAGAGCACAAAGCGAAGATTGGCAAAATAGACTGTTCAGTAGAGATTACGTCACTTAAGt  <  5:127829/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
AGAATTCGCCTCGAGCTTAGGTTCGCTAAGCATCCGGAGTTAAAGAATGGATGGACTGCCGGGCTTCGAGATGAGGTGTGTCGAGAAGTTATAGAACTGACTGGCCACGCAGATCCAGCTTCACTGCATCATTACTTTCACAGTGAATATACAGCCGTCACTTCAGATGCTACCTACAGCTCTATGCTTAAACTGAAGGATGAAATCGATGCTGCAAATGAGACTATCACAGCCTTAGAGCACAAAGCGAAGATTGGCAAAATAGACTGTTCAGTAGAGATTACGTCACTTAAGT  >  NC_002947/4181255‑4181549

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