Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,995,682 C→T 7.7% R469H (CGC→CAC)  PP_5252 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,995,6820CT 7.7% 144.3 / 3.0 52R469H (CGC→CAC) PP_5252hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (26/22);  new base T (2/2);  total (28/24)
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.15e-01

TGTCACGGCTCAGGCCTTCCGGGTACAGCTCCATGCCACCGACCGTCTTGCCGCTGACCAGCCAGTACAGTGAGGTCCATGGGTTGTAGCTGGACACCCGAGTGGCATCGGTACCCGCGCCAACCGGCACGCCCATCTCGAGCATGCGTTTGATCGGCGGGGTCTGCTCGGCGGCCTTGGCGCCGTAGCGGTCGACGAAGTACTCACCCTGGAAGGCCATGCGGTCCTGGATGGCAATACCACCGCCCAGCGCACGCACGCGCTCGATGTTCTGCGGGGTGATGGTTTCGGC  >  NC_002947/5995536‑5995827
                                                                                                                                                  |                                                                                                                                                 
tgtCACGGCTCAGGCCTTCCGGGTACAGCTCCATGCCACCGACCGTCTTGCCGCTGACCAGCCAGTACAGTGAGGTCCATGGGTTGTAGCTGGACACCCGAGTGGCATCGGTACCCGCGCCAACCGGCACTCCCATCTCGAGCATGCGtt                                                                                                                                                <  2:943104/150‑1 (MQ=255)
    aCGGCTCAGGCCTTCCGGGTACAGCACCATGCCACCGACCGTCTTGCCGCTGACCAGCCATTACAGTGAGGTCCATGGGTTGTAGCTGGACACCCGAGTGGCATCGGTACCCGAGCCAACCGGGTCGCCCATCTCGAGCATGCGTTTGAt                                                                                                                                            >  1:261339/1‑150 (MQ=255)
     cGGCTCAGGCCTTCCGGGTACAGCTCCATGCCACCGACCCTCTTGCCGCTGACCAGCCAGTACAGTGAGGTCCATGGGTTGTAGCTGGACACCCGAGTGGCATCGGTACCCGCGCCAACTGGCACGCCCATCTCGAGCATGCGTTTGATc                                                                                                                                           <  2:654467/150‑1 (MQ=255)
            ggCCTTCCGGGTACAGCTCCATGCCACCGACCGTCTTGCCGCTGACCAGCCAGTACAGTGAGGTCCATGGGTTGTAGCTGGACACCCGAGTGGCATCGGTACCCGCGCCAACCGGCACGCCCATCTCGAGCATGCGTTTGATCGGCgggg                                                                                                                                    <  1:692619/150‑1 (MQ=255)
                       tACAGCTCCATGCCACCGACCGTCTTGCCGCTGACCAGCAAGTACAGTGAGGTCCATGGGTTGTAGCTGGACACCCGAGTGGCATCGGTACCCGCGCCAACCGGCACGCCCATCTCGAGCATGCGTTTGATCGGCGGGGTCTGCTcggcg                                                                                                                         <  2:790118/150‑1 (MQ=255)
                               cATGCCACCGACCGTCTTGCCGCTGACCAGCCAGTACAGTGAGGTCCATGGGTTGTAGCTGGACACCCGAGTGGCATCGGTACCCGCGCCAACCGGCACGCCCATCTCGAGCATGCGTTTGATCGGCGGGGTCTGCTCGGCGGCCTTGgc                                                                                                                 <  2:642662/150‑1 (MQ=255)
                                aTGCCACCGACCGTCTTGCCGCTGACCAGCCAGTACAGTGAGGTCCATGGGTTGTAGCTGGACACCCGAGTGGCATCGGTACCCGCGCCAACCGGCACGCCCATCTCGAGCATGCGTTTGATCGGCggg                                                                                                                                     >  2:969998/1‑129 (MQ=255)
                                aTGCCACCGACCGTCTTGCCGCTGACCAGCCAGTACAGTGAGGTCCATGGGTTGTAGCTGGACACCCGAGTGGCATCGGTACCCGCGCCAACCGGCACGCCCATCTCGAGCATGCGTTTGATCGGCggg                                                                                                                                     <  1:969998/129‑1 (MQ=255)
                                aTGCCACCGACCGTCTTGCCGCTGACCAGCCAGTACAGTGAGGTCCATGGGTTGTAGCTGGACACCCGAGTGGCATCGGTACCCGCGCCAACCGGCACGCCCATCTCGAGCATGCGTTTGATCGGCGGGGTCTGCTCGGCGGCCTTGgcg                                                                                                                >  2:562974/1‑150 (MQ=255)
                                         accgTCTTGCCGCTGACCAGCCAGTACAGTGAGGTCCATGGGTTGTAGCTGGACACCCGAGTGGCATCGGTACCCGCGCCAACCGGCACGCCCATCTCGAGCATGCGTTTGATcggcgg                                                                                                                                      >  1:13322/1‑119 (MQ=255)
                                         accgTCTTGCCGCTGACCAGCCAGTACAGTGAGGTCCATGGGTTGTAGCTGGACACCCGAGTGGCATCGGTACCCGCGCCAACCGGCACGCCCATCTCGAGCATGCGTTTGATcggcgg                                                                                                                                      <  2:13322/119‑1 (MQ=255)
                                             tCTTGCCGCTGACCAGCCAGTACAGTGAGGTCCATGGGTTGTAGCTGGACACCCGAGTGGCATCGGTACCCGCGCCAACCGGCACGCCCATCTCGAGCATGCGTTTGATCGGCGGGGTCTGCTCGGCGGCCTTGGCGCCGTAGCGGTcga                                                                                                   <  2:464732/150‑1 (MQ=255)
                                                    gCTGACCAGCCAGTACAGTGAGGTCCATGGGTTGTAGCTGGACACCCGAGTGGCATCGGTACCCGCGCCAACCGGCACGCCCATCTCGAGCATGCGTTTGATCGGCGGGGTCTGCTCGGCGGCCTTGGCGCCGTAGCGGTCGACGAAGTa                                                                                            >  2:272146/1‑150 (MQ=255)
                                                          cagccagTACAGTGAGGTCCATGGGTTGTAGCTGGACACCCGAGTGGCATCGGTACCCGCGCCAACCGGCACGCCCATCTCGAGCATGCGTTTGATCGGCGGGGTCTGCTCGGCGGCCTTGGCGCCGTAGCGGTCGACGAAGTACTCAcc                                                                                      >  2:21488/1‑150 (MQ=255)
                                                          cagccagTACAGTGAGGTCCATGGGTTGTAGCTGGACACCCGAGTGGCATCGGTACCCGCGCCAACCGGCACGCCCATCTCGAGCATGCGTTTGATCGGCGGGGTCTGCTCGGCGGCCTTGGCGCCGTAGCGGTCGACGAAGTACTCAcc                                                                                      >  1:648245/1‑150 (MQ=255)
                                                                 tACAGTGAGGTCCATGGGTTGTAGCTGGACACCCGAGTGGCATCGGTACCCGCGCCAACCGGCACGCCCATCTCGAGCATGCGTTTGATCGGCGGGGTCTGCTCGGCGGCCTTGGCGCCGTAGCGGTCGACGAAGTACTCACCCTGCAAg                                                                               <  2:706091/150‑1 (MQ=255)
                                                                          gTCCATGGGTTGTAGCTGGACACCCGAGTGGCATCGGTACCCGCGCCAACCGGCACGCCCATCTCGAGCATGCGTTTGATCGGCGGGGTCTGCTCGGCGGCCTTGGCGCCGTAGCGGTCGACGAAg                                                                                              <  1:219279/126‑1 (MQ=255)
                                                                          gTCCATGGGTTGTAGCTGGACACCCGAGTGGCATCGGTACCCGCGCCAACCGGCACGCCCATCTCGAGCATGCGTTTGATCGGCGGGGTCTGCTCGGCGGCCTTGGCGCCGTAGCGGTCGACGAAg                                                                                              >  2:219279/1‑126 (MQ=255)
                                                                            ccATGGGTTGTAGCTGGACACCCGAGTGGCATCGGTGCCCGCGCCAACCGGCACGCCCATCTCGAGCATGCGTTTGATCGGCGGGGTCTGCTCGGCGGCCTTGGCGCCGTAGCGGTCGACGAAGTACTCACCCTGGAAGGCCATGCGGTc                                                                    >  1:1058072/1‑150 (MQ=255)
                                                                              aTGGGTTGTAGCTGGACACCCGAGTGGCATCGGTACCCGCGCCAACCGGCACGCCCATCTCGAGCATGCGTTTGATCGGCGGGGTCTGCTCGGCGGCCTTGGCGCCGTAGCGGTCGACGAAGTACTCACCCTGGAAGGCCATGCGGTCCt                                                                  >  1:1038229/1‑150 (MQ=255)
                                                                                   ttGTAGCTGGACACCCGAGTGGCATCGGTACCCGCGCCAACCGGCACGCCCATCTCGAGCATGCGTTTGATCGGCGGGGTCTGCTCGGCGGCCTTGGCGCCGTAGCGGTCGACGAAGTACTCACCCTGGAAGGCCATGCGGTCCTGGATg                                                             >  2:722459/1‑150 (MQ=255)
                                                                                    tGTAGCTGGACACCCGAGTGGCATCGGTACCCGCGCCAACCGGCACGCCCATCTCGAGCATGCGTTTGATCGGCGGGGTCTGCTCGGCGGCCTTGGCGCCGTAGCGGTCGACGAAGTACTCACCCTGGAAGGCCATGCGGTCCTGGATgg                                                            <  1:562974/150‑1 (MQ=255)
                                                                                     gTAGCTGGACACCCGAGTGGCATCGGTACCCGCGCCAACCGGCACGCCCATCTCGAGCATGCGTTTGATCGGCGGGGTCTGCTCGGCGGCCTTGGCGCCGTAGCGGTCGACGAAGTACTCAccc                                                                                     >  1:660238/1‑124 (MQ=255)
                                                                                     gTAGCTGGACACCCGAGTGGCATCGGTACCCGCGCCAACCGGCACGCCCATCTCGAGCATGCGTTTGATCGGCGGGGTCTGCTCGGCGGCCTTGGCGCCGTAGCGGTCGACGAAGTACTCAccc                                                                                     <  2:660238/124‑1 (MQ=255)
                                                                                       aGCTGGACACCCGAGTGGCATCGGTACCCGCGCCAACCGGCACGCCCATCTCGAGCATGCGTTTGATCGGCGGGGTCTGCTCGGCGGCCTTGGCGCCGTAGCGGTCGACGAAGTACTCACCCTGGAAGGCCATGCGGTCCTGGATGGCaa                                                         <  2:1038229/150‑1 (MQ=255)
                                                                                        gCTGGACACCCGAGTGGCATCGGTACCCGCGCCAACCGGCACGCCCATCTCGAGCATGCGTTTGATCGGCGGGGTCTGCTCGGCGGCCTTGGCGCCGTAGCGGTCGACGAAGTACTCACCCTGGAAGGCCATGCGGTCCTGGATGGCAAt                                                        >  1:81815/1‑150 (MQ=255)
                                                                                            gACACCCGAGTGGCATCGGTACCCGCGCCAACCGGCACGCCCATCTCGAGCATGCGTTTGATCGGCGGGGTCTGCTCGGCGGCCTTGGCGCCGTAGCGGTCGACGAAGTACTCACCCTGGAAGGCCATGCGGTCCTGGATGGCAATacca                                                    <  2:81815/150‑1 (MQ=255)
                                                                                              cacCCGAGTGGCATCGGTACCCGCGCCAACTGGCAGGCCCATCTCGTGCATGCGTTTGATCGGCGGGGTCTGCTCGGCGGCCTTGGCGCCTTATCGGTCGACGAAGTACTCACCCTGGAAGGCCTTGCTGTTAATGATGCCACTaccacc                                                  <  2:363228/150‑1 (MQ=255)
                                                                                                cccGAGTGGCATCGGTACCCGCGCCAACCGGCACGCCCATCTCGAGCATGCGTTTGATCGGCGGGGTCTGCTCGGCGGCCTTGGCGCCGTAGCGGTCGACGAAGTACTCACCCTGGAAGGCCATGCGGTCCTGGATGGCAATACCACCGc                                                >  2:71076/1‑150 (MQ=255)
                                                                                                 ccGAGTGGCATCGGTACCCGCGCCAACCGGCACGCCCATCTCGAGCATGCGTTTGATCGGCGGGGTCTGCTCGGCGGCCTTGGCGCCGTAGCGGTCGACGAAGTACTCACCCTGGAAGGCCATGCGGTCCTGGATGGCAATACCACCGcc                                               >  2:863727/1‑150 (MQ=255)
                                                                                                 ccGAGTGGCATCGGTACCCGCGCCAACCGGCACGCCCATCTCGAGCATGCGTTTGATCGGCGGGGTCTGCTCGGCGGCCTTGGCGCCGTAGCGGTCGACGAAGTACTCACCCTGGAAGGCCATGCGGTCCTGGATCGCAATACCACAGcc                                               >  2:382759/1‑150 (MQ=255)
                                                                                                 ccGAGTGGCATCAGTACACGCGCCAACCTGCACTCCCATCTCGAGCATGCGTTTGATCGGCGGGGTCTGCTCGGCGGCCTTGGCGCCGACAAGGTCAATGTCTGACTCACCCTTCTATTCCATGCGGTCCTTGATGGCAATACCACCGcc                                               >  2:242914/1‑150 (MQ=255)
                                                                                                  cGAGTGGCATCGGTACCCGCTCCAACCGGCACGCCCATCTCGAGCATGCGGTTGATCGGCGGGGTCTGCTCGGCGGCCTTGGCGCCGTAGCGGTCGACGAAGTACTCACCCTGGAAGGCCATGCGGTCCTGGATGGCAATACCACCGccc                                              >  2:156052/1‑150 (MQ=255)
                                                                                                  cGAGTGGCATCGGTACCCGCGCCAACCGGCACGCCCATCTCGAGCATGCGTTTGATCGGCGGGGTCTGCTCGGCGGCCTTGGCGCCGTAGCGGTCGACGAAGTACTCACCCTGGAAGGCCATGCGGTCCTGGATGGCAATACCACCGccc                                              <  2:72583/150‑1 (MQ=255)
                                                                                                             ggTACCCGCGCCAACCGGCACGCCCATCTCGAGCATGCGTTTGATCGGCGGGGTCTGCTCGGCGGCCTTGGCGCCGTAGCGGTCGACGAAGTACTCACCCTGGAAGGCCATGCGGTCCTGGATGGCAATACCACCGCCCAGCgcacgcac                                   <  1:71076/150‑1 (MQ=255)
                                                                                                              gTACCCGCGCCAACCGGCACGCCCATCTCGAGCATGCGTTTGATCGGCGGGGTCTGCTCGGCGGCCTTGGCGCCGTAGCGGTCGACGAAGTACTCACCCTGGAAGGCCATGCGGTCCTGGATGGCAATACCACCGCCCa                                             >  2:373035/1‑139 (MQ=255)
                                                                                                              gTACCCGCGCCAACCGGCACGCCCATCTCGAGCATGCGTTTGATCGGCGGGGTCTGCTCGGCGGCCTTGGCGCCGTAGCGGTCGACGAAGTACTCACCCTGGAAGGCCATGCGGTCCTGGATGGCAATACCACCGCCCa                                             <  1:373035/139‑1 (MQ=255)
                                                                                                                  ccGCGCCAACCGGCACGCCCATCTCGAGCATGCGTTTGATCGGCGGGGTCTGCTCGGCGGCCTTGGCGCCGTAGCGGTCGACGAAGTACTCACCCTGGAAGGCCATGCGGTCCTGGATGGCAATACCACCGCCCAGCGCACGCACGCGCt                              >  1:379701/1‑150 (MQ=255)
                                                                                                                             ggCACGCCCATCTCGAGCATGCGTTTGATCGGCGGGGTCTGCTCGGCGGCCTTGGCGCCGTAGCGGTCGACGAAGTACTCACCCTGGAAGGCCATGCGGTCCTGGATGGCAATACCACCGCCCAGCGCACGCACGCGCTCGATGTTCTGc                   >  2:388026/1‑150 (MQ=255)
                                                                                                                                aCGCCCATCTCGAGCATGCGTTTGATCGGCGGGGTCTGCTCGGCTGCCTTGGCGCCGTAGCGGTCGACGAAGTACTCACCCTGGAAGGCCATGCGGTCCTGGATGGCAATACCACCGCCCAGCGCACGCACGCGCTCGATGTTCTGCggg                >  2:157193/1‑150 (MQ=255)
                                                                                                                                  gCCCATCTCGAGCATGCGTTTGATCGGCGGGGTCTGCTCGGCGGCCTTGGCGCCGTGGCGGTCGACGAAGTACTCACCCTGGGAGGCCATGCGGTCCTGGATGGCAATACCACCGCCCAGCGCACGCACGCGCTCGATGTTCTGCGGGGt              >  2:1120617/1‑150 (MQ=255)
                                                                                                                                    ccATCTCGAGCATGCGTTTGATCGGCGGGGTCTGCTCGGCGGCCTTGGCGCCGTAGCGGTCGACGAAGTACTCACCCTGGAAGGCCATGCGGTCCTGGATGGCAATACCACCGCCCAGCGCACGCACGCGCTCGATGTTCTGCGGGGTGa            <  1:485715/150‑1 (MQ=255)
                                                                                                                                      aTCTCGAGCATGTGTTTGATCGGCGGGGTCTGCTCGGCGGCCTTGGCGCCGTAGCGGTCGACGAAGTACTCACCCTGGAAGGCCATGCGGTCCTGGATGGCaa                                                         >  2:626793/1‑103 (MQ=255)
                                                                                                                                      aTCTCGAGCATGTGTTTGATCGGCGGGGTCTGCTCGGCGGCCTTGGCGCCGTAGCGGTCGACGAAGTACTCACCCTGGAAGGCCATGCGGTCCTGGATGGCaa                                                         >  2:626765/1‑103 (MQ=255)
                                                                                                                                      aTCTCGAGCATGTGTTTGATCGGCGGGGTCTGCTCGGCGGCCTTGGCGCCGTAGCGGTCGACGAAGTACTCACCCTGGAAGGCCATGCGGTCCTGGATGGCaa                                                         <  1:626765/103‑1 (MQ=255)
                                                                                                                                      aTCTCGAGCATGTGTTTGATCGGCGGGGTCTGCTCGGCGGCCTTGGCGCCGTAGCGGTCGACGAAGTACTCACCCTGGAAGGCCATGCGGTCCTGGATGGCaa                                                         <  1:626793/103‑1 (MQ=255)
                                                                                                                                      aTCTCGAGCATGCGTTTGATCGGCGGGGTCTGCTCGGCGGCCTTGGCGCCGTAGCGGTCGACGAAGTACTCACCCTGGAAGGCCATGCGGTCCTGGATGGCAATAACACCGCCCAGCGCACGCACGCGCTCGATGTTCTGCGGGGTGATg          >  2:361724/1‑150 (MQ=255)
                                                                                                                                       tctcGAGCATGCGTTTGATCGGCGGGGTCTGCTCGGCGGCCTTGGCGCCGTAGCGGTCGACGAAGTACTCACCCTGGAAGGCCATGCGGTCCTGGATGGCAATACCACCGCCCAGCGCACGCACGCGCTCGATGTTCTGCGGGGTGATgg         <  2:786938/150‑1 (MQ=255)
                                                                                                                                         tcGAGCATGCGTTTGATCGGCGGGGTCTGCTCGGCGGCCTTGGCGCCGTAGCGGTCGACGAAGTACTCACCCTGGAAGGCCATGCGGTCCTGGATGGCAATACCACCGCCCAGCGCACGCACGCGCTCGATGTTCTGCGGGGTGATGGtt       >  2:176744/1‑150 (MQ=255)
                                                                                                                                          cGAGCATGCGTTTGATCGGCGGGGTCTGCTCGGCGGCCTTGGCGCCGTAGCGGTCGACGAAGTACTCACCCTGGAAGGCCATGCGGTCCTGGATGGCAATACCACCGCCCAGCGCACGCACGCGCTCGATGTTCTGCGGGGTGATGGttt      <  1:94806/150‑1 (MQ=255)
                                                                                                                                              cATGCGTTTGATCGGCGGGGTCTGCTCGGCGGCCTTGGCGCCGTAGCGGTCGACGAAGTACTCACCCTGGAAGGCCATGCGGTCCTGGATGGCAATACCACCGCCCAGCGCACGCACGCGCTCGATGTTCTGCGGGGTGATGGTTTCGGc  >  1:921965/1‑150 (MQ=255)
                                                                                                                                              cATGCGTTTGATCGGCGGGGTCTGCTCGGCGGCCTTGGCGCCGTAGCGGTCGACGAAGTACTCACCCTGGAAGGCCATGCGGTCCTGGATGGCAATACCACCGCCCAGCGCACGCACGCGCTCGATGTTCTGCGGGGTGATGGTTTCGGc  <  1:863727/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TGTCACGGCTCAGGCCTTCCGGGTACAGCTCCATGCCACCGACCGTCTTGCCGCTGACCAGCCAGTACAGTGAGGTCCATGGGTTGTAGCTGGACACCCGAGTGGCATCGGTACCCGCGCCAACCGGCACGCCCATCTCGAGCATGCGTTTGATCGGCGGGGTCTGCTCGGCGGCCTTGGCGCCGTAGCGGTCGACGAAGTACTCACCCTGGAAGGCCATGCGGTCCTGGATGGCAATACCACCGCCCAGCGCACGCACGCGCTCGATGTTCTGCGGGGTGATGGTTTCGGC  >  NC_002947/5995536‑5995827

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