Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,112,033 T→A 11.1% H118Q (CAT→CAA PP_5362 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,112,0330TA11.1% 97.6 / 5.8 36H118Q (CAT→CAAPP_5362hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (19/13);  new base A (2/2);  total (21/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.41e-01

GGGCACAAACCAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGAAAAGCTATGCC  >  NC_002947/6111885‑6112175
                                                                                                                                                    |                                                                                                                                              
gggCACAAACCAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCAtt                                                                                                                                               >  2:377328/1‑150 (MQ=255)
      aaaCCAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCt                                                                                                                                         <  1:179530/150‑1 (MQ=255)
       aaCCAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCtt                                                                                                                                        >  2:280440/1‑150 (MQ=255)
                      tGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCg                                                                                                                         >  2:110581/1‑150 (MQ=255)
                              tGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGcc                                                                                                                 <  1:377328/150‑1 (MQ=255)
                                         aCCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCAttt                                                                                                      <  2:478324/150‑1 (MQ=255)
                                                    cAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGttttt                                                                                           >  2:28318/1‑150 (MQ=255)
                                                          tGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTctg                                                                                     <  1:110581/150‑1 (MQ=255)
                                                                 gCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTa                                                                              >  1:60047/1‑150 (MQ=255)
                                                                   cTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCa                                                                                                                                            >  2:409803/1‑86 (MQ=255)
                                                                   cTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCa                                                                                                                                            <  1:409803/86‑1 (MQ=255)
                                                                       gtgggCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCGGCGGTGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGt                                                                        <  2:500420/150‑1 (MQ=255)
                                                                               cAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGc                                                                >  1:344674/1‑150 (MQ=255)
                                                                                             gCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGTGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTATTCCAGTCGGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCt                                                  >  1:29179/1‑150 (MQ=255)
                                                                                             gCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACgcgc                                                   <  2:536311/149‑1 (MQ=255)
                                                                                             gCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACgcgc                                                   >  1:536311/1‑149 (MQ=255)
                                                                                                 gACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGaaa                                              >  2:376351/1‑150 (MQ=255)
                                                                                                 gACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGaaa                                              >  2:272568/1‑150 (MQ=255)
                                                                                                  accaccTTGAACCGAGTTTCCAGTTCCAGCAACATTCCCAGCAAGGTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAACACGCGCTGaaaa                                             <  2:344674/150‑1 (MQ=255)
                                                                                                    caccTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCtt                                                              >  2:367068/1‑131 (MQ=255)
                                                                                                    caccTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCtt                                                              <  1:367068/131‑1 (MQ=255)
                                                                                                       cTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCa                                        >  1:167956/1‑150 (MQ=255)
                                                                                                                    tCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTtctc                                                                                       <  1:130375/90‑1 (MQ=255)
                                                                                                                    tCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTtctc                                                                                       >  2:130375/1‑90 (MQ=255)
                                                                                                                    tCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTg                                                                         <  1:508409/104‑1 (MQ=255)
                                                                                                                    tCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTg                                                                         >  2:508409/1‑104 (MQ=255)
                                                                                                                      cAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGATTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCt                         >  2:210296/1‑150 (MQ=255)
                                                                                                                         ttCCAGCGGGGAGCGCTTGAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGATCATCCGCTGAAAACACCAGCCTTGAACACCACttct                      >  2:119247/1‑148 (MQ=255)
                                                                                                                            cAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCa                   <  1:280440/150‑1 (MQ=255)
                                                                                                                                     ctcgatCAAGCTGCAATGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGAc                                <  2:21363/123‑1 (MQ=255)
                                                                                                                                     ctcgatCAAGCTGCAATGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGAc                                >  1:21363/6‑128 (MQ=255)
                                                                                                                                         ttCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGAAAAGCTa      >  1:506897/1‑150 (MQ=255)
                                                                                                                                          tCAAGCTGCAATGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGa                                                           <  1:388456/96‑1 (MQ=255)
                                                                                                                                          gcAAGCTGCAATGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGa                                                           >  2:388456/2‑96 (MQ=255)
                                                                                                                                             aGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGAAAAGCTATGcc  >  1:126398/1‑150 (MQ=255)
                                                                                                                                                 gCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAAc                            >  1:52649/1‑120 (MQ=255)
                                                                                                                                                 gCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAAc                            <  2:52649/120‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
GGGCACAAACCAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGAAAAGCTATGCC  >  NC_002947/6111885‑6112175

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