Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,785,761 C→T 7.5% A319V (GCC→GTC)  PP_3348 → GGDEF domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,785,7610CT 7.5% 150.7 / 2.2 53A319V (GCC→GTC) PP_3348GGDEF domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (22/27);  new base T (2/2);  total (24/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.57e-01

TGCAGCGCATCGGTGCAGCGATCAGTGCCAGCCTGCGTGACAACGACATTGCCGGCAGGCTGGGCGGCGAGGAGTTCGCCATCCTGCTGCCCAATACGCCGCCAGACACGGCGATTGAGGTTGCCGAGCGCCTGCGGGCCCTGATAGCCAAACTGCGTGTGAATGATGAACACCCGGTCACCGCAAGCATAGGTTTAGCGTTCGCGAAGGCCCCTGCCACTGAGCTCGGCGCCTTGCTCAGCAGTGCAGACAGCGCCATGTACCTTTCCAAGGTATCGGGCCGTAACCGGGTCAC  >  NC_002947/3785614‑3785908
                                                                                                                                                   |                                                                                                                                                   
tgcAGCGCATCGGTGCAGCGATCAGTGCCAGCCTGCGTGACAACGACATTGCCGGCAGGCTGGGCGGCGAGGAGTTCGCCATCCTGCTGCCCAATACGCCGCCAGACACGGCGATTGAGGTTGCCGAGCGCCTGCGGGCCCTGATAGc                                                                                                                                                     >  4:127864/1‑148 (MQ=255)
 gcAGCGCATCGGTGCAGCGATCAGTGCCAGCCTGCGTGACAACGACATTGCCGGCAGGCTGGGCGGCGAGGAGTTCGCCATCCTGCTGCCCAATACGCCGCCAGACACGGCGATTGAGGTTGCCGAGCGCCTGCGGGCCCTGATAGcc                                                                                                                                                    >  5:255859/1‑148 (MQ=255)
           ggTGCAGCGATCAGTGCCAGCCTGCGTGACAACGACATTGCCGGCAGGCTGGGCGGCGAGGAGTTCGCCATCCTGCTGCCCAATACGCCGCCAGACACGGCGATTGAGGTTGCCGAGCGCCTGCGGGCCCTGATAGCCAAACTGCgtg                                                                                                                                          <  8:201858/148‑1 (MQ=255)
            ctGCAGCGATCAGTGCCAGCCTGCGTGACAACGACATTGCCGGCAGGCTGGGCGGCGAGGAGTTCGCCATCCTGCTGCCCAATACGCCGCCAGACACGGCGATTGAGGTTGCCGAGCGCCTGCGGGCCCTGATAGCCAAACTGCgtgt                                                                                                                                         <  1:200443/147‑1 (MQ=255)
               cAGCGATCAGTGCCAGCCTGCGTGACAACGACATTGCCGGCAGGCTGGGCGGCGAGGAGTTCGCCATCCTGCTGCCCAATACGCCGCCAGACACGGCGATTGAGGTTGCCGAGCGCCTGCGGGCCCTGATAGCCAAACTGCgtgt                                                                                                                                         >  2:200443/1‑145 (MQ=255)
                aGCGATCAGTGCCAGCCTGCGTGACAACGACATTGCCGGCAGGCTGGGAGGCGAGGAGTTCGCCATCCTGCTGCCCAATACGCCGCCAGACACGGCGATTGAGGTTGCCGAGCGCCTGCGGGCCCTGATAGTCAAACTGCGTGTGAat                                                                                                                                     >  4:80962/1‑148 (MQ=255)
                aGCGATCAGTGCCAGCCTGCGTGACAACGACATTGCCGGCAGGCTGGGAGGCGAGGAGTTCGCCATCCTGCTGCCCAATACGCCGCCAGACACGGCGATTGAGGTTGCCGAGCGCCTGCGGGCCCTGATAGTCAAACTGCGTGTGAat                                                                                                                                     >  4:85673/1‑148 (MQ=255)
                 gCGATCAGTGCCAGCCTGCGTGACAACGACATTGCCGGCAGGCTGGGCGGCGAGGAGTTCGCCATCCTGCTGCCCAATACGCCGCCAGACACGGCGATTGAGGTTGCCGAGCGCCTGCGGGCCCTGATAGCCAAACTGCGTGTGAatga                                                                                                                                   <  5:212524/149‑1 (MQ=255)
                  cGATCAGTGCCAGCCTGCGTGACAACGACATTGCCGGCAGGCTGGGCGGCGAGGAGTTCGCCATCCTACTGCCCAATACGCCGCCAGACACGGCGATTGAGGTTGCCGAGCGCCTGCGGGCCCTGATAGCCAAACTGCGTGTGAatgat                                                                                                                                  <  3:278501/149‑1 (MQ=255)
                       aGTGCCAGCCTGCGTGACAACGACATTGCCGGCAGGCTGGGCGGCGAGGAGTTCGCCATCCTGCTGCCCAATACGCCGCCAGACACGGCGATTGAGGTTGCCGAGCGCCTGCGGGCCCTGATAGCCAAACTGCGTGTGAATGATGAaca                                                                                                                             >  8:200561/1‑149 (MQ=255)
                         tGCCAGCCTGCGTGACAACGACATTGCCGGCAGGCTGGGCGGCGAGGAGTTCGCCATCCTGCTGCCCAATACGCCGCCAGACACGGCGATTGAGGTTGCCGAGCGCCTGCGGGCCCTGATAGCCAAACTGCgt                                                                                                                                           <  6:38265/133‑1 (MQ=255)
                         tGCCAGCCTGCGTGACAACGACATTGCCGGCAGGCTGGGCGGCGAGGAGTTCGCCATCCTGCTGCCCAATACGCCGCCAGACACGGCGATTGAGGTTGCCGAGCGCCTGCGGGCCCTGATAGCCAAACTGCgt                                                                                                                                           >  5:38265/1‑133 (MQ=255)
                         tGCCAGCCTGCGTGACAACGACATTGCCGGCAGGCTGGGCGGCGAGGAGTTCGCCATCCTGCTGCCCAATACGCCGCCAGACACGGCGATTGAGGTTGCCGAGCGCCTGCGGGCCCTGATAGCCAAACTGCGTGTGAAtgatga                                                                                                                                <  7:12622/144‑1 (MQ=255)
                         tGCCAGCCTGCGTGACAACGACATTGCCGGCAGGCTGGGCGGCGAGGAGTTCGCCATCCTGCTGCCCAATACGCCGCCAGACACGGCGATTGAGGTTGCCGAGCGCCTGCGGGCCCTGATAGCCAAACTGCGTGTGAAtgatga                                                                                                                                >  8:12622/1‑144 (MQ=255)
                                      gACAACGACATTGCCGGCAGGCTGGGTGGCGAGGAGTTCGCCATCCTGCTGCCCAATACGCCGCCAGACACGGCGATTGAGGTTGCCGAGCGCCTGCGGGCCCTGATAGTCAAACTGCGTGTGAATGATGAACACCCGGTCACCGCAAg                                                                                                              <  3:85673/149‑1 (MQ=255)
                                       aCAACGACATTGCCGGCAGGCTGGGGGGCGAGGAGTTCGCCATCCTGCTGCCCAATACGCCGCCAGACACGGCGATTGAGGTTGCCGAGCGCCTGCGGGCCCTGATAGTCAAACTGCGTGTGAATGATGAACACCCGGTCACCGCAAg                                                                                                              <  3:80962/148‑1 (MQ=255)
                                        cAACGACATTGCCGGCAGGCTGGGCGGCGAGGAGTTCGCCATCCTGCTGCCCAATACGCCGCCAGACACGGCGATTGAGGTTGCCGAGCGCCTGCGGGCCCTGATAGCCAAACTGCGTGTGAATGATGAacac                                                                                                                            <  4:182255/133‑1 (MQ=255)
                                        cAACGACATTGCCGGCAGGCTGGGCGGCGAGGAGTTCGCCATCCTGCTGCCCAATACGCCGCCAGACACGGCGATTGAGGTTGCCGAGCGCCTGCGGGCCCTGATAGCCAAACTGCGTGTGAATGATGAacac                                                                                                                            >  3:182255/1‑133 (MQ=255)
                                         aaCGACATTGCCGGCAGGCTGGGCGGCGAGGAGTTCGCCATCCTGCTGCCCAATACGCCGCCAGACACGGCGATTGAGGTTGCCGAGCGCCTGCGGGCCCTGATAGCCAAACTGAGTGTGAATGATGAACACCCGGTCACCGCAAGCa                                                                                                            >  1:98086/1‑148 (MQ=255)
                                          aCGACATTGCCGGCAGGCTGGGCGGCGAGGAGTTCGCCATCCTGCTGCCCAATACGCCGCCAGACACGGCGATTGAGGTTGCCGAGCGCCTGCGGGCCCTGATAGCCAAACTGCGTGTGAATGATGAACACCCGGTCACCGCAAGCa                                                                                                            >  4:151205/1‑147 (MQ=255)
                                              cATTGCCGGCAGGCTGGGCGGCGAGGAGTTCGCCATCCTGCTGCCCAATACGCCGCCAGACACGGCGATTGAGGTTGCCGAGCGCCTGCGGGCCCTGATAGCCAAACTGCGTGTGAATGATGAACACCCGGTCACCGCAAGCATAGGtt                                                                                                      >  8:135539/1‑149 (MQ=255)
                                                      gCAGGCTGGGCGGCGAGGAGTTCGCCATCCTGCTGCCCAATACGCCGCCAGACACGGCGATTGAGGTTGCCGAGCGCCTGCGGGCCCTGATAGCCAAACTGCGTGTGAATGATGAACACCCGGTCACCGCAAGCATAGGTTTAGCGTTc                                                                                              <  4:253687/149‑1 (MQ=255)
                                                              ggcggcGAGGAGTTCGCCATCCTGCTGCCCAATACGCCGCCAGACACGGCGATTGAGGTTGCCGAGCGCCTGCGGGCCCTGATAGCCAAACTGCGTGTGAATGATGAACACCCGGTCACCGCAAGCATAGGTTTAGCGTTCGCGAAGGc                                                                                      <  7:146490/149‑1 (MQ=255)
                                                                    gaggagTTCGCCATCCTGCTGCCCAATACGCCGCCAGACACGGCGATTGAGGTTGCCGAGCGCCTGCGGGCCCTGATAGCCAAACTGCGTGTGAATGATGAACACCCGGTCACCGCAAGCATAGGTTTAGCGTTCGCGAAGGCCCCTGc                                                                                >  5:167768/1‑149 (MQ=255)
                                                                       gagTTCGCCATCCTGCTGCCCAATACGCCGCCAGACACGGCGATTGAGGTTGCCGAGCGCCTGCGGGCCCTGATAGCCAAACTGCGTGTGAATGATGAACACCCGGTCACCGCAAGCATAGGTTTAGCGTTCGCGAAGGCCCCTGCCAc                                                                             <  7:200561/149‑1 (MQ=255)
                                                                          ttCGCCATCCTGCTGCCCAATACGCCGCCAGACACGGCGATTGAGGTTGCCGAGCGCCTGCGGGCCCTGGTAGCCAAACTGCGTGTGAATGATGAACACCCGGTCACCGCAAGCATAGGTTTAGCGTTCGCGAAGGCCCCTGCCACt                                                                            >  2:327977/1‑147 (MQ=255)
                                                                          ttCGCCATCCTGCTGCCCAATACGCCGCCAGACACGGCGATTGAGGTTGCCGAGCGCCTGCGGGCCCTGATAGCCAAACTGCGTGTGAATGATGAACACCCGGTCACCGCAAGCATAGGTTTAGCGTTCGCGAAGGCCCCTGCCACTGa                                                                          <  3:151205/149‑1 (MQ=255)
                                                                            cGCCATCCTGCTGCCCAATACGCCGCCAGACACGGCGATTGAGGTTGCCGAGCGCCTGCGGGCCCTGATAGCCAAACTGCGTGTGAATGATGAACACCCGGTCACCGCAAGCATAGGTTTAGCGTTCGCGAAGGCCCCTGCCACTGAg                                                                         <  1:130829/148‑1 (MQ=255)
                                                                            cGCCATCCTGCTGCCCAATACGCCGCCAGACAAGGCGATTGAGGTTGCCGAGCGCCTGCGGGCCCTGATAGCCAAACTGCGTGTGAATGATGAACACCCGGTCACCGCAAGCATAGGTTTAGCGTTCGCGAAGGCCCCTGCCACTGAg                                                                         >  6:259926/1‑148 (MQ=255)
                                                                              ccATCCTGCTGCCCAATACGCCGCCAGACACGGCGATTGAGGTTGCCGAGCGCCTGCGGGCCCTGATAGCCAAACTGCGTGTGAATGATGAACACCCGGTCACCGCAAGCATAGGTTTAGCGTTCGCGAAGGCCCCTGCCACTGAGCTc                                                                      >  6:286808/1‑149 (MQ=255)
                                                                                  ccTGCTGCCCAATACGCCGCCAGTCACGGCGATTGAGGTTGCCGAGCGCCTGCGGGCCCTGATAGCCAAACTGCGTGTGAATGATGAACACCCGGTCACCGCAAGCATAGGTTTAGCGTTCGCGAAGGCCCCt                                                                                  <  6:40239/133‑1 (MQ=255)
                                                                                   ctgctgCCCAATACGCCGCCAGACACGGCGATTGAGGTTGCCGAGCGCCTGCGGGCCCTGATAGCCAAACTGCGTGTGAATGATGAACACCCGGTCACCGCAAGCATAGGTTTAGCGTTCGCGAAGGCCCCTGCCACTGAGCTCGgcgc                                                                 >  7:94084/1‑149 (MQ=255)
                                                                                   ctgctgCCCAATACGCCGCCAGACACGGCGATTGAGGTTGCCGAGCGCCTGCGGGCCCTGATAGCCAAACTGCGTGTGAATGATGAACACCCGGTCACCGCAAGCATAGGTTTAGCGTTCGCGAAGGCCCCTGCCACTGAGCTCGgcg                                                                  <  7:135539/148‑1 (MQ=255)
                                                                                         cccAATACGCCGCCAGACACGGCGATTGAGGTTGCCGAGCGCCTGCGGGCCCTGATAGCCAAACTGCGTGTGAATGATGAACACCCGGTCACCGCAAGCATAGGTTTAGCGTTCGCGAAGGCCCCTGCCACTGAGCTCGGCGCCTTg                                                             <  2:153919/147‑1 (MQ=255)
                                                                                         cccAATACGCCGCCAGACACGGCGATTGAGGTTGCCGAGCGCCTGCGGGCCCTGATAGCCAAACTGCGTGTGAATGATGAACACCCGGTCACCGCAAGCATAGGTTTAGCGTTCGCGAAGGCCCCTGCCACTGAGCTCGGCGCCTTGCt                                                           >  1:4683/1‑149 (MQ=255)
                                                                                              tACGCCGCCAGACACGGCGATTGAGGTTGCCGAGCGCCTGCGGGCCCTGATAGCCAAACTGAGTGTGAATGATGAACACCCGGTCACCGCAAGCATAGGTTTAGCGTTCGCGAAGGCCCCTGCCACTGAGCTCGGCGCCTTGCTcagca                                                      <  2:98086/149‑1 (MQ=255)
                                                                                                   cgccAGACACGGCGATTGAGGTTGCCGTGCGCCTGCGGGCCCTGATAGCCAAACTGCGTGTGATTGATGAACACCCGGTCACCGCAAGCATAGGTTTAGCGTTCGCGAAGGCCCCTGCCACTGAGCTCGGCGCCTTGCTCAGCAGTGCa                                                 <  7:94785/149‑1 (MQ=255)
                                                                                                   cgccAGACACGGCGATTGAGGTTGCCGAGCGCCTGCGGGCCCTGATAGCCAAACTGCGTGTGAATGATGAACACCCGGTCACCGCAAGCATAGGTTTAGCGTTCGCGAAGGCCCCTGCCACTGAGCTCGGCGCCTTGCTCAGCAGTGCa                                                 <  6:255859/149‑1 (MQ=255)
                                                                                                    gccAGACACGGCGATTGAGGTTGCCGAGCGCCTGCGGGCCCTGATAGCCAAACTGCGTGTGAATGATGAACACCCGGTCACCGCAAGCATAGGTTTAGCGTTCGCGAAGGCCCCTGCCACTGAGCTCGGCGCCTTGCTcagca                                                      <  1:320925/143‑1 (MQ=255)
                                                                                                    gccAGACACGGCGATTGAGGTTGCCGAGCGCCTGCGGGCCCTGATAGCCAAACTGCGTGTGAATGATGAACACCCGGTCACCGCAAGCATAGGTTTAGCGTTCGCGAAGGCCCCTGCCACTGAGCTCGGCGCCTTGCTcagca                                                      >  2:320925/1‑143 (MQ=255)
                                                                                                       aGACACGGCGATTGAGGTTGCCGAGCGCCTGCGGGCCCTGATAGCCAAACTGCGTGTGAATGATGAACACCCGGTCACCGCAAGCATAGGTTTa                                                                                                    >  1:101333/1‑94 (MQ=255)
                                                                                                       aGACACGGCGATTGAGGTTGCCGAGCGCCTGCGGGCCCTGATAGCCAAACTGCGTGTGAATGATGAACACCCGGTCACCGCAAGCATAGGTTTa                                                                                                    <  2:101333/94‑1 (MQ=255)
                                                                                                         acacGGCGATTGAGGTTGCCGAGCGCCTGCGGGCCCTGATAGCCAAACTGCGTGTGAATGATGAACACCCGGTCACCGCAAGCATAGGTTTAGCGTTCGCGAAGGCCCCTGCCACTGAGCTCGGCGCCTTGCTCAGCAGTGCAGACAgc                                           >  5:135774/1‑149 (MQ=255)
                                                                                                               cGATTGAGGTTGCCGAGCGCCTGCGGGCCCTGATAGCCAAACTGCGTGTGAATGATGAACACCCGGTCACCGCAAGCATACGTTTACCGTTCGCGAATGCCCCTGCCACTGAGCTCGGCGCCTTGCTCAGCAGTGCAGACAGCGCCATg                                     <  4:53565/149‑1 (MQ=255)
                                                                                                                  ttGAGGTTGCCGAGCGCCTGCGGGCCCTGATAGCCAAACTGCGTGTGAATGATGAACACCCGGTCACCGCAAGCATAGGTTTAGCGTTCGCGAAGGCCCCTGCCACTGAGCTCGGCGCCTTGCTCAGCAGTGCAGACAGCGCCATGTAc                                  <  2:40194/149‑1 (MQ=255)
                                                                                                                     aGGTTGCCGAGCGCCTGCGGGCCCTGATAGCCAAACTGCGTGTGAATGATGAACACCCGGTCACCGCAAGCATAGGTTTAGCGTTCGCGAAGGCCCCTGCCACTGAGCTCGGCGCCTTGCTCAGCAGTGCAGACAGCGCCATGTACCtt                               >  4:83124/1‑149 (MQ=255)
                                                                                                                       gTTGCCGAGCGCCTGCGGGCCCTGATAGCCAAACTGCGTGTGAATGATGAACACCCGGTCACCGCAAGCATAGGTTTAGCGTTCGCGAAGGCCCCTGCCACTGAGCTCGgcgc                                                                 >  7:85802/1‑113 (MQ=255)
                                                                                                                       gTTGCCGAGCGCCTGCGGGCCCTGATAGCCAAACTGCGTGTGAATGATGAACACCCGGTCACCGCAAGCATAGGTTTAGCGTTCGCGAAGGCCCCTGCCACTGAGCTCGgcgc                                                                 <  8:85802/113‑1 (MQ=255)
                                                                                                                        ttGCCGAGCGCCTGCGGGCCCTGATAGCCAAACTGCGTGTGAATGATGAACACCCGGTCACCGCAAGCATAGGTTTAGCGTTCGCGAAGGCCCCTGCCACTGAGCTCGGCGCCTTGCTCAGCAGTGCAGACCGCGCCATGTACCTTTc                             >  5:287530/1‑148 (MQ=255)
                                                                                                                           ccGAGCGCCTGCGGGCCCTGATAGCCAAACTGCGTGTGAATGATGAACACCCGGTCACCGCAAGCATAGGTTTAGCGTTCGCGAAGGCCCCTGCCACTGAGCTCGGCGCCTTGCTCAGCAGTGCAGACAGCGCCATGTACCTTTCCAAg                         <  5:286808/149‑1 (MQ=255)
                                                                                                                                   cTGCGGGCCCTGATAGCCAAACTGCGTGTGAATGATGAACACCCGGTCACCGCAAGCATAGGTTTAGCGTTCGCGAAGGCCCCTGCCACTGAGCTCGGCGCCTTGCTCAGCAGTGCAGACAGCGCCATGTAc                                  <  2:234662/132‑1 (MQ=255)
                                                                                                                                   cTGCGGGCCCTGATAGCCAAACTGCGTGTGAATGATGAACACCCGGTCACCGCAAGCATAGGTTTAGCGTTCGCGAAGGCCCCTGCCACTGAGCTCGGCGCCTTGCTCAGCAGTGCAGACAGCGCCATGTAc                                  >  1:234662/1‑132 (MQ=255)
                                                                                                                                             tGATAGCCAAACTGCGTGTGAATGATGAACACCCGGTCACCGCAAGCATAGGTTTAGCGTTAGCGAAGGCCCCTGCCACTGACCTCGGCGCCTTGCTCAGCAGTGCAGACAGCGCCATGTACCTTTCCAAGGTATCGGGCCGTAACCgg       >  1:17678/1‑149 (MQ=255)
                                                                                                                                               aTAGCCAAACTGCGTGTGAATGATGAACACCCGGTCACCGCAAGCATAGGTTTAGCGTTCGCGAAGGCCCCTGCCACTGAGCTCGGCGCCTTGCTCAGCAGTGCAGACAGCGCCATGTACCTTTCCAAGGTATCGGGCCGTAACCgg       <  8:94084/147‑1 (MQ=255)
                                                                                                                                                  gCCAAACTGCGTGTGAATGATGAACACCCGGTCACCGCAAGCATAGGTTTAGCGTTCGCGAAGGCCCCTGCCACTGAGCTCGGCGCCTTGCTCAGCAGTGCAGACAGCGCCATGTACCTTTCCAAGGTATCGGGCCGTAACCGGGTCAc  <  1:327977/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TGCAGCGCATCGGTGCAGCGATCAGTGCCAGCCTGCGTGACAACGACATTGCCGGCAGGCTGGGCGGCGAGGAGTTCGCCATCCTGCTGCCCAATACGCCGCCAGACACGGCGATTGAGGTTGCCGAGCGCCTGCGGGCCCTGATAGCCAAACTGCGTGTGAATGATGAACACCCGGTCACCGCAAGCATAGGTTTAGCGTTCGCGAAGGCCCCTGCCACTGAGCTCGGCGCCTTGCTCAGCAGTGCAGACAGCGCCATGTACCTTTCCAAGGTATCGGGCCGTAACCGGGTCAC  >  NC_002947/3785614‑3785908

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