Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,416,652 G→T 6.9% intergenic (+13/+193) PP_3903 → / ← PP_3904 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,416,6520GT 6.9% 184.7 / 3.4 58intergenic (+13/+193)PP_3903/PP_3904hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (28/26);  new base T (2/2);  total (30/28)
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.82e-01

GGCTTCGATGTAAAAGGCGCGATGTGCCGGGACGTAGCCGTTGTAGCTGCGGCCTGGGCGATCGGTGAACTGCAGCGCGAGATGCTGAAGACCATCCAGAAGCCCGGCGGCGGAAAAATCAGCGTCGATTAATCAACTGCCCTAGAGGGCACGCAATGTTCAACATGGCAACCATGGCGGCAGATGAATGCCGCGAATAGGAGAAAGGCTATGAGGCGCTTAATTTGTGCTGGTGTCTATTGGCTCATTGAGCATCTTCTTGACCAACTAGACCGCAGAAAAGCAGCTCGC  >  NC_002947/4416508‑4416798
                                                                                                                                                |                                                                                                                                                  
ggCTTCGATGTAAAAGGCGCGATGTGCCGGGACGTAGCCGTTGTAGCTGCGGCCTGGGCGATCGGTGAACTGCAGCGCGAGATGCTGAAGACCATCCAGAAGCCCGGCGGCGGAAAAATCAGCGTCGATTAATCAACTGCCCTAGAggg                                                                                                                                                >  2:275674/1‑149 (MQ=255)
ggCTTCGATGTAAAAGGCGCGATGTGCCGGGACGTAGCCGTTGTAGCTGCGGCCTGGGCGATCGGTGAACTGCAGCGCGAGATGCTGAAGACCATCCAGAAGCCCGGCGGCGGAAAAATCAGCGTCGATTAATCAACTGCCCTAGAggg                                                                                                                                                <  2:68063/149‑1 (MQ=255)
     cGATGTAAAAGGCGCGATGTGCCGGGACGTACCCGTTGTAGCTGCGGCCTGGGCGATCGGTGAACTGCAGCGCGAGATGCTGAAGACCATCCAGAAGCCCGGCGGCGGAAAAATCAGCGTCGATTAATCAACTGCCCTCTAGGGCagt                                                                                                                                            <  7:280675/148‑3 (MQ=255)
       aTGTAAAAGGCGCGATGTGCCGGGACGTAGCCGTTGTAGCTGCGGCCTGGGCGATCGGTGAACTGCAGCGCGAGATGCTGAAGACCATCCAGAAGCCCGGCGGCGGAAAAATCAGCGTCGATTAATCAACTGCCCTAGAGGGCACGCaa                                                                                                                                         >  7:253427/1‑149 (MQ=255)
        tGTAAAAGGCGCGATGTGCCGGGACGTAGCCGTTGTAGCTGCGGCCTGGGCGATCGGTGAACTGCAGCGCGAGATGCTGAAGACCATCCAGAAGCCCGGCGGCGGAAAAATCAGCGTCGATTAATCAACTGCCCTAGAGGGCACGCAAt                                                                                                                                        <  4:202824/149‑1 (MQ=255)
        agTAAAAGGCGCGATGTGCCGGGACGTAGCCGATGTAGCTCCGGCCTGGCAGATCCCTGAACAGCAGCGACAGATCCTGAAGACCAACCAGAAGCCCGGCGGCGGAAAAATCAGCGTCGATTAATCAACTGCACTAGACGACACCCAAt                                                                                                                                        >  1:276742/2‑149 (MQ=255)
          tAAAAGGCGCGATGTGCCGGGACGTAGCCGTTGTAGCTGCGGCCTGGGCGATCGGTGAACTGCAGCGCGAGATGCTGAAGACCATCCAGAAGCCCGGCGGCGGAAAAATCAGCGTCGATTAATCAACTGCCCTAGAGGGCACGCAATg                                                                                                                                       <  7:34685/148‑1 (MQ=255)
          tAAAAGGCGCGATGTGCCGGGACGTAGCCGTTGTAGCTGCGGCCTGGGCGATCGGTGAACTGCAGCGCGAGATGCTGAAGACCATCCAGAAGCCCGGCGGCGGAAAAATCAGCGTCGATTAATCAACTGCCCTAGAGGGCACGCAATg                                                                                                                                       <  1:275674/148‑1 (MQ=255)
          tAAAAGGCGCGATGTGCCGGGACGTAGCCGTTGTAGCTGCGGCCTGGGCGATCGGTGAACTGCAGCGCGAGATGCTGAAGACCATCCAGAAGCCCGGCGGCGGAAAAATCAGCGTCGATTAATCAACTGCCCTAGAGGGCACGCAATGt                                                                                                                                      >  6:183164/1‑149 (MQ=255)
            atagtggcgATGTGCCGGGACGTAGCCGTTGTAGCTGCGGCTTGGGCGATGGGTTAACTGCAGCGCGAGATGCTGAAGACCATCCAGAAGCCCGGCGTCGGAAAATTCAGCGTCGATTAATCAACTGCCCTAGAGGGCACGCAATGTTc                                                                                                                                    <  2:276742/143‑1 (MQ=255)
              aGGCGCGATGTGCCGGGACGTAGCCGTTGTAGCTGCGGCCTGGGCGATCGGTGAACTGCAGCGCGAGATGCTGAAGACCATCCAGAAGCCCGGCGGCGGAAAAATCAGCGTCGATTAATCAACTGCCCTAGAGGGCACGCAATGTTCaa                                                                                                                                  <  4:157176/149‑1 (MQ=255)
              aGGCGCGATGTGCCGGGACGTAGCCGTTGTAGCGGCGGCCTGGGCGATCGGTGAACTGCTGCGCGAGATGCTGAAGACCATCCAGAAGCCCGGCGTCGGAAAAATCAGCGTCGATTAATCAACTGCCCTAGAGTTCACGCAATTTTCaa                                                                                                                                  <  4:249295/149‑1 (MQ=255)
                   cgATGTGCCGGGACGTAGCCGTTGTAGCTGCGGCCTGGGCGATCGGTGAACTGCAGCGCGAGATGCTGAAGACCATCCAGAAGCCCGGCGGCGGAAAAATCAGCGTCGATTAATCAACTGCCCTAGAGGGCACGCAATGTTCAACATgg                                                                                                                             <  3:283046/149‑1 (MQ=255)
                    gATGTGCCGGGACGTAGCCGTTGTAGCTGCGGCCTGGGCGATCGGTGAACTGCAGCGCGAGATGCTGAAGACCATCCAGAAGCCCGGCGGCGGAAAAATCAGCGTCGATTAATCAACTGCCCTAGAGGGCACGCAATGTTCAACATgg                                                                                                                             >  5:157118/1‑148 (MQ=255)
                              gACGTAGCCGTTGTAGCTGCGGCCTGGGCGATCGGTGAACTGCAGCGCGAGATGCTGAAGACCATCCAGAAGCCCGGCGGCGGAAAAATCAGCGTCGATTAATCAACTGCCCTAGAGGGCACGCAATGTTCAACATGGCAACCATggcg                                                                                                                  <  7:251508/149‑1 (MQ=255)
                              gACGTAGCCGTTGTAGCTGCGGCCTGGGCGATCGGTGAACTGCAGCGCGAGATGCTGAAGACCATCCAGAAGCCCGGCGGCGGAAAAATCAGCGTCGATTAATCAACTGCCCTAGAGGGCACGCAATGTTCAACATGGCAACCATggcg                                                                                                                  <  7:271338/149‑1 (MQ=255)
                               aCGTAGCCGTTGTAGCTGCGGCCTGGGCGATCGGTGAACTGCAGCGCGAGATGCTGAAGACCATCCAGAAGCCCGGCGGCGGAAAAATCAGCGTCGATTAATCAACTGCCCTAGAGGGCACGCAATGTTCAACATGGCAACCATggcgg                                                                                                                 <  8:253427/149‑1 (MQ=255)
                                  tAGCCGTTGTAGCTGCGGCCTGGGCGATCGGTGAACTGCAGCGCGAGATGCTGAAGACCATCCAGAAGCCCGGCGGCGGAAAAATCAGCGTCGATTAATCAACTGCCCTAGAGGGCACGCAATGTTCAACATGGCAACCATggcgg                                                                                                                 >  5:107076/1‑146 (MQ=255)
                                              cTGAGGCCTGGGCGATCGGTGAACTGCAGCGCGAGATGCTGAAGACCATCCAGAAGCCCGGCGGCGGAAAATCAGCGTCGATTAATCAACTGCCCTCTAGGGCagt                                                                                                                                            <  5:191485/106‑3 (MQ=255)
                                                 cGGCCTGGGCGATCGGTGAACTGCAGCGCGAGATGCTGAAGACCATCCAGAAGCCCGGCGGCGGAAAAATCAGCGTCGATTAATCAACTGCCCTAGAGGGCACGCAATGTTCAACATGGCAACCATGGCGGCAGATGAATGCCGCGAAt                                                                                               <  4:301740/149‑1 (MQ=255)
                                                           gATCGGTGAACTGCAGCGCGAGATGCTGAAGACCATCCAGAAGCCCGGCGGCGGAAAAATCAGCGTCGATTAATCAACTGCCCTAGAGGGCACGCAATGTTCAACATGGCAACCATGGCGGCAGATGAATGCCGCGAATAGGAGAAAgg                                                                                     <  2:5877/149‑1 (MQ=255)
                                                               ggTGAACTGCAGCGCGAGATGCTGAAGACCATCCAGAAGCCCGGCGGCGGAAAAATCAGCGTCGATTAATCAACTGCCCTAGAGGGCACGCAATGTTCAACATGGCAACCATGGCGGCAGATGAATGCCGCGAATAGGAGAAAGGCTAt                                                                                 <  6:157118/149‑1 (MQ=255)
                                                                gTGAACTGCAGCGCGAGATGCTGAAGACCATCCAGAAGCCCGGCGGCGGAAAAATCAGCGTCGATTAATCAACTGCCCTAGAGGGCACGCAATGTTCAACATGGCAACCATGGCGGCAGATGAATGCCGCGAATAGGAGAAAGGCTAt                                                                                 >  4:146976/1‑148 (MQ=255)
                                                                  gAACTGCAGCGCGAGATGCTGAAGACCATCCAGAAGCCCGGCGGCGGAAAAATCAGCGTCGATTAATCAACTGCCCTAGAGGGCACGCAATGTTCAACATGGCAACCATGGCGGCAGATGAATGCCGCGAATAGGAGAAAGGCTATGAg                                                                              >  8:35702/1‑149 (MQ=255)
                                                                   aaCTGCAGCGCGAGATGCTGAAGACCATCCAGAAGCCCGGCGGCGGAAAAATCAGCGTCGATTAATCAACTGCCCTAGAGGGCACGCAATGTTCAACATGGCAACCATGGCGGCAGATGAATGCCGCGAATAGGAGAAAGGCTATGAgg                                                                             >  4:142537/1‑149 (MQ=255)
                                                                     cTGCAGCGCGAGATGCTGAAGACCATCCAGAAGCCCGGCGGCGGAAAAATCAGCGTCGATTAATCAACTGCCCTAGAGGGCACGCAATGTTCAACATGGCAACCATGGCGGCAGATGAATGCCGCGAATAGGAGAAAGGCTATGAGgcg                                                                           >  6:148387/1‑149 (MQ=255)
                                                                                gaTGCTGAAGACCATCCAGAAGCCCGGCGGCGGAAAAATCAGCGTCGATTAATCAACTGCCCTAGAGGGCACGCAATGTTCAACATGGCAACCATGGCGGCAGATGAATGCCGCGAATAGGAGAAAGGCTATGAGGCGCTTAATTtgtg                                                                <  6:107076/149‑1 (MQ=255)
                                                                                    cTGAAGACCATCCAGAAGCCCGGCGGCGGAAAAATCAGCGTCGATTAATCAACTGCCCTAGAGGGCACGCAAAGTTCAACATGGCAACCATGGCGGCAGATGAATGCCGCGAATAGGAGAAAGGCTATGAGGCGCTTAATTTGTGCTgg                                                            >  7:51466/1‑149 (MQ=255)
                                                                                      gAAGACCATCCAGAAGCCCGGCGGCGGAAAAATCAGCGTCGATTAATCAACTGCCCTAGAGGGCACGCAATGTTCAACATGGCAACCATGGCGGCAGATGAATGCCGCGAATAGGAGAAAGGCTATGAGGCGCTTAATTTGTGCTGgtg                                                          >  7:91131/1‑149 (MQ=255)
                                                                                       aaGACCATCCAGAAGCCCGGCGGCGGAAAAATCAGCGTCGATTAATCAACTGCCCTAGAGGGCACGCAATGTTCAACATGGCAACCATGGCGGCAGATGAATGCCGCGAATAGGAGAAAGGCTATGAGGCGCTTAATTTGTGCTGgt                                                           <  6:165550/147‑1 (MQ=255)
                                                                                          aCCATCCAGAAGCCCGGCGGCGGAAAAATCAGCGTCGATTAATCAACTGCCCTAGAGGGCACGCAATGTTCAACATGGCAACCATGGCGGCAGATGAATGCCGCGAATAGGAGAAAGGCTATGAGGCGCTTAATTTGTGCTGGTGTCTa                                                      >  2:50344/1‑149 (MQ=255)
                                                                                            cATCCAGAAGCCCGGCGGCGGAAAAATCAGCGTCGATTAATCAACTGCCCTAGAGGGCACGCAATGTTCAACATGGCAACCATGGCGGCAGATGAATGCCGCGAATAGGAGAAAGGCTTTGAGGCGCTTAATTTGTGCTGGTGTCTAtt                                                    >  8:273449/1‑149 (MQ=255)
                                                                                              tCCAGAAGCCCGGCGGCGGAAAAATCAGCGTCGATTAATCAACTGCCCTAGAGGGCACGCAATGTTCAACATGGCAACCATGGCGGCAGATGAATGCCGCGAATAGGAGAAAGGCTATGAGGCGCTTAATTTGTGCTGGTGTCTATTgg                                                  >  3:129367/1‑149 (MQ=255)
                                                                                                cAGAAGCCCGGCGGCGGAAAAATCAGCGTCGATTAATCAACTGCCCTAGAGGGCACGCAATGTTCAACATGGCAACCATGGCGGCAGATGAATGCCGCGAATAGGAGAAAGGCTATGAGGCGCTTAATTTGTGCTGGTGTCTATTGGCt                                                >  4:264530/1‑149 (MQ=255)
                                                                                                 aGAAGCCCGGCGGCGGAAAAATCAGCGTCGATTAATCAACTGCCCTAGAGGGCACGCAATGTTCAACATGGCAACCATGGCGGCAGATGAATGCCGCGAATAGGAGAAAGGCTATGAGGCGCTTAATTTGTGCTGGTGTCTATTGGCTc                                               <  8:91131/149‑1 (MQ=255)
                                                                                                        cggcggcggAAAAATCAGCGTCGATTAATCAACTGCCCTAGAGGGCACGCAATGTTCAACATGGCAACCATGGCGGCAGAt                                                                                                            >  4:109161/1‑81 (MQ=255)
                                                                                                        cggcggcggAAAAATCAGCGTCGATTAATCAACTGCCCTAGAGGGCACGCAATGTTCAACATGGCAACCATGGCGGCAGAt                                                                                                            <  3:109161/81‑1 (MQ=255)
                                                                                                         ggcggcggAAAAATCAGCGTCGATTAATCAACTGCCCTAGAGGGCACGCAATGTTCAACATGGCAACCATGGCGGCAGATGAATGCCGCGAATAGGAGAAAGGCTATGAGGCGCTTAATTTGTGCTGGTGTCTATTGGCTCATTGAGCa                                       >  6:18791/1‑149 (MQ=255)
                                                                                                               ggAAAAATCAGCGTCGATTAATCAACTGCCCTAGAGGGCACGCAATGTTCAACATGGCAACCATGGCGGCAGATGAATGCCGCGAATAGGAGAAAGGCTATGAGGCGCTTAATTTGTGCTGGTGTCTATTGGCTCATTGAg                                         >  4:257076/1‑141 (MQ=255)
                                                                                                               ggAAAAATCAGCGTCGATTAATCAACTGCCCTAGAGGGCACGCAATGTTCAACATGGCAACCATGGCGGCAGATGAATGCCGCGAATAGGAGAAAGGCTATGAGGCGCTTAATTTGTGCTGGTGTCTATTGGCTCATTGAg                                         <  3:257076/141‑1 (MQ=255)
                                                                                                                gAAAAATCAGCGTCGATTAATCAACTGCCCTAGAGGGCACGCAATGTTCAACATGGCAACCATGGCGGCAGATGAATGCCGCGAATAGGAGAAAGGCTATGAGGCGCTTAATTTGTGCTGGTGTCTATTGGCTCATTGAGCATcttctt                                >  4:202238/1‑149 (MQ=255)
                                                                                                                 aaaaaTCAGCGGCGATTAATCAACTGCCCTAGAGGGCACGCAATGTTCAACATGGCAACCATGGCGGCAGATGAATGCCGCGAATAGGAGAAAGGCTATGAGGCGCTTAATTTGTGCTGGTGTCTATTGGCTCATTGAGCATCTTCTTg                               <  6:274884/149‑1 (MQ=255)
                                                                                                                  aaaaTCAGCGTCGATTAATCAACTGCCCTAGAGGGCACGCAATGTTCAACATGGCAACCATGGCGGCAGATGAATGCCGCGAATAGGAGAAAGGCTATGAGGCGCTTAATTTGTGCTGGTGTCTATTGGCTCATTGAGCATCTTCTTGa                              >  5:205944/1‑149 (MQ=255)
                                                                                                                     aTCAGCGTCGATTAATCAACTGCCCTAGAGGGCACGCAATGTTCAACATGGCAACCATGGCGGCAGATGAATGCCGCGAATAGGAGAAAGGCTATGAGGCGCTTAATTTGTGCTGGTGTCTATTGGCTCATTGAGCATCTTCTTGACCa                           <  7:35702/149‑1 (MQ=255)
                                                                                                                        aGCGTCGATTAATCAACTGCCCTAGAGGGCACGCAATGTTCAACATGGCAACCATGGCGGCAGATGAATGCCGCGAATAGGAGAAAGGCTATGAGGCGCTTAATTTGTGCTGGTGTCTATTGGCTCATTGAGCATCTTCTTGACCAAc                         <  3:142537/148‑1 (MQ=255)
                                                                                                                         gCGTCGATTAATCAACTGCCCTAGAGGGCACGCAATGTTCAACATGGCAACCATGGCGGCAGATGAATGCCGCGAATAGGAGAAAGGCTATGAGGCGCTTAATTTGTGCTGGTGTCTATTGGCTCATTGAGCATCTTCTTGACCAACTa                       <  6:30529/149‑1 (MQ=255)
                                                                                                                           gTCGATTAATCAACTGCCCTAGAGGGCACGCAATGTTCAACATGGCAACCATGGCGGCAGATGAATGCCGCGAATAGGAGAAAGGCTATGAGGCGCTTAATTTGTGCTGGTGTCTATTGGCTCAtt                                            >  8:232658/1‑126 (MQ=255)
                                                                                                                           gTCGATTAATCAACTGCCCTAGAGGGCACGCAATGTTCAACATGGCAACCATGGCGGCAGATGAATGCCGCGAATAGGAGAAAGGCTATGAGGCGCTTAATTTGTGCTGGTGTCTATTGGCTCAtt                                            <  7:232658/126‑1 (MQ=255)
                                                                                                                            tCGATTAATCAACTGCCCTAGAGGGCACGCAATGTTCAACATGGCAACCATGGCGGCAGATGAATGCCGCGAATATGAGAAAGGCTATGAGGCGCTTAATTTGTGCTGGTGTCCATTGGCTCATTGAGCATCTTCTTGACCAACTAGAc                    >  3:216340/1‑149 (MQ=255)
                                                                                                                              gATTAATCAACTGCCCTAGAGGGCACGCAATGTTCAACATGGCAACCCTGGCGGCAGCTGAATGCCGCGAATAGGAGAAAGGCTATGAGGCGCTTAATTTGTGCTGGGGTCTATTGGCTAATTGAGCATCTTCTTGACCAACTATAcac                  >  4:36158/1‑147 (MQ=255)
                                                                                                                               aTTAATCAACTGCCCTAGAGGGCACGCAATGTTCAACATGGCAACCATGGCGGCAGATGAATGCCGCGAATAGGAGAAAGGCTATGAGGCGCTTAATTTGTGCTGGTGTCTATTGGCTCATTGAGCATCTTCTTGACCAACTAGACCg                  >  4:126851/1‑148 (MQ=255)
                                                                                                                                ttAATCAACTGCCCTAGAGGGCACGCAATGTTCAACATGGCAACCATGGCGGCAGATGAATGCCGCGAATAGGAGAAAGGCTATGAGGCGCTTAATTTGTGCTGGTGTCTATTGGCTCATTGAGCATCTTCTTGACCAACTAGACCGCa                <  4:129367/149‑1 (MQ=255)
                                                                                                                                  aaTCAACTGCCCTAGAGGGCACGCAATGTTCAACATGGCAACCATGGCGGCAGATGAATGCCGCGAATAGGAGAAAGGCTATGAGGCGCTTAATTTGTGCTGGTGTCTATTGGCTCATTGAGCATCTTCTTGa                              >  7:80310/1‑133 (MQ=255)
                                                                                                                                  aaTCAACTGCCCTAGAGGGCACGCAATGTTCAACATGGCAACCATGGCGGCAGATGAATGCCGCGAATAGGAGAAAGGCTATGAGGCGCTTAATTTGTGCTGGTGTCTATTGGCTCATTGAGCATCTTCTTGa                              <  8:80310/133‑1 (MQ=255)
                                                                                                                                   attgcgtGCCCTCTAGGGCACGCAATGTTCAACATGGCAACCATGGCGGCAGATGAATGCCGCGAATAGGAGAAAGGCTATGAGGCGCTTAATTTGTGCTGGTGTCTATTGGCTCATTGAGCATCTTCTTGACCAACTAGACCGCAGa              >  3:294007/7‑148 (MQ=255)
                                                                                                                                     tgcgtGCCCTCTAGGGCACGCAATGTTCAACATGGCAACCATGGCGGCAGATGAATGCCGCGAATAGGAGAAAGGCTATGAGGCGCTTAATTTGTGCTGGTGTCTATTGGCTCATTGAGCATCTTCTTGACCAACTAGACCGCAGaa             >  7:142230/5‑147 (MQ=255)
                                                                                                                                            ccTAGAGGGCACGCAATGTTCAACATGGCAACCATGGCGGCAGATGAATGCCGCGAATAGGAGAAAGGCTATGAGGCGCTTAATTTGTGCTGGTGTCTATTGGCTCATTGAGCATCTTCTTGACCAACTAGACCGCAGAAAAGCAGCTc    >  7:229142/1‑149 (MQ=255)
                                                                                                                                              tAGAGGGCACGCAATGTTCAACATGGCAACCATGGCGGCAGATGAATGCCGCGAATAGGAGAAAGGCTATGAGGCGCTTAATTTGTGCTGGTGTCTATTGGCTCATTGAGCATCTTCTTGACCAACTAGACCGCAGAAAAGCAGCTcgc  >  3:71996/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
GGCTTCGATGTAAAAGGCGCGATGTGCCGGGACGTAGCCGTTGTAGCTGCGGCCTGGGCGATCGGTGAACTGCAGCGCGAGATGCTGAAGACCATCCAGAAGCCCGGCGGCGGAAAAATCAGCGTCGATTAATCAACTGCCCTAGAGGGCACGCAATGTTCAACATGGCAACCATGGCGGCAGATGAATGCCGCGAATAGGAGAAAGGCTATGAGGCGCTTAATTTGTGCTGGTGTCTATTGGCTCATTGAGCATCTTCTTGACCAACTAGACCGCAGAAAAGCAGCTCGC  >  NC_002947/4416508‑4416798

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: