Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,970,525 A→G 8.9% W257R (TGG→CGG)  PP_5236 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,970,5250AG 8.9% 152.3 / 5.2 56W257R (TGG→CGG) PP_5236hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (24/27);  new base G (3/2);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

ACGGCCTGCACCTTGTCAGCCTGTGGCAGCTTCAGCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGCCCAGATAGCGCGGCAGCGACTGTGGCAGC  >  NC_002947/5970383‑5970673
                                                                                                                                              |                                                                                                                                                    
aCGGCCTGCACCTTGTCAGCCTGTGGCAGCTTCAGCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCgggg                                                                                                                                                <  1:19343/149‑1 (MQ=255)
    ccTGCACCTTGTCAGCCTGTGGCAGCTTCAGCCGGCCATCGACAGCCATCCAGCGCAGCGAATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGgcag                                                                                                                                             >  4:84636/1‑148 (MQ=255)
     cTGCACCTTGTCAGCCTGTGGCAGCTTCAGCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGgcag                                                                                                                                             >  3:52144/1‑147 (MQ=255)
              gTCAGCCTGTGGCAGCTTCAGCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACACTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTacac                                                                                                                                    <  8:281349/147‑1 (MQ=255)
                 aGCCTGTGGCAGCTTCAGCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGcagcag                                                                                                                                          <  6:7012/138‑1 (MQ=255)
                 aGCCTGTGGCAGCTTCAGCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGcagcag                                                                                                                                          >  5:7012/1‑138 (MQ=255)
                         gCAGCTTCAGCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCGCTGGGgcggca                                                                                                                                           >  6:194173/1‑129 (MQ=255)
                         gCAGCTTCAGCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCGCTGGGgcggca                                                                                                                                           <  5:194173/129‑1 (MQ=255)
                          cAGCTTCAGCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGa                                                                                                                      <  7:16472/149‑1 (MQ=255)
                                cAGCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCtt                                                                                                                >  8:88431/1‑149 (MQ=255)
                                cAGCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCt                                                                                                                 <  5:197585/148‑1 (MQ=255)
                                  gCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGACCGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTc                                                                                                               >  3:255308/1‑148 (MQ=255)
                                  gCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTc                                                                                                               >  3:140294/1‑148 (MQ=255)
                                  gCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTc                                                                                                               >  4:171479/1‑148 (MQ=255)
                                    cGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGaa                                                                                                            <  4:140294/149‑1 (MQ=255)
                                     ggCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAag                                                                                                           <  2:170863/149‑1 (MQ=255)
                                               aGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTAcaca                                                                                                                                   <  8:103612/115‑1 (MQ=255)
                                               aGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTAcaca                                                                                                                                   >  7:103612/1‑115 (MQ=255)
                                                 ccATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGc                                                                                                                                 >  8:259045/1‑115 (MQ=255)
                                                       agcgcagcgcATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATc                                                                                         <  7:99560/149‑1 (MQ=255)
                                                        gcgcagcgcaTTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGgggcgg                                                                                               <  8:52404/142‑1 (MQ=255)
                                                        gcgcagcgcaTTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGgggcgg                                                                                               >  7:52404/1‑142 (MQ=255)
                                                        gcgcagcgcaTTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCa                                                                                        <  3:171479/149‑1 (MQ=255)
                                                             gcgcaTTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCa                                                                                   <  7:88431/149‑1 (MQ=255)
                                                               gcaTTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGaa                                                                                                            <  4:146373/122‑1 (MQ=255)
                                                               gcaTTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGaa                                                                                                            >  3:146373/1‑122 (MQ=255)
                                                                  ttACACAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGgcg                                                                              >  1:323138/1‑149 (MQ=255)
                                                                            gCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTaa                                                                    <  3:84636/149‑1 (MQ=255)
                                                                                   ccAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGGATGGCGCAGCCAGTAATAGCGAg                                                             <  6:109728/149‑1 (MQ=255)
                                                                                   ccAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCAGCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGa                                                              <  4:255308/148‑1 (MQ=255)
                                                                                            gcACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGt                                                    >  3:64207/1‑149 (MQ=255)
                                                                                                   gggCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATg                                             <  2:323138/149‑1 (MQ=255)
                                                                                                    ggCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATgg                                            <  4:64207/149‑1 (MQ=255)
                                                                                                     gCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGt                                           >  8:204414/1‑149 (MQ=255)
                                                                                                       cATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAgg                                                                                                                             <  2:323742/65‑1 (MQ=255)
                                                                                                       cATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAgg                                                                                                                             >  1:323742/1‑65 (MQ=255)
                                                                                                       cATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGa                                                              <  5:277530/128‑1 (MQ=255)
                                                                                                       cATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGa                                                              <  1:139699/128‑1 (MQ=255)
                                                                                                       cATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGa                                                              >  6:277530/1‑128 (MQ=255)
                                                                                                        aTTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGa                                        >  3:242873/1‑149 (MQ=255)
                                                                                                        aTTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATCGTCGa                                        >  6:75102/1‑149 (MQ=255)
                                                                                                               aaCAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCtt                                 <  7:173657/149‑1 (MQ=255)
                                                                                                               aaCAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCtt                                 <  4:52144/149‑1 (MQ=255)
                                                                                                                    cTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGTGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAgg                                                            <  4:291064/117‑1 (MQ=255)
                                                                                                                    cTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGTGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAgg                                                            >  3:291064/1‑117 (MQ=255)
                                                                                                                              tCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGCCCAGATAGCGCgg                  >  2:55144/1‑149 (MQ=255)
                                                                                                                                 aGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGc                                                                          >  4:254597/1‑90 (MQ=255)
                                                                                                                                 aGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGc                                                                          <  3:254597/90‑1 (MQ=255)
                                                                                                                                    ctgctgCCCCGCTGGGGCGGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGccc                             >  5:194273/1‑132 (MQ=255)
                                                                                                                                    ctgctgCCCCGCTGGGGCGGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGccc                             <  6:194273/132‑1 (MQ=255)
                                                                                                                                    ctgctgCCCCGCTGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGCCCAGATAGCGCGGCAGCGa            >  3:149826/1‑149 (MQ=255)
                                                                                                                                     tGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGCCCAGATAGCGCGGCAGCGAc           >  4:107601/1‑149 (MQ=255)
                                                                                                                                     tGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGCCCAGATAGCGCGGCAGCGAc           >  8:284424/1‑149 (MQ=255)
                                                                                                                                         gccCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGCCCAGATAGCGCGGCAGCGACtg         <  8:161224/147‑1 (MQ=255)
                                                                                                                                         gccCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGCCCAGATAGCGCGGCAGCGACtg         >  7:161224/1‑147 (MQ=255)
                                                                                                                                          ccccAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGCCCAGATAGCGCGGCAGCGACTGTgg      <  7:204414/149‑1 (MQ=255)
                                                                                                                                              aGCGCGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGATGCAGATCAGAACATGGCGCAGACAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGCCCAGATAGCGCGGCAGCGACTGTGGCAGc  >  5:247493/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
ACGGCCTGCACCTTGTCAGCCTGTGGCAGCTTCAGCCGGCCTTCGACAGCCATCCAGCGCAGCGCATTACGCAAACGCTCGTCCCAGGTTGCGCACAGCGGGCCATTGACCAACAGCTCGAACGCATCAAGGCTGCCGCCCCAGCGGGGCAGCAGGTACACAGCAAGGGCCTGGATCGCTTCGAAGCAGCCGGGGCGGGCGATCAGCGCATGGCGCAGCCAGTAATAGCGAGGGTCCTCGTTTTCATGGTCGACCCGCTTGCCCAGATAGCGCGGCAGCGACTGTGGCAGC  >  NC_002947/5970383‑5970673

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