Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,016,199 C→T 8.3% T14T (ACG→ACA)  PP_0875 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,016,1990CT 8.3% 135.4 / 3.1 50T14T (ACG→ACA) PP_0875hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (24/20);  new base T (2/2);  total (27/23)
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.32e-01

GGGGCCTGGCGCTTGAGGTAGTACGCGGTCCAAAGGCCGGTATAGCCGGCACCGATGATGCACACGTCAACGTCGAGGTCATTGCGCAGTGCCGGGCGCGCGCACAGCGGCTCGTCGAGCTGGTCCATCCACAGGCTGAGCGTGCGCAGGGGTTGCATCGGGTTTGGCTCCACATCCGTCAAGGTCTGTGGGCGATGCTAGTGGTGTCAGTGGGCGGCTGTCTTACGTGCGTGCACGCAGGGAAATTTGCTTCAGGTAGGCCTTGGGCGTGAACCCGGTGTGTTCGCGG  >  NC_002947/1016059‑1016347
                                                                                                                                            |                                                                                                                                                    
ggggCCTGGCGCTTGAGGTAGTACGCGGTCCAAAGGCCGGTATAGCCGGCACCGATGATGCACACGTCAACGTCGAGGTCATTGCGCAGTGCCGGGCGCGCGCACAGCGGCTCGTCGAGCTGGTCCATCCACAGGCTGAGCGTGCGCAg                                                                                                                                              <  6:87129/149‑1 (MQ=255)
      tGGCGCTTGAGGTAGTACGCGGTCCAAAGGCCGGTATAGCCGGCACCGATGATGCACACATCAACGTCGAGGTCATTGCGCAGTGCCGGGCGCGCGCACAGCGGCTCGTCGAGCTGGTCCATCCACAGGCTGAGCGTGCGCAGGGGTTg                                                                                                                                        >  1:105326/1‑149 (MQ=255)
          gcTTGAGGTAGTACGCGGTCCAAAGGCCGGTATAGCCGGCACCGATGATGCACACGTCAACGTCGAGGTCATTGCGCAGTGCCGGGCGCGCGCACAGCGGCTCGTCGAGCTGGTCCATCCACAGGCTGAGCGTGCGCAGGGGTTGCATc                                                                                                                                    >  8:134760/1‑149 (MQ=255)
                    gtaCGCGGTCCAAAGGCCGGTATAGCCGGCACCGATGATGCACACGTCAACGTCGAGGTCATTGCGCAGTGCCGGGCGCGCGCACAGCGGCTCGTCGAGCTGGTCCATCCACAGGCTGAGCGTGCGCAGGGGTTGCATCGGGTTTGGc                                                                                                                           <  6:374407/148‑1 (MQ=255)
                      aCGCGGTCCAAAGGCCGGTATAGCCGGCACCGATGATGCACACGTCAACGTCGAGGTCATTGCGCAGTGCCGGGCGCGCGCACAGCGGCTCGTCGAGCTGGTCCATCCACAGGCTGAGCGTGCGCAGGGGTTGCATCGGGTTTGGCTc                                                                                                                         <  7:342260/148‑1 (MQ=255)
                            tCCAAAGGCCGGTATAGCCGGCACCGATGATGCACACGTCAACGTCGAGGTCATTGCGCAGTGCCGGGCGCGCGCACAGCGGCTCGTCGAGCTGGTCCATCCACAGGCTGAGCGTGCGCAGGGGTTGCATCGGGTTTGGCTCCACATcc                                                                                                                  >  4:129583/1‑149 (MQ=255)
                               aaaGGCCGGTATAGCCGGCACCGATGATGCACACGTCAACGTCGAGGTCATTGCGCAGTGCCGGGCGCGCGCACAGCGGCTCGTCGAGCTGGTCCATCCACAGGCTGAGCGTGCGCAGGGGTTGCATCGGGTTTGGCTCCACATCCGt                                                                                                                >  1:234803/1‑148 (MQ=255)
                                  ggCCGGTATAGCCGGCACCGATGATGCACACGTCAACGTCGAGGTCATTGCGCAGTGCCGGGCGCGCGCACAGCGGCTCGTCGAGCTGGTCCATCCACAGGCTGAGCGTGCGCAGGGGTTGCATCGGGTTTGGCTCCACATCCGTCAAg                                                                                                            <  2:336426/149‑1 (MQ=255)
                                       gTATAGCCGGCACCCATGAACCACACGACAACCTCGAGGTCAATGCGCAGTCCCGGGCGCGCGAACAGCGGCTCGTCGAGCTCGTCCATCCACAGCCTGAGACAGCCCACGCGTTGCATCGGGATTGGCTCCACATCCGTCAAgcactg                                                                                                       >  4:85707/1‑144 (MQ=255)
                                                             acacGTCAACGTCGAGGTCATTGCGCAGTGCCGGGCGCGCGCACAGCGGCTCGTCGAGCTGGTCCATCCACAGGCTGAGCGTGCGCAGGGGTTGCATCGGGTTTGGCTCCACATCCGTCAAGGTCTGTGGGCGATGCTAGTGGTGTCAg                                                                                 >  4:276676/1‑149 (MQ=255)
                                                              cacGTCAACGTCGAGGTCATTGCGCAGTGCCGGGCGCGCGCACAGCGGCTCGTCGAGCTGGTCCATCCACAGGCTGAGCGTGCGCAGGGGTTGCATCGGGTTTGGCTCCACATCCGTCAAGGTCTGTGGGCGATGCTAGTGGTGTCAGt                                                                                <  8:121392/149‑1 (MQ=255)
                                                                cGTCAACGTCGAGGTCATTGCGCAGTGCCGGGCGCGCGCACAGCGGCTCGTCGAGCTGGTCCATCCACAGGCTGAGCGTGCGCAGGGGTTGCATCGGGTTTGGCTCCACATCCGTCAAGGTCTGTGGGCGATGc                                                                                             <  2:252465/134‑1 (MQ=255)
                                                                cGTCAACGTCGAGGTCATTGCGCAGTGCCGGGCGCGCGCACAGCGGCTCGTCGAGCTGGTCCAACCACAGGCTGAGCGTGCGCAGGGGTTGCATCGGGTTTGGCTCCACATCCGTCAAGGTCTGTGGGCGATGc                                                                                             >  1:252465/1‑134 (MQ=255)
                                                                 gTCAACGTCGAGGTCATTGCGCAGTGCCGGGCGCGCGCACAGCGGCTCGTCGAGCTGGTCCATCCACAGGCTGAGCGTGCGCAGGGGtt                                                                                                                                         <  4:31832/89‑1 (MQ=255)
                                                                 gTCAACGTCGAGGTCATTGCGCAGTGCCGGGCGCGCGCACAGCGGCTCGTCGAGCTGGTCCATCCACAGGCTGAGCGTGCGCAGGGGtt                                                                                                                                         >  3:31832/1‑89 (MQ=255)
                                                                 gTCAACGTCGAGGTCATTGCGCAGTGCCGGGCGCGCGCACAGCGGCTCGTCGAGCTGGTCCATCCACAGGCTGAGCGTGCGCAGGGGTTGCATCGGGTTTGGCTCCACATCCGTCAAGGTCTGTGGGCGATGCTAGTGGTGTCAGTggg                                                                             <  3:208660/149‑1 (MQ=255)
                                                                       gTCGAGGTCATTGCGCAGTGCCGGGCGCGCGCACAGCGGCTCGTCGAGCTGGTCCATCCACAGGCTGAGCGTGCGCAGGGGTTGCATCGGGTTTGGc                                                                                                                           <  7:347670/97‑1 (MQ=255)
                                                                       gTCGAGGTCATTGCGCAGTGCCGGGCGCGCGCACAGCGGCTCGTCGAGCTGGTCCATCCACAGGCTGAGCGTGCGCAGGGGTTGCATCGGGTTTGGc                                                                                                                           >  8:347670/1‑97 (MQ=255)
                                                                       gTCGAGGTCATTGCGCAGTGCCGGGCGCGCGCACAGCGGCTCGTCGAGCTGGTCCATCCACAGGCTGAGCGTGCGCAGGGGTTGCATCGGGTTTGGCTCCACATCCGTCAAGGTCtgt                                                                                                      <  7:29499/118‑1 (MQ=255)
                                                                       gTCGAGGTCATTGCGCAGTGCCGGGCGCGCGCACAGCGGCTCGTCGAGCTGGTCCATCCACAGGCTGAGCGTGCGCAGGGGTTGCATCGGGTTTGGCTCCACATCCGTCAAGGTCtgt                                                                                                      >  8:29499/1‑118 (MQ=255)
                                                                          gAGGTCATTGCGCAGTGCCGGGCGCGCGCACAGCGGCTCGTCGAGCTGGTCCATCCACAGGCTGAGCGTGCGCAGGGGTTGCATCGGGTTTGGCTCCACATCCGTCAAGGTCTGTGGGCGATGCTAGTGGTGTCAGTGGGCGGCTGTCt                                                                    >  5:320624/1‑149 (MQ=255)
                                                                          gAGGTCATTGCGCAGTGCCGGGCGCGCGCACAGCGGCTCGTCGAGCTGGTCCATCCACAGGCTGAGCGTGCGCAGGGGTTGCATCGGGTTTGGCTCCACATCCGTCAAGGTCTGTGGGCGATGCTAGTGGTGTCAGTGGGCGGCTGTCt                                                                    >  8:223236/1‑149 (MQ=255)
                                                                            ggTCATTGCGCAGTGCCGGGCGCGCGCACAGCGGCTCGTCGAGCTGGTCCATCCACAGGCTGAGTGTGCGCAGGGGTTGCATCGGGTTTGGCTCCACATCCGTCAAGGTCTGTGGGCGATGCTAGTGGTGTCAGTGGGCGGCTGTCtt                                                                   >  5:30099/1‑148 (MQ=255)
                                                                            ggTCATTGCGCAGTGCCGGGCGCGCGCACAGCGGCTCGTCGAGCTGGTCCATCCACAGGCTGAGTGTGCGCAGGGGTTGCATCGGGTTTGGCTCCACATCCGTCAAGGTCTGTGGGCGATGCTAGTGGTGTCAGTGGGCGGCTGTCtt                                                                   >  5:35412/1‑148 (MQ=255)
                                                                            ggTCATTGCGCAGTGCCGGGCGCGCGCACAGCGGCTCGTCGAGCTGGTCCATCCACAGGCTGAGTGTGCGCAGGGGTTGCATCGGGTTTGGCTCCACATCCGTCAAGGTCTGTGGGCGATGCTAGTGGTGTCAGTGGGCGGCTGTCtt                                                                   <  6:30099/148‑1 (MQ=255)
                                                                            ggTCATTGCGCAGTGCCGGGCGCGCGCACAGCGGCTCGTCGAGCTGGTCCATCCACAGGCTGAGTGTGCGCAGGGGTTGCATCGGGTTTGGCTCCACATCCGTCAAGGTCTGTGGGCGATGCTAGTGGTGTCAGTGGGCGGCTGTCtt                                                                   <  6:35412/148‑1 (MQ=255)
                                                                               cATTGCGCAGTGCCGGGCGCGCGCACAGCGGCTCGTCGAGCTGGTCCATCCACAGGCTGAGCGTGCGCAGGGGTTGCATCGGGTTTGGCTCCACATCCGTCTAGGGCTGTGGGCGATGCTAGTGGGGTCAGTCGGCGGCTGTCTTAcgt                                                               >  2:296053/1‑149 (MQ=255)
                                                                                    cgcAGTGCCGGGCGCGCGCACAGCGGCTCGTCGAGCTGGTCCTTCCACAGGCTGAGGGTGCGCAGGGGTTGCATCGGGTTTGGCTCCACATCCGTCAAGGTCTGTGGGCGATGGTAGTGGTGTCAGTGGGCGGCTGTCTTAcgtgcgt                                                           <  3:85707/148‑1 (MQ=255)
                                                                                           ccGGGCGCGCGCACAGCGGCTCGTCGAGCTGGTCCATCCACAGGCTGAGCGTGCGCAGGGGTTGCATCGGGTTTGGCTCCACATCCGTCAAGGTCTGTGGGCGATGCTAGTGGTGTCAGTGGGCGGCTGTCTTACGTGCGTGCACGCAg                                                   >  8:275237/1‑149 (MQ=255)
                                                                                           ccGGGCGCGCGCAAAGCGGCTCGTCGAGCTGGTCCATCCCCAGGCTGAGCGTGCGCAGGGGTTGCATCGGGTTTGGCTCCACATCCGTCAAGGTGTGTGGGCGATGCTAGTGGTGTCAGTGGGCGGCTGTCTTACGTGCGTGCACGCAg                                                   >  6:192257/1‑149 (MQ=255)
                                                                                               gcgcgcgcACAGCGGCTCGTCGAGCTGGTCCATCCACAGGCTGAGCGTGCGCAGGGGTTGCATCGGGTTTGGCTCCACATCCGTCAAGGTCTGTGGGCGATGCTAGTGGTGTCAGTGGGCGGCTGTCTTACGTGCGTGCACGCAGGGa                                                >  5:35854/1‑148 (MQ=255)
                                                                                                   gcgcACAGCGGCTCGTCGAGCTGGTCCATCCACAGGCTGAGCGTGCGCAGGGGTTGCATCGGGTTTGGCTCCACATCCGTCAAGGTCTGTGGGCGATGCTAGTGGTGTCAGTGGGCGGCTGTCTTACGTGCGTGCACGCAGGGAAAttt                                           >  8:208514/1‑149 (MQ=255)
                                                                                                     gcACAGCGGCTCGTCGAGCTGGTCCATCCANAGGCTGAGCGTGCGCAGGGGTTGCATCGGGTTTGGCTCCACATCCGTCAAGGTCTGTGGGCGATGCTAGTGGTGTCAGTGGGCGGCTGTCTTACGTGCGTGCACGCAGGGAAATTTGc                                         >  1:31997/1‑149 (MQ=255)
                                                                                                     gcACAGCGGCTCGTCGAGCTGGTCCATCCACAGGCTGAGCGTGCGCAGGGGTTGCATCGGGTTTGGCTCCACAGCCGTCAAGGTCTGTGGGCGATGCTAGTGGTGTCAGTGGGCGGCTGTCTTACGTGCGTGCCCGCAGGGAAATTTGc                                         >  6:280104/1‑149 (MQ=255)
                                                                                                       acaGCGGCTCGTCGAGCTGGTCCATCCACAGGCTGAGCGTGCGCAGGGGTTGCATCGGGTTTGGCTCCACATCCGTCAAGGTCTGTGGGCGATGCTAGTGGTGTCAGTGGGCGGCTGTCTTACGTGCGTGCACGCAGGGAAATTTGCtt                                       <  3:129583/149‑1 (MQ=255)
                                                                                                       acaGCGGCTCGTCGAGCTGGTCCATCCACAGGCTGAGCGTGCGCAGGGGTTGCATCGGGTTTGGCTCCACATCCGTCAAGGTCTGTGGGCGATGCTAGTGGTGTCAGTGGGCGGCTGTCTTACGTGCGTGCACGCAGGGAAATTTGCtt                                       <  7:105719/149‑1 (MQ=255)
                                                                                                          gCGGCTCGTCGAGCTGGTCCATCCACCGGCTGAGCGTGCGCAGGGGTTGCATCGGGTTTGGCTCCACATCCGTCAAGGTCTGTGGGCGATGCTAGTGGTGTCAGTGGGCGGCTGTCTTACGTGCGTGCACGCAGGGAAATTTGCTTCAg                                    <  6:196991/149‑1 (MQ=255)
                                                                                                                    gAGCTGGTCCATCCACAGGCTGAGCGTGCGCAGGGGTTGCATCGGGTTTGGCTCCACATCCGTCAAGGTCTGTGGGCGATGCTAGTGGTGTCAGTGGGCGGCTGTCTTACGTGCGTGCACGCAGGGAAATTTGCTTCAGGTAGGCCtt                           <  7:134760/148‑1 (MQ=255)
                                                                                                                     aGCTGGTCCATCCACAGGCTGAGCGTGCGCAGGGGTTGCATCGGGTTTGGCTCCACATCCGTCAAGGTCTGTGGGCGATGCTAGTGGTGTCAGTGGGCGGCTGTCTTACGTGCGTGCACGCAGGGAAATTTGCTTCAGGTAGGCCTTgg                         <  1:296053/149‑1 (MQ=255)
                                                                                                                          gTCCATCCACAGGCTGAGCGTGCGCAGGGGTTGCATCGGGTTTGGCTCCACATCCGTCAAGGTCTGTGGGCGATGCTAGTGGTGTCAGTGGGCGGCTGTCTTACGTGCGTGCACGCAGGGAAATTTGCTTCAGGTAGGCCTTGGGCGt                     >  4:189483/1‑148 (MQ=255)
                                                                                                                             catacaCAGCCTGAGCGTGCGCAGGGGTAGCATCGGTTTTGGCTCCACATCCGTCAAGGTCTGTGGGCGATGCTAGTGGTGTCAGTGGGCGGCTGTCTTACGTGCGTGCACCCAGGGAAATTTGCt                                        >  3:277478/1‑126 (MQ=255)
                                                                                                                                 cacaGGCTGAGCGTGCGCAGGGGTTGCATCGGGTTTGGCTCCACATCCGTCAAGGTCTGTGGGCGATGCTAGTGGTGTCAGTGGGCGGCTGTCTTACGTGCGTGCACGCAGGGAAATTTGCTTCAGGTAGGCCTTGGGCGTGAACCCg              <  3:189483/148‑1 (MQ=255)
                                                                                                                                 cacaGGCTGAGCGTGCGCAGGGGTTGCATCGGGTTTGGCTCCACATCCGTCAAGGTCTGTGGGCGATGCTAGTGGTGTCAGTGGGCGGCGGTCTTACGTGAGTGCACGCAGGGAAATTTGCTTCAGGTAGGCGTTGGGCGTGACCCCgg             >  2:35262/1‑149 (MQ=255)
                                                                                                                                  acaGGCTGAGCGTGCGCAGGGGTTGCATCGGGTTTGGCTCCACATCCGTCAAGGTCTGTGGGCGATGCTAGTGGTGTCAGTGGGCGGCTGTCTTACGTGCGTGCACGCAGGGAAATTTGCTTCAGGTAGGCCTTGGGCGTGAACCCGgt            >  2:52319/1‑149 (MQ=255)
                                                                                                                                  acaGGCTGAGCGTGCGCAGGGGTTGCATCGGGTTTGGCTCCACATCCGTCAAGGTCTGTGGGCGATGCTAGTGGTGTCAGTGGGCGGCTGTCTTACGTGCGTGCACGCAGGGAAATTTGCTTCAGGTAGGCCTTGGGCGTGAACCCGgt            >  4:21112/1‑149 (MQ=255)
                                                                                                                                  acaGGCTGAGCCTGCGCAGGGGTTGCATCGGGTTTGGCTCCACATCCGTCAAGGTCTGTGGGCGATGCTAGTGGGGTCAGTGGGCGGCGGTCTTACGTGCGTGCACGCAGGGAAATTTGCTTCAGGTAGGCCTTGGGCGTGAACCCGgt            <  8:198332/149‑1 (MQ=255)
                                                                                                                                    aGGCTGAGCGTGCGCAGGGGTTGCATCGGGTTTGGCTCCACATCCGTCAAGGTCTGTGGGCGATGCTAGTGGTGTCAGTGGGCGGCTGTCTTACGTGCGTGCACGCAGGGAAATTTGCTTCAGGTAGGCCTTGGGCGTGAACCCGgtgt          <  7:275237/149‑1 (MQ=255)
                                                                                                                                       cTGAGCGTGCGCAGGGGTTGCATCGGGTTTGGCTCCACATCCGTCAAGGTCTGTGGGCGATGCTAGTGGTGTCAGTGGGCGGCTGTCTTACGTGCGTGCACGCAGGGAAATTTGCTTCAGGTAGGCCTTGGGCGTGAACCCGGTGTGtt       >  2:226610/1‑149 (MQ=255)
                                                                                                                                          aGCGTGCGCAGGGGTTGCATCGGGTTTGGCTCCACATCCGTCAAGGTCTGTGGGCGATGCTAGTGGTGTCAGTGGGCGGCTGTCTTACGTGCGTGCACGCAGGGAAATTTGCTTCAGGTAGGCCTTGGGCGTGAACCCGGTGTGTTcgc    <  7:223236/149‑1 (MQ=255)
                                                                                                                                           gCGTGCGCAGGGGTTGCATCGGGTTTGGCTCCACATCCGTCAAGGTCTGTGGGCGATGCTAGGGGGGTCAGTGGGCGGCTGTCTTACGTGCGTGCACGCAGGGAAATTTGCTTCAGGTAGGCCTTGGGCGTGAACCCGGTGTGTTcgcg   <  5:280104/149‑1 (MQ=255)
                                                                                                                                            cGTGCGCAGGGGTTGCATCGGGTTTGGCTCCACATCCGTCAAGGTCTGTGGGCGATGCTAGTGGTGTCAGTGGGCGGCTGTCTTACGTGCGTGCACGCAGGGAAATTTGCTTCAGGTAGGCCTTGGGCGTGAACCCGGTGTGTTCGCgg  <  7:208514/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                    
GGGGCCTGGCGCTTGAGGTAGTACGCGGTCCAAAGGCCGGTATAGCCGGCACCGATGATGCACACGTCAACGTCGAGGTCATTGCGCAGTGCCGGGCGCGCGCACAGCGGCTCGTCGAGCTGGTCCATCCACAGGCTGAGCGTGCGCAGGGGTTGCATCGGGTTTGGCTCCACATCCGTCAAGGTCTGTGGGCGATGCTAGTGGTGTCAGTGGGCGGCTGTCTTACGTGCGTGCACGCAGGGAAATTTGCTTCAGGTAGGCCTTGGGCGTGAACCCGGTGTGTTCGCGG  >  NC_002947/1016059‑1016347

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