Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,085,456 G→C 9.2% intergenic (+165/‑430) PP_1864 → / → PP_1865 hypothetical protein/transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,085,4560GC 9.2% 104.1 / 2.2 43intergenic (+165/‑430)PP_1864/PP_1865hypothetical protein/transposase
Reads supporting (aligned to +/- strand):  ref base G (21/16);  new base C (2/2);  total (24/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGT  >  NC_002947/2085318‑2085601
                                                                                                                                          |                                                                                                                                                  
ccAGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGTATGCCCGTc                                                                                                                                           >  5:308337/1‑148 (MQ=255)
        tgcttctgcttctgcttctgcttctgctCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTgg                                                                                                                                 <  3:10966/149‑1 (MQ=255)
          cttctgcttctgcttctgcttctgctCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTgg                                                                                                                                 <  4:149173/147‑1 (MQ=255)
           ttctgcttctgcttctgcttctgctCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATCCATGCCCGTCa                                                                                                                                          >  1:47416/1‑138 (MQ=255)
           ttctgcttctgcttctgcttctgctCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATCCATGCCCGTCa                                                                                                                                          <  2:47416/138‑1 (MQ=255)
           ttctgcttctgcttctgcttctgctCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTAAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGTAGGGAGGTGACGGGCATGGATGACAGTCAAGCGCttgtcc                                                                                                                               >  5:168668/1‑144 (MQ=17)
              tgcttctgcttctgcttctgctCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGcc                                                                                                                              >  6:256093/1‑146 (MQ=255)
              tgcttctgcttctgcttctgctCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGcc                                                                                                                              <  5:256093/146‑1 (MQ=255)
                cttctgcttctgcttctgctCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAg                                                                                                                          >  8:241465/1‑148 (MQ=255)
                     gcttctgcttctgctCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATggg                                                                                                                    >  7:80695/1‑149 (MQ=40)
                         ctgcttctgctCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGt                                                                                                                                            <  2:224735/122‑1 (MQ=39)
                         ctgcttctgctCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGt                                                                                                                                            >  1:224735/1‑122 (MQ=39)
                            cttctgctCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCa                                                                                                             >  6:82678/1‑149 (MQ=35)
                             ttctgctCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAg                                                                                                            >  3:204134/1‑149 (MQ=34)
                               ctgctCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAgcg                                                                                                          >  6:34247/1‑149 (MQ=32)
                                  ctctTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAg                                                                                                                          <  5:130012/130‑1 (MQ=32)
                                  ctctTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAg                                                                                                                          >  6:130012/1‑130 (MQ=32)
                                       gCTTTGGCTTTTGATTATAAGCGCGCGGTAGTTCAGGCGATGCATATTGAGACTTCTGGAGGCCGAGCGCAGGGCTTGCTGAGGGAGGTGTCGGGCATGGATGCCCGTCATGCGCTGGGCCCCAGGATGGGGCCTGCTGCGCGGtcct                                                                                                   <  2:139498/148‑1 (MQ=11)
                                        cTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGtcctcc                                                                                                 <  6:250764/149‑1 (MQ=32)
                                         tttGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTccc                                                                                                >  7:206528/1‑149 (MQ=32)
                                             tcTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATTCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCTGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCATGATGGGGCCTGCAGCGCGGTCCTCCCggg                                                                                             <  6:168668/147‑1 (MQ=12)
                                                 ttGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGc                                                                                                                                                >  3:267093/1‑94 (MQ=33)
                                                 ttGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGc                                                                                                                                                <  4:267093/94‑1 (MQ=33)
                                                   gCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAg                                                                                       <  8:206528/148‑1 (MQ=32)
                                                     ttCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGccc                                                                                    >  5:313575/1‑149 (MQ=32)
                                                            cgcgcgGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGt                                                                               >  7:185262/1‑147 (MQ=32)
                                                                                               aggaggCCGAGCGCAGGGCTTGCGGAGGGAGGCGACGGGCATGGAGGCCCGTCAAGCGCGGGGCCCCAGGAGGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAg                                          <  7:124089/149‑1 (MQ=2)
                                                                                                ggccgCCGCGCGCAGGGGGGGCGGCGGGAGGTGCCGGGCAGGCATCCCCGTCAAGCGCGGGGCCCCAGGAGGGGGCCGGCAGCGCGGTCCGCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAgg                                         <  5:34247/145‑1 (MQ=11)
                                                                                                ggaggCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAgg                                         <  6:313575/149‑1 (MQ=31)
                                                                                                      cGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGcc                                   >  7:164122/1‑149 (MQ=32)
                                                                                                           gcAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACgg                              >  8:222407/1‑149 (MQ=32)
                                                                                                           gcAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACgg                              >  7:57697/1‑149 (MQ=32)
                                                                                                               ggCTTGGGGTGGGTGGTGCCGGGCTTGTAGGCCCGTAAAGCGCTGGGCCCCTGTATGGGGCCTGCTGGGCGGTCCTCCCGGGAGAAAGCCCGTAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCttt                          <  8:185262/149‑1 (MQ=11)
                                                                                                               ggCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCttt                          <  5:82678/149‑1 (MQ=34)
                                                                                                                 cTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTgg                        >  1:29681/1‑149 (MQ=34)
                                                                                                                  ttGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGt                       >  1:144657/1‑149 (MQ=34)
                                                                                                                          ggaggTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTgg               <  2:29681/149‑1 (MQ=34)
                                                                                                                          ggaggTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAAGGGAGCCGACGGCTTTGGTTACTTTgg               >  5:303141/1‑149 (MQ=17)
                                                                                                                            aggTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACttt                 >  2:90072/1‑145 (MQ=34)
                                                                                                                            aggTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACttt                 <  1:90072/145‑1 (MQ=34)
                                                                                                                              ttGACGGGCATCCATGCCCGTCAAGCGCTGGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGcaa            >  1:183947/2‑147 (MQ=11)
                                                                                                                                      cATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCATCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAg   >  1:214740/1‑149 (MQ=21)
                                                                                                                                       aTGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGTCGTGGCCAAAGt  <  4:204134/149‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                  
CCAGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGT  >  NC_002947/2085318‑2085601

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: