Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,079,650 A→G 9.8% intergenic (‑385/+49) PP_2688 ← / ← PP_2689 hypothetical protein/endoribonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,079,6500AG 9.8% 138.0 / 5.4 51intergenic (‑385/+49)PP_2688/PP_2689hypothetical protein/endoribonuclease
Reads supporting (aligned to +/- strand):  ref base A (22/24);  new base G (3/2);  total (25/26)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.41e-01

CGGCGCATCAAAAAAGGGCGCTTTTGACCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGA‑CGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTCCACCAGCATTTCTTCCCGGGCGAGGCCAGCCACCACTACGGCCGCGCGGTTGGGGTAGGGTGCCTG  >  NC_002947/3079502‑3079791
                                                                                                                                                    |                                                                                                                                              
cGGCGCATCAAAAAAGGGCGCTTTTGACCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCa                                                                                                                                                >  8:4315/1‑149 (MQ=255)
cGGCGCATCAAAAAAGGGCGCTTTTGACCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCa                                                                                                                                                >  5:269129/1‑149 (MQ=255)
     cATCAAAAAAGGGCGCTTTTGACCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGa                                                                                                                                           <  7:194621/149‑1 (MQ=255)
      aTCAAAAAAGGGCGCTTTTGACCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACt                                                                                                                                              <  4:131342/145‑1 (MQ=255)
      aTCAAAAAAGGGCGCTTTTGACCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACt                                                                                                                                              >  3:131342/1‑145 (MQ=255)
      aTCAAAAAAGGGCGCTTTTGACCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGAt                                                                                                                                          >  1:38364/1‑149 (MQ=255)
      aTCAAAAAAGGGCGCTTTTGACCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGAt                                                                                                                                          >  2:86898/1‑149 (MQ=255)
        cAAAAAAGGGCGCTTTTGACCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTg                                                                                                                                        >  8:42667/1‑149 (MQ=255)
        cAAAAAAGGGCGCTTTTGACCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTg                                                                                                                                        >  1:130115/1‑149 (MQ=255)
          aaaaaGGGCGCTTTTGACCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTgcg                                                                                                                                      >  5:27161/1‑149 (MQ=255)
           aaaaGGGCGCTTTTGACCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTgcgc                                                                                                                                     <  6:269129/149‑1 (MQ=255)
                          aCCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGATTATCg                                                                                                                      >  1:280142/1‑149 (MQ=255)
                          aCCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCg                                                                                                                      <  7:290266/149‑1 (MQ=255)
                           cccGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGt                                                                                                                     <  5:124831/149‑1 (MQ=255)
                           cccGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGt                                                                                                                     >  5:59614/1‑149 (MQ=255)
                           cccGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGt                                                                                                                     <  3:56418/149‑1 (MQ=255)
                            ccGGGCGCGGCTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTg                                                                                                                    <  8:152355/149‑1 (MQ=255)
                             cGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTga                                                                                                                   <  7:42667/149‑1 (MQ=255)
                                               cGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGA‑CGGGCGCCTTCAGGGc                                                                                                <  6:27161/149‑1 (MQ=255)
                                                      gggCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGA‑CGGGCGCCTTCAGGGCGTCaga                                                                                          <  1:131749/148‑1 (MQ=255)
                                                          gcAGGCCGGTGCCTGGAGGGAGAAAAAGCGAGTATTCCGGGTGCATGACGTGCGGAGGAGGGGGCATGCATCTAGCCCCGCGAGGGCGGCACACGATTGCGCGAACGCTCAATATCTTGAGA‑CGGGCGCCTTCAGCGCGTCAGaggcag                                                                                     >  7:294660/1‑149 (MQ=255)
                                                            aGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGA‑CGGGCGCCTTCAGGGCGTCAGaggcaggc                                                                                   >  4:103713/1‑149 (MQ=255)
                                                            aGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGA‑CGGGCGCCTTCAGGGCGTCAGaggcaggc                                                                                   <  2:192161/149‑1 (MQ=255)
                                                            aGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGA‑CGGGCGCCTTCAGGGCGTCAGaggcaggc                                                                                   <  4:247864/149‑1 (MQ=255)
                                                                     ccTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGA‑CGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGaccac                                                                          <  1:86898/149‑1 (MQ=255)
                                                                             gagaAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGA‑CGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTcca                                                                  >  2:21745/1‑149 (MQ=255)
                                                                                  aaaGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGA‑CGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTCCACCAGc                                                             <  6:59614/149‑1 (MQ=255)
                                                                                      ggAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAat                                                                                                                          <  6:109671/85‑1 (MQ=255)
                                                                                      ggAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAat                                                                                                                          >  5:109671/1‑85 (MQ=255)
                                                                                       gAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGA‑CGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTCCACCAGCATttc                                                        <  1:238643/149‑1 (MQ=255)
                                                                                           aTTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGA‑CGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTCCACCAGCATTTCTTcc                                                    >  6:6548/1‑149 (MQ=255)
                                                                                                  gTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGt                                                                                                                     >  4:212271/1‑78 (MQ=255)
                                                                                                        gACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGA‑CGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTCCACCAGCATTTCTTCCCGGGCGAGgccag                                       >  8:199847/1‑149 (MQ=255)
                                                                                                           gTGGGGCGGAGGGGGCATGCATCGAGTGCCGCCCGCGCGGCGCt                                                                                                                                              >  6:158120/1‑44 (MQ=17)
                                                                                                           gTGGGGCGGAGGGGGCATGCATCGAGTGCCGCCCGCGCGGCGCt                                                                                                                                              <  5:158120/44‑1 (MQ=17)
                                                                                                              gggCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGA‑CGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTCCACCAGCATTTCTTCCCGGGCGAGGCCAGccacc                                  <  2:130115/148‑1 (MQ=255)
                                                                                                                  ggAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGA‑CGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTCCACCAGCATTTCTTCCCGGGCGAGGCCAGCCACCACTAc                             >  8:31844/1‑149 (MQ=255)
                                                                                                                      ggggCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGA‑CGGGCGCCTTCAGGGCGGCAGAGGCAGTCGTAGACCACCAGTTCCACCAGCATTTCTTCCCGGGCGACGCCAGCCACCACTACGGcc                         <  1:78216/149‑1 (MQ=255)
                                                                                                                         gcatgcatCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGA‑CGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTCCACCAGCATTTCTTCCCGGGCGAGGCCAGCCACCACTACGGCcgcg                      <  7:31844/149‑1 (MQ=255)
                                                                                                                           cgcaatCGAGTGCCGCCCGCGCGGCGCTCGATTGCGCGAACGCTGAATATCGTGAGA‑CGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTCCACCAGCATTTCTTCCCGGGCGAGGCCAGCCACCACTACGGCcgcgcg                    >  2:115122/5‑149 (MQ=255)
                                                                                                                           atgcatCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGA‑CGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGCCCACCAGTTCCACCAGCATTTCTTCCCGGGCGAGGCCAGCCACCCCTACGGCcgcgcg                    >  7:18498/1‑149 (MQ=255)
                                                                                                                           atgcatCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGA‑CGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTCCACCAGCATTTCTTCCCGGGCGAGGCCAGCCACCACTACGGCcgcgcg                    <  1:21745/149‑1 (MQ=255)
                                                                                                                              aatCGAGTGCCGCCCGCGCGGCGCTCGATTGCGCGAACGCTGAATATCGTGAGA‑CGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTCCACCAGCATTTCTTCCCGGGCGAg                                            <  8:229075/121‑1 (MQ=255)
                                                                                                                              aatCGAGTGCCGCCCGCGCGGCGCTCGATTGCGCGAACGCTGAATATCGTGAGA‑CGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTCCACCAGCATTTCTTCCCGGGCGAg                                            >  7:229075/2‑122 (MQ=255)
                                                                                                                                 cGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGAGCGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTCCACCAGCATTTCTTCCCGGGCGAGGCCAGCCACCACTACGGCCGCGCGGTTgg               >  1:241725/1‑149 (MQ=255)
                                                                                                                                  gAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGA‑CGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTCCACCAGCATTTCTTCCCGGGCGAGGCCAGCCACCACTACGGCCGCGCGGTTggg              <  2:38364/148‑1 (MQ=255)
                                                                                                                                   aGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGA‑CGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTCCACCAGCATTTCTTCCCGGGCGAGGCCAGCCACCACTACGGCCGCGCGGTTgggg             >  1:63530/1‑148 (MQ=255)
                                                                                                                                         gcgcgGGCGGCACTCGATTGCGCGATCGCTGAATATCGTGAGA‑CGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTCCACCAGCATTTCTTCCCGGGCGAGGCCAGCCACCACTACGGCCGCGCGGTTGGGGTAGGGTg      >  3:148736/1‑149 (MQ=255)
                                                                                                                                             gggcggCACTCGATTGCGCGAACGCTGAATATCGTGAGA‑CGGGCGCCTTCAGGGCGGCAGAGGCAGGCGTAGACCACCAGTTCCACCAGCATTTCTTCCCGGGCGAGGCCAGCCACCACTACGGCCGCGCGGTTGGGGTAGGGTGCCTg  <  5:6548/149‑1 (MQ=255)
                                                                                                                                                 ggcACTCGATTGCGCGAACGCTGAATATCGTGAGA‑CGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTCCACCAGCATTTCTTCCCGGGCGAGGCCAGCCACCACTACGGCCGCGCGGTTGGGGTAGGGTGCCTg  >  4:196353/1‑145 (MQ=255)
                                                                                                                                                 ggcACTCGATTGCGCGAACGCTGAATATCGTGAGA‑CGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTCCACCAGCATTTCTTCCCGGGCGAGGCCAGCCACCACTACGGCCGCGCGGTTGGGGTAGGGTGCCTg  <  3:196353/145‑1 (MQ=255)
                                                                                                                                                   cACTCGATTGCGCGAACGCTGAATATCGTGAGA‑CGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTCCACCAGCATTTCTTCCCGGGCGAGGCCAGCCACCACTACGGCCGCGCGGTTGGGGTAGGGTGc     >  6:97494/1‑140 (MQ=255)
                                                                                                                                                   cACTCGATTGCGCGAACGCTGAATATCGTGAGA‑CGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTCCACCAGCATTTCTTCCCGGGCGAGGCCAGCCACCACTACGGCCGCGCGGTTGGGGTAGGGTGc     <  5:97494/140‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
CGGCGCATCAAAAAAGGGCGCTTTTGACCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGA‑CGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTCCACCAGCATTTCTTCCCGGGCGAGGCCAGCCACCACTACGGCCGCGCGGTTGGGGTAGGGTGCCTG  >  NC_002947/3079502‑3079791

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