Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,822,424 C→T 8.3% G128G (GGC→GGT)  PP_5099 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,822,4240CT 8.3% 137.4 / 2.5 48G128G (GGC→GGT) PP_5099hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (23/21);  new base T (2/2);  total (25/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 = 9.10e-01

TCAGCGGTAAGGTCAGTGACCTGACCGGGCGCAGCAGCCCGCTGTCGTTCAAGCAGATTACCGAACAGGGCGCGGTGTATTACATCGCCCAGTTCAAGATCGATCAGCCAGAAACCGTCACCTTCGACTTGAACATCGAAACCGGCGGCATCAGCAACTCTCTCAGCTTCAACCAGGAAGTGTTCCCAGGCGAATGATGAATTTCCAGCAACTCGTATTGGCCAGCCACAATGCCGGCAAACTGAAGGAACTCCAGGCCATGCTCGGCCATTCCGTGCAGCTGCACTCGATTG  >  NC_002947/5822279‑5822571
                                                                                                                                                 |                                                                                                                                                   
tCAGCGGTAAGGTCAGTGACCTGACCGGGCGCAGCAGCCCGCTGTCGTTCAAGCAGATTACCGAACAGGGCGCGGTGTATTACATCGCCCAGTTCAAGATCGATCAGCCAGAAACCGTCACCTTCGACTTGATCATCGAAACcggcgg                                                                                                                                                   >  8:240308/1‑148 (MQ=255)
     ggTAAGGTCAGTGACCTGACCGGGCGCAGCAGCCCGCTGTCGTTCAAGCAGATTACCGAACAGGGCGCGGTGTATTACATCGCCCAGTTCAAGATCGATCAGCCAGAAACCGTCACCTTCGACTTGAACATCGAAACCGGCGGCATCAg                                                                                                                                             >  7:88098/1‑149 (MQ=255)
     ggTAAGGTCAGTGACCTGACCGGGCGCAGCAGCCCGCTGTCGTTCAAGCAGATTACCGAACAGGGCGCGGTGTATTACATCGCCCAGTTCAAGATCGATCAGCCAGAAACCGTCACCTTCGACTTGAACATCGAAACCGGCGGCATCAg                                                                                                                                             <  4:226421/149‑1 (MQ=255)
         aGGTCAGTGACCTGACCGGGCGCAGCAGCCCGCTGTCGTTCAAGCAGATTACCGAACAGGGCGCGGTGTATTACATCGCCCAGTTCAAGATCGATCAGCCAGAAACCGTCACCTTCGACTTGAACATCGAAACCGGCGGCATCAGCaa                                                                                                                                          <  3:302008/148‑1 (MQ=255)
           gTCAGTGACCTGACCGGGCGCAGCAGCCCGCTGTCGTTCAAGCAGATTACCGAACAGGGCGCGGTGTATTACATCGCCCAGTTCAAGATCGATCAGCCAGAAACCGTCACCTTCGACTTGAACATCGAAACcggc                                                                                                                                                     <  8:81016/135‑1 (MQ=255)
           gTCAGTGACCTGACCGGGCGCAGCAGCCCGCTGTCGTTCAAGCAGATTACCGAACAGGGCGCGGTGTATTACATCGCCCAGTTCAAGATCGATCAGCCAGAAACCGTCACCTTCGACTTGAACATCGAAACcggc                                                                                                                                                     >  7:81016/1‑135 (MQ=255)
              aGTGACCTGACCGGGCGCAGCAGCCCGCTGTCGTTCAAGCAGATTACCGAACAGGGCGCGGTGTATTACATCGCCCAGTTCAAGATCGATCAGCCAGAAACCGTCACCTTCGACTTGAACATCGAAACCGGCGGCATCAGCAActctct                                                                                                                                    <  2:100217/149‑1 (MQ=255)
               gTGACCTGACCGGGCGCAGCAGCCCGCTGTCGTTCAAGCAGATTACAGAACAAGGCGCGGTGTATTACATCGCCCAGTTCAAGATCGATCAGCCAGAAACCGTCACCTTCGACTTGAACATCGAAACCGGCGGCATCAGCAActctct                                                                                                                                    >  1:238385/1‑148 (MQ=255)
                    ctgaccGGGCGCAGCAGCCCGCTGTCGTTCAAGCAGATTACCGAACAGGGCGCGGTGTATTACATCGCCCAGTTCAAGATCGATCAGCCAGAAACCGTCACCTTCGACTTGAACATCGAAACCGGCGGCATCAGCAACTCTCTCAGCt                                                                                                                               >  6:18574/1‑148 (MQ=255)
                         cGGGCGCAGCAGCCCGCTGTCGTTCAAGCAGATTACCGAACAGGGCGCGGTGTATTACATCGCCCAGTTCAAGATCGATCAGCCAGAAACCGTCACCTTCGACTTGAACATCGAAACCGGCGGCATCAGCAACTCTCTCAGCTTCAAcc                                                                                                                         <  8:88098/149‑1 (MQ=255)
                              gcagcaGCCCGCTGTCGTTCAAGCAGATTACCGAACAGGGCGCGGTGTATTACATCGCCCAGTTCAAGATCGATCAGCCAGAAACCGTCACCTTCAACTTGAACATCGAAACCGGCGGCATCAGCAACTCTCTCAGCTACAACCAgg                                                                                                                      >  6:195841/1‑147 (MQ=255)
                                agcagcCCGCTGTCGTTCAAGCAGATTACCGAACATGGCGCGGTGTATTACATCGCCCAGTTCAAGATCGATCAGCCAGAAACCGTCACCTTCGACTTGAACATCGAAACCGGCGGCATCAGCAACTCTCTCAGCTTCAACCAGTAAgt                                                                                                                  <  2:238385/149‑1 (MQ=255)
                                  cagcCCGCTGTCGGTCAAGCAGATTACCGAACAGGGCGCGGTGTATTACATCGCCCAGTTCAAGATCGATCAGCCAGAAACCGTCACCTTCGACTTGAACATCGAAACCGGCGGCATCAGCAACTCTCTCAGCTTCAACCAGGAAgtgt                                                                                                                <  8:280007/149‑1 (MQ=255)
                                     cccGCTGTCGTTCAAGCAGATTACCGAACAGGGCGCGGTGTATTACATCGCCCAGTTCAAGATCGATCAGCCAGAAACCGTCACCTTCGACTTGAACATCGAAACCGGCGGCATCAGCAACTCTCTCAGCTTCAACCAGGAAGTGTTcc                                                                                                             >  4:46613/1‑149 (MQ=255)
                                      ccGCTGTCGTTCAAGCAGATTACCGAACAGGGCGCGGTGTATTACATCGCCCAGTTCAAGATCGATCAGCCAGAAACCGTCACCTTCGACTTGAACATCGAAACCGGCGGCATCAGCAACTCTCTCAGCTTCAACCAGGAAGTGTTcc                                                                                                             <  3:46613/148‑1 (MQ=255)
                                       cGCTGTCGTTCAAGCAGATTACCGAACAGGGCGCGGTGTATTACATCGCCCAGTTCAAGATCGATCAGCCAGAAACCGTCACCTTCGACTTGAACATCGAAACCGGCGGCATCAGCAACTCTCTCAGCTTCAACCAGGAAGTGTTcc                                                                                                             >  2:190746/1‑147 (MQ=255)
                                                  aaGCAGATTACCGAACAGGGCGCGGTGTATTACATCGCCCAGTTCAAGATCGATCAGCCAGAAACCGTCACCTTCGACTTGAACATCGAAACCGGCGGCATCAGCAACTCTCTCAGCTTCAACCAGGAAGTGTTCCCAGGCGAatgat                                                                                                 <  1:84926/148‑1 (MQ=255)
                                                    gCAGATTACCGAACAGGGCGCGGTGTATTACATCGCCCAGTTCAAGATCGATCAGCCAGAAACCGTCACCTTCGACTTGATCATCGAAACCGGCGGCATCAGCAACTCTCTCAGCTTCAACCAGGAAGTGTTCCCAGGCGAATGATGaa                                                                                              <  7:240308/149‑1 (MQ=255)
                                                             cGAACAGGGCGCGGTGTATTACATCGCCCAGTTCAAGATCGATCAGCCAGAAACCGTCACCTTCGACTTGAACATCGAAACCGGCGGCATCAGCAACTCTCTCAGCTTCAACCAGGAAGTGTTCCCAGGCGAATGATGAATTTCCAGc                                                                                      >  7:7642/1‑148 (MQ=255)
                                                               aaCAGGGCGCGGTGTATTACATCGCCCAGTTCAAGATCGATCAGCCAGAAACCGTCACCTTCGACTTGAACATCGAAACCGGCGGCATCAGCAACTCTCTCAGCTTCAAc                                                                                                                          >  3:244600/1‑110 (MQ=255)
                                                               aaCAGGGCGCGGTGTATTACATCGCCCAGTTCAAGATCGATCAGCCAGAAACCGTCACCTTCGACTTGAACATCGAAACCGGCGGCATCAGCAACTCTCTCAGCTTCAAc                                                                                                                          <  4:244600/110‑1 (MQ=255)
                                                                   gggCGCGGTGTATTACATCGCCCAGTTCAAGATCGATCAGCCAGAAACCGTCACCTTCGACTTGAACATCGAAACCGGCGGCATCAGCAACTCTCTCAGCTTCAACCAGGAAGTGTTCCCAGGCGAATGATGAATTTCCAGCAACTCGt                                                                               <  4:306993/149‑1 (MQ=255)
                                                                         ggTGTATTACATCGCCCAGTTCAAGATCGATCAGCCAGAAACCGTCACCTTCGACTTGAACATCGAAACCGGCGGCATCAGCAACTCTCTCAGCTTCAACCAGGAAGTGTTCCCAGGCGAATGATGAATTTCCAGCAACTCGTATTGgc                                                                         >  6:114528/1‑149 (MQ=255)
                                                                               ttACATCGCCCAGTTCAAGATCGATCAGCCAGAAACCGTCACCTTCGACTTGAACATCGAAACCGGCGGCATCAGCAACTCTCTCAGCTTCAACCAGGAAGTGTTCCCAGGCGAATGATGAATTTCCAGCAACTCGTATTGgccagcca                                                                   >  6:52909/1‑149 (MQ=255)
                                                                                       cccAGTTCAAGATCGATCAGCCAGAAACCGTCACCTTCGACTTGAACATCGAAACCGGCGGCATCAGCAACTCTCTCAGCTTCAACCAGGAAGTGTTCCCAGGCGAATGATGAATTTCCAGCAACTCGTATTGGCCAGCCACAATGCCg                                                           <  1:190746/149‑1 (MQ=255)
                                                                                        ccAGTTCAAGATCGATCAGCCAGAAACCGTCACCTTCGACTTGAACATCGAAACCGGCGGCATCAGCAACTCTCTCAGCTTCAACCAGGAAGTGTTCCCAGGCGAATGATGAATTTCCAGCAACTCGTATTGGCCAGCCACAATGCCgg                                                          >  3:37816/1‑149 (MQ=255)
                                                                                         cAGTTCAAGATCGATCAGCCAGAAACCGTCACCTTCGACTTGAACATCGAAACCGGCGGCATCAGCAACTCTCTCAGCTTCAACCAGGAAGTGTTCCCAGGCGAatgat                                                                                                 >  1:44545/1‑109 (MQ=255)
                                                                                         cAGTTCAAGATCGATCAGCCAGAAACCGTCACCTTCGACTTGAACATCGAAACCGGCGGCATCAGCAACTCTCTCAGCTTCAACCAGGAAGTGTTCCCAGGCGAatgat                                                                                                 <  2:44545/109‑1 (MQ=255)
                                                                                                aGATCGATCAGCCAGAAACCGTCACCTTCGACTTGAACATCGAAACCGCCGGCATCAGCAACTCTCTCAGCTTCAACCAGGAAGTGTTCCCAGGCGAATGATGAATTTCCAGCAACTCGTATTGGCCAGCCACAATGCCGGCAAACTg                                                   >  4:16111/1‑148 (MQ=255)
                                                                                                aGATCGATCAGCCAGAAACCGTCACCTTCAACTTGAACATCGAAACCGCCGTCATCAGCAACTCTCTCAGCTTCAACCAGGAAGTGTTCCCAGGCGAATGATGAATTTCCAGCAACTCGTATTGGCCAGCCACAATGCCGGCAAACTGa                                                  <  5:195841/149‑1 (MQ=255)
                                                                                                 gatcgatcAGCCAGAAACCGTCACCTTCGACTTGAACATCGAAACCGGCGGCATCAGCAACTCTCTCAGCTTCAACCAGGAAGTGTTCCCAGGCGAATGATGAATTTCCAGCAACTCGTATTGGCCAGCCACAATGCCGGCAAACTGaa                                                 <  4:37816/149‑1 (MQ=255)
                                                                                                   tcgatcAGCCAGAAACCGTCACCTTCGACTTGAACATCGAAACCGGCGGCATCAGCAACTCTCTCAGCTTCAACCAGGAAGTGTTCCCAGGCGAATGATGAATTTCCAGCAACTCGTATTGGCCAGCCACAATGCCGGCAAACTGAAgg                                               <  6:80822/149‑1 (MQ=255)
                                                                                                        cAGCCAGAAACCGTCACCTTCGACTTGAACATCGAAACCGGCGGCATCAGCAACTCTCTCAGCTTCAACCAGGAAGTGTTCCCAGGCGAATGATGAATTTCCAGCAACTCGTATTGGCCAGCCACAATGCCGGCAAACTGAAGGAACTc                                          >  1:113349/1‑149 (MQ=255)
                                                                                                           ccAGAAACCGTCACCTTCGACTTGAACATCGAAACCGGCGGCATCAGCAACTCTCTCAGCTTCAACCAGGAAGTGTTCCCAGGCGAATGATGAATTTCCCGCAACTCGTATTGGCCCGCCACAATGCGGGCAAAGTGAGGGAAGTCCAg                                       >  3:220835/1‑149 (MQ=255)
                                                                                                                   cGTCACCTTCGACTTGAACATCGAAACCGGCGGCATCAGCAACTCTCTCAGCTTCAACCAGGAAGTGTTCCCAGGCGAATGATGAATTTCCAGCAACTCGTATTGGCCAGCCACAATg                                                              <  8:24700/118‑1 (MQ=255)
                                                                                                                   cGTCACCTTCGACTTGAACATCGAAACCGGCGGCATCAGCAACTCTCTCAGCTTCAACCAGGAAGTGTTCCCAGGCGAATGATGAATTTCCAGCAACTCGTATTGGCCAGCCACAATg                                                              >  7:24700/1‑118 (MQ=255)
                                                                                                                         cTTCGACTTGAACATCGAAACCGGCGGCATCAGCAACTCTCTCAGCTTCAACCAGGAAGTGTTCCCAGGCGAATGATGAATTTCCAGCAACTCGTATTGGCCAGCCACAATGCCGGCAAACTGAAGGAACTCCAGGCCATGCTCGGCCa                         <  5:114528/149‑1 (MQ=255)
                                                                                                                                ttGAACATCGAAACCGGCGGCATCAGCAACTCTCTCAGCTTCAACCAGGAAGTGTTCCCAGGCGAATGATGAATTTCCAGCAACTCGTATTGGCCAGCCACAATGCCGGCAAACTGAAGGAACTCCAGGCCATGCTCGGCCATTCCGTg                  >  7:178378/1‑149 (MQ=255)
                                                                                                                                   aaCATCGAAACCGGCGGCATCAGCAACTCTCTCAGCTTCAACCAGGAAGTGTTCCCAGGCGAATGATGAATTTCCAGCAACTCGTATTGGCCAGCCACAATGCCGGCAAACTGAAGGAACTCCAGGCCATGCTCGGCCATTCCGTGCAg               >  4:41962/1‑149 (MQ=255)
                                                                                                                                      aTCGAAACCGGCGGCATCAGCAACTCTCTCAGCTTCAACCAGGAAGTGTTCCCAGGCGAATGATGAATTTCCAGCAACTCGTATTGGCCAGCCACAATGCCGGCAAACTGAAGGAACTCCAg                                       >  3:21460/1‑122 (MQ=255)
                                                                                                                                      aTCGAAACCGGCGGCATCAGCAACTCTCTCAGCTTCAACCAGGAAGTGTTCCCAGGCGAATGATGAATTTCCAGCAACTCGTATTGGCCAGCCACAATGCCGGCAAACTGAAGGAACTCCAg                                       <  4:21460/122‑1 (MQ=255)
                                                                                                                                      aTCGAAACCGGCGGCATCAGCAACTCTCTCAGCTTCAACCAGGAAGTGTTCCCAGGCGAATGATGAATTTCCAGCAACTCGTATTGGCCAGCCACAATGCCGGCAAACTGAAGGAACTCCAGGCCATGCTCGGCCATTCCGTGCAGCTg            >  5:228020/1‑149 (MQ=255)
                                                                                                                                        cGAAACCGGCGGCATCAGCAACTCTCTCAGCTTCAACCAGGAAGTGTTCCCAGGCGAATGATGAATTTCCAGCAACTCGTATTGGCCAGCCACAATGCCGGCAAACTGAAGGAACTCCAGGCCATGCTCGGCCATTcc                     >  3:62908/1‑138 (MQ=255)
                                                                                                                                        cGAAACCGGCGGCATCAGCAACTCTCTCAGCTTCAACCAGGAAGTGTTCCCAGGCGAATGATGAATTTCCAGCAACTCGTATTGGCCAGCCACAATGCCGGCAAACTGAAGGAACTCCAGGCCATGCTCGGCCATTcc                     <  4:62908/138‑1 (MQ=255)
                                                                                                                                         gAAACCGGCGGCATCAGCAACTCTCTCAGCTTCAACCAGGAAGTGTTCCCAGGCGAATGATGAATTTCCAGCAACTCGTATTGGCCAGCCACAATGCCGGCAAACTGaa                                                 <  1:79698/109‑1 (MQ=255)
                                                                                                                                         gAAACCGGCGGCATCAGCAACTCTCTCAGCTTCAACCAGGAAGTGTTCCCAGGCGAATGATGAATTTCCAGCAACTCGTATTGGCCAGCCACAATGCCGGCAAACTGaa                                                 >  2:79698/1‑109 (MQ=255)
                                                                                                                                          aaaCCGGTGGCATCAGCAACTCTCTCAGCTTCAACCAGGAAGTGTTCCCAGGCGAATGATGAATTTCCAGCAACTCGTATTGGCCAGCCACAATGCCGGCAAACTGAAg                                                <  7:121084/109‑1 (MQ=255)
                                                                                                                                          aaaCCGGTGGCATCAGCAACTCTCTCAGCTTCAACCAGGAAGTGTTCCCAGGCGAATGATGAATTTCCAGCAACTCGTATTGGCCAGCCACAATGCCGGCAAACTGAAg                                                >  6:211522/1‑109 (MQ=255)
                                                                                                                                          aaaCCGGTGGCATCAGCAACTCTCTCAGCTTCAACCAGGAAGTGTTCCCAGGCGAATGATGAATTTCCAGCAACTCGTATTGGCCAGCCACAATGCCGGCAAACTGAAg                                                >  8:121084/1‑109 (MQ=255)
                                                                                                                                          aaaCCGGTGGCATCAGCAACTCTCTCAGCTTCAACCAGGAAGTGTTCCCAGGCGAATGATGAATTTCCAGCAACTCGTATTGGCCAGCCACAATGCCGGCAAACTGAAg                                                <  5:211522/109‑1 (MQ=255)
                                                                                                                                           aaCCGGCGGCATCAGCAACTCTCTCAGCTTCAACCAGGAAGTGTTCCCAGGCGAATGATGAATTTCCAGCAACTCGTATTGGCCAGCCACAATGCCGGCAAACTGAAGGAACTCCAGGCCATGCTCGGCCATTCCGTGCAGCTGCACTc       >  2:211535/1‑149 (MQ=255)
                                                                                                                                               ggcggcATCAGCAACTCTCTCAGCTTCAACCAGGAAGTGTTCCCAGGCGAATGATGAATTTCCAGCAACTCGTATTGGCCAGCCACAATGCCGGCAAACTGAAGGAACTCCAGGCCATGCTCGGCCATTCCGTGCAGCTGCACTCGAtt   <  5:52909/149‑1 (MQ=255)
                                                                                                                                                gcggcATCAGCAACTCTCTCAGCTTCAACCAGGAAGTGTTCCCAGGCGAATGATGAATTTCCAGCAACTCGTATTGGCCAGCCACAATGCCGGCAAACTGAAGGAACTCCAGGCCATGCTCGGCCATTCCGTGCAGCTGCACTCGATTg  >  6:131768/1‑149 (MQ=255)
                                                                                                                                                gcggcATCAGCAACTCTCTCAGCTTCAACCAGGAAGTGTTCCCAGGCGAATGATGAATTTCCAGCAACTCGTATTGGCCAGCCACAATGCCGGCAAACTGAAGGAACTCCAGGCCATGCTCGGCCATTCCGTGCAGCTGCACTCGATTg  <  2:113349/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TCAGCGGTAAGGTCAGTGACCTGACCGGGCGCAGCAGCCCGCTGTCGTTCAAGCAGATTACCGAACAGGGCGCGGTGTATTACATCGCCCAGTTCAAGATCGATCAGCCAGAAACCGTCACCTTCGACTTGAACATCGAAACCGGCGGCATCAGCAACTCTCTCAGCTTCAACCAGGAAGTGTTCCCAGGCGAATGATGAATTTCCAGCAACTCGTATTGGCCAGCCACAATGCCGGCAAACTGAAGGAACTCCAGGCCATGCTCGGCCATTCCGTGCAGCTGCACTCGATTG  >  NC_002947/5822279‑5822571

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