Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,530,959 A→C 9.7% T66T (ACA→ACC)  PP_2219 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,530,9590AC 9.7% 123.7 / 2.8 48T66T (ACA→ACC) PP_2219hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (23/20);  new base C (3/2);  total (26/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

AGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACAGCAGCTGAGCTGGGGTGCCGACCGGCGGATTGTGAGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCGCCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGAT  >  NC_002947/2530814‑2531107
                                                                                                                                                 |                                                                                                                                                    
agatcatgatcaagatcaagagcaagatgaagatcaagatcaagatcaagatctACGGCTAGGGCAACGGCAACGGCAACGGCAACGTCAAGCGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCAcagca                                                                                                                                                   <  3:325804/149‑1 (MQ=25)
    caagatcaagatcaagatcaagatcaagatcaagatcaagatcaagatcaaCGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCACCAGGAGGGGGATGGGCACAGAATCTCAGTTTACTACTGGCACAGCagct                                                                                                                                                >  6:134822/1‑148 (MQ=255)
       gatcaagatcaagatcaagatcaagatcaagatcaagatcaagatcaaCGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACAGCagctgagc                                                                                                                                            >  2:25700/1‑149 (MQ=255)
          caagatcaagatcaagatcaagatcaagatcaagatcaagctcaaCGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACAGCAGCTGAGCTgg                                                                                                                                         >  1:17790/1‑149 (MQ=255)
          caagatcaagatcaagatcaagatcaagatcaagatcaagatcaaCGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACAGCAGCTGAGCTgg                                                                                                                                         >  2:77294/1‑149 (MQ=255)
           aagatcaagatcaagatcaagatcaagatcaagatcaagatcaaCGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACAGCAGCTGAGCTggg                                                                                                                                        <  8:247026/149‑1 (MQ=255)
               tcaagatcaagatcaagatcaagatcaagatcaagatcaaCGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACAGCAGCTGAGCTGGGGTGc                                                                                                                                    <  5:183441/149‑1 (MQ=255)
                    atcaagatcaagatcaagatcaagatcaagatcaaCGGCAACGGCAACGGCAACGGCAACGGCAACGGCACCGGCAACGGCAACAGCAGGTGGAGGGGCACAGACTCTCAGTTGACTCCTGGCCCCGCAGCTGAGCGGGGGTGCCGAcc                                                                                                                               >  4:160808/1‑149 (MQ=21)
                     tcaagatcaagatcaagatcaagatcaagatcaaCGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACAGCAGCTGAGCTGGGGTGCCGACcg                                                                                                                              >  6:155561/1‑149 (MQ=255)
                        agatcaagatcaagatcaagatccagatcaaCGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACAGCAGCTGAGCTGGGGTGCCGACcggcg                                                                                                                           <  3:160808/149‑1 (MQ=255)
                          atcaagctcaagatcaagatcaagatcaaCGGCAACCGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACAGCAGCTGAGCTGGGGTGCCGACCGGCGGa                                                                                                                         <  6:304792/149‑1 (MQ=255)
                                atcaagatcaagatcaagatcaaCGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACAGCAGCTGAGCTGGGGTGCCGACCGGCGGATtgt                                                                                                                     >  4:136149/1‑147 (MQ=255)
                                 tcaagatcaagatcaagatcaaCGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACCCCAGCTCAGCTGCTGTGccag                                                                                                                                 >  2:212673/1‑132 (MQ=25)
                                 tcaagatcaagatcaagatcaaCGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACCCCAGCTCAGCTGCTGTGccag                                                                                                                                 <  1:212673/134‑3 (MQ=25)
                                      atcaagatcaagatcaaCGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGCATCTCAGTTTACTACTGGCACAGCAGCTGAGCTGGGGTGCCGACCGGCGGATTGTGAGCAGTg                                                                                                             >  4:38955/1‑149 (MQ=255)
                                         aagatcaagatcaaCGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACAGCAGCTGAGCTGGGGTGCCGACCGGCGGATTGTGAGCAGTGGaa                                                                                                          >  7:110967/1‑149 (MQ=255)
                                                agatcaaCGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACCCCagctcagctgctg                                                                                                                                       <  1:229563/113‑4 (MQ=37)
                                                agatcaaCGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACCCCagctcagctgctg                                                                                                                                       >  2:229563/1‑110 (MQ=37)
                                                agatcaaCGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACAGCAGCTGAGCTGGGGTGCCGACCGGCGGATTGTGAGCAGTGGAAGCAAAtt                                                                                                   <  5:134822/149‑1 (MQ=255)
                                                   tcaaCGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACAGCAGCTGAGCTGGGGTGCCGACCGGCGGATTGTGAGCAGTGGAAGCAAATTCGc                                                                                                <  5:155561/149‑1 (MQ=255)
                                                   tcaaCGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAAAGCAAACAGAAAGTGGAAGGGCAAAGAATCTCAGTTTAATACTGGCACAGGAGCTGAGGTGGGGTGca                                                                                                                                   >  6:195289/1‑113 (MQ=21)
                                                   tcaaCCGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACAGCAGCTGAGCTGGGGTGcc                                                                                                                                   <  5:195289/114‑1 (MQ=255)
                                                     aacggcaacggcaacggcaacggcaacggcaacggcaaccgcaacggcAACAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACAGCAGCTGAGCTGGGGTGCCGACCGGCGGATTGTGAGCAGTGGAAGCAAATTCGCa                                                                                               <  3:38955/148‑1 (MQ=255)
                                                        ggcaaaggcaacggcaaaggcaacggcaacggcaaaggcaacggccacAGCAGGTGGATGGGCACAGCATCTAAGGTTACTACTGGCACAGCAGCTGATCTGGGGTGCAGGCCGGAGGATTGTGAGCAGTGGAAGGAAATTCTCAGCCa                                                                                           >  1:107356/1‑149 (MQ=37)
                                                                 aacggcaacggcaacggcaacggcaacggcaacggcaacAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACAGCAGCTGAGCTGGGGTGCCGACCGGCGGATTGTGAGCAGTGGAAGCAAATTCGCACCCATTTGCAGtt                                                                                  >  3:305324/1‑149 (MQ=255)
                                                                     gcaacggcaacggcaacggcaacggcaacggcaacAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACAGCAGCTGAGCTGGGGTGCCGACCGGCGGATTGTGAGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTc                                                                              <  1:25700/149‑1 (MQ=255)
                                                                            caacggcaacggcaacggcaacggcaacAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACAGCAGCTGAGCTGGGGTGCCGACCGGCGGATTGTGAGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACgcgg                                                                       <  1:77294/149‑1 (MQ=255)
                                                                             aacggcatcggcaacgtcaacggcaacAGCAGGTGGATGGGCACATAATCTCAGTTTACTACTGGCACAGCAGCTGAGCTGGGGTGCCGACCGGCGGATTGTGAGCAGTGGAAGCGAATTCGCACCCATTTGCAGTTGCTCCACgcggc                                                                      <  3:379987/149‑1 (MQ=255)
                                                                             aacggcaacggcaacggcaacggcaacAGCAGGTGGATGGGCACAGCATCTCAGTTTACTACTGGCACAGCAGCTGAGCTGGGGTGACGACCGGCGGATTGTGAGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACgcggc                                                                      >  2:226256/1‑149 (MQ=255)
                                                                             aacggcaacggcaacggcaacggcaacAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACAGCAGCTGAGCTGGGGTGCCGACCGGCGGATTGTGAGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACgcggc                                                                      >  2:356137/1‑149 (MQ=255)
                                                                                   aagggctacggcaatggcaacAGCAGTTGGATGGGCACAGAATCTCAGTTTCCTAGTGGCACAGCAGCTGAGCTGGGGTGCCGACCGGCGGATTGTGAGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGgc                                                                 <  3:321276/145‑1 (MQ=255)
                                                                                         aacggcaacggcaacAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACAGCAGCTGAGCTGGGGTGCCGACCGGCGGATTGTGAGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAg                                                          <  2:17790/149‑1 (MQ=255)
                                                                                             gcaacggcaacAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACAGCAGCTGAGCTGGGGTGCCGACCGGCGGATTGTGAGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCACCTTAGACTg                                                      <  3:136149/149‑1 (MQ=255)
                                                                                                       aagcagGTGGATGGTCACAGAATCTCATTTTACTACTGGAACAGCAGCTGATCGTGGGTGCCGACCTGAGGATTGTGAGCAGTGGAATCAAATTCGCACCCATTTGCAGTTTCTCCACGCGGCGTGGCGCCTTATACTGCCGCCCCCTg                                            <  2:107356/148‑1 (MQ=255)
                                                                                                        agcagGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACAGCAGCTGAGCTGGGGTGCCGACCGGCGGATTGTGAGCAGTGGAAGCAAATTCGCACCCATTTGc                                                                                      >  6:64916/1‑106 (MQ=255)
                                                                                                        agcagGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACAGCAGCTGAGCTGGGGTGCCGACCGGCGGATTGTGAGCAGTGGAAGCAAATTCGCACCCATTTGc                                                                                      <  5:64916/106‑1 (MQ=255)
                                                                                                         gcagGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACAGCAGCTGAGCTGGGGTGCCGACCGGCGGATTGTGAGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCGCCCCCTGt                                           >  6:197346/1‑148 (MQ=255)
                                                                                                             gTGGATGGGCACAGAATCTCAGTTTACTACTGGCACAGCAGCTGAGCTGGGGTGCCGACCGGCGGATTGTGAGCAGTGGAAGCAAATTAGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCGCCCCCTGTaa                                         >  3:64671/1‑146 (MQ=255)
                                                                                                                gATGGGCACAGAATCTCAGTTTACTACTGGCACAGCAGCTGAGCTGGGGTGCCGACCGGCGGATTGTGAGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCGCCCCCTGTAAAAGAGt                                   >  7:274474/1‑149 (MQ=255)
                                                                                                                     gCACAGAATCTCAGTTTACTACTGGCACAGCAGCTGAGCTGGGGTGCCGACCGGCGGATTGTGAGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCGCCCCCTGTAAAAGAGTGCCgg                              >  6:269604/1‑149 (MQ=255)
                                                                                                                                         actGGCACAGCAGCTGAGCTGGGGTGCCGACCGGCGGATTGTGAGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCGCCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCg          >  5:171078/1‑149 (MQ=255)
                                                                                                                                          ctGGCACAGCAGCTGAGCTGGGGTGCCGACCGGCGGATTGTGAGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCGCCCCCTGTAAAAGTGTGCCGGTTGGCGCTTGAAGAATTCCGc         <  1:226256/149‑1 (MQ=255)
                                                                                                                                          ctGGCACAGCAGCTGAGCTGGGGTGCCGACCGGCGGATTGTGAGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCGCCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCg          <  4:64671/148‑1 (MQ=255)
                                                                                                                                          ctGGCACAGCAGCTGAGCTGGGGTGCCGACCGGCGGATTGTGAGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCGCCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGc         >  2:295271/1‑149 (MQ=255)
                                                                                                                                           tGGCACAGCAGCTGAGCTGGGGTGCCGACCGGCGGATTGTGAGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCGCCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCt        <  8:110967/149‑1 (MQ=255)
                                                                                                                                             gCACAGCAGCTGAGCTGGGGTGCCGACCGGCGGATTGTGAGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCGCCCCCTGTAAAAGAGTGCCGGTTGGCGCTTg                    >  7:156967/1‑135 (MQ=255)
                                                                                                                                             gCACAGCAGCTGAGCTGGGGTGCCGACCGGCGGATTGTGAGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCGCCCCCTGTAAAAGAGTGCCGGTTGGCGCTTg                    <  8:156967/135‑1 (MQ=255)
                                                                                                                                             gCACAGCAGCTGAGCTGGGGTGCCGACCGGCGGATTGTGAGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCGCCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTgg      >  2:181991/1‑149 (MQ=255)
                                                                                                                                             gCACAGCAGCTGAGCTGGGGTGCAGACCGGCGGATTGTGAGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGATCCACGAGGCGTGGCGCATTAGAATGCAGCCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTg       >  2:26177/1‑148 (MQ=255)
                                                                                                                                                 agcagcTGAGCTGGGGTGCCGACCGGCGGATTGTGAGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCGCCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGAt  <  1:356137/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
AGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACGGCAACAGCAGGTGGATGGGCACAGAATCTCAGTTTACTACTGGCACAGCAGCTGAGCTGGGGTGCCGACCGGCGGATTGTGAGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCGCCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGAT  >  NC_002947/2530814‑2531107

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