Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,275,093 G→C 9.3% intergenic (‑284/+28) PP_5462 ← / ← PP_5463 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,275,0930GC 9.3% 128.6 / 4.1 43intergenic (‑284/+28)PP_5462/PP_5463hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (14/25);  new base C (2/2);  total (16/27)
Fisher's exact test for biased strand distribution p-value = 6.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

GTGCGAATGGGCATTTTGGTCTATCCGGCGGAGAGTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATCGCACTTTCAGGTGACTGCCTTTGGATAAGTTGTATCCAGCCTGTTCGTCTG  >  NC_002947/1274945‑1275242
                                                                                                                                                    |                                                                                                                                                     
gTGCGAATGGGCATTTTGGTCTATCCGGCGGAGAGTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAGGGTGAAAGGAGAATGAGCAAAGCCCGCCGATgg                                                                                                                                                      <  1:885009/150‑1 (MQ=255)
     aaTGGGCATTTTGGTCTATCCGGCGGAGAGTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATgggcggg                                                                                                                                                 <  2:977163/150‑1 (MQ=255)
           cATTTTGGTCTATCCGGCGGAGAGTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAg                                                                                                                                           <  2:977951/150‑1 (MQ=255)
            aTTTTGGTCTATCCGGCGGAGAGTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCCATCGGCGGGCtttgct                                                                                                                                          >  2:698244/1‑146 (MQ=255)
            aTTTTGGTCTATCCGGCGGAGAGTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCCATCGGCGGGCtttgct                                                                                                                                          <  1:698244/150‑5 (MQ=255)
                 ggTCTATCCGGCGGAGAGTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAggggggg                                                                                                                                     <  2:455907/150‑1 (MQ=255)
                       tCCGGCGGAGAGTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGaa                                                                                                                               >  1:668755/1‑150 (MQ=255)
                           gcggAGAGTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGc                                                                                                                           <  2:668755/150‑1 (MQ=255)
                               agagTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATc                                                                                                                       <  2:198957/150‑1 (MQ=255)
                                gagTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCg                                                                                                                      >  2:306684/1‑150 (MQ=255)
                                  gTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGta                                                                                                                    <  1:295133/150‑1 (MQ=255)
                                              gCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATgggcgg                                                                                                                                                  >  2:261756/1‑108 (MQ=255)
                                              gCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATgggcgg                                                                                                                                                  <  1:261756/108‑1 (MQ=255)
                                                 tGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTtctc                                                                                                     <  2:624077/150‑1 (MQ=255)
                                                 tGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTtctc                                                                                                     <  2:1099217/150‑1 (MQ=255)
                                                 tGCTTCCCGGTTACGTCGGCGCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTtctc                                                                                                     <  2:1178022/150‑1 (MQ=255)
                                                       ccGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGTTGGGGGTGATGGTGGTAAAGGTGAAAGGAGAAGGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCAc                                                                                               <  2:250691/150‑1 (MQ=255)
                                                           ttACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGaa                                                                                           <  2:708386/150‑1 (MQ=255)
                                                            tACGTAGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTTATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGaaa                                                                                          <  1:974262/150‑1 (MQ=255)
                                                                tCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAg                                                                                                                                           <  1:796943/97‑1 (MQ=255)
                                                                tCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAg                                                                                                                                           >  2:796943/1‑97 (MQ=255)
                                                                tCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTa                                                                                      <  2:371237/150‑1 (MQ=255)
                                                                acGGCTACCATAAGGTGGGGGAGAGCGGTAATGAGGGGGGTGCGGTTGGTCAAGGTGAAAGGAGCGTGAGCAAAGCCCGCCGATGGGCGGGCTTCATGGGGGGCCGGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTa                                                                                      >  2:74493/2‑150 (MQ=255)
                                                                     tCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTtctc                                                                                                     <  1:776761/130‑1 (MQ=255)
                                                                     tCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTtctc                                                                                                     >  2:776761/1‑130 (MQ=255)
                                                                     tCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAg                                                                                 <  1:344782/150‑1 (MQ=255)
                                                                       cAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAggg                                                                                                                            >  2:110898/1‑105 (MQ=255)
                                                                       cAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAggg                                                                                                                            <  1:110898/105‑1 (MQ=255)
                                                                       cAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTcc                                                                                                 <  2:505968/132‑1 (MQ=255)
                                                                       cAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTcc                                                                                                 >  1:505968/1‑132 (MQ=255)
                                                                        aaTAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCAc                                                                              >  1:972949/1‑150 (MQ=255)
                                                                            aGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAggggggg                                                                                                                                     <  2:96895/91‑1 (MQ=255)
                                                                            aGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAggggggg                                                                                                                                     >  1:96895/1‑91 (MQ=255)
                                                                             ggTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCAcccgcc                                                                         >  1:1075933/1‑150 (MQ=255)
                                                                             ggTGAGGGAGAGCGTTAATGATGGGGGTACTGGTGGTAAATGTGAAAGGGGAATGACAAATGCCCGCCAATGGGCGGGTTTTGTGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAGAGGTAGGAAGCAaccgcc                                                                         <  2:16058/150‑1 (MQ=255)
                                                                                     agagCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCActcc                                                                 <  2:972949/150‑1 (MQ=255)
                                                                                           gTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCtcatca                                                           <  1:306684/150‑1 (MQ=255)
                                                                                                     gggTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATgg                                                                                                                                                      >  2:296270/1‑49 (MQ=255)
                                                                                                     gggTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATgg                                                                                                                                                      <  1:296270/49‑1 (MQ=255)
                                                                                                           tggtggTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTCCTACGCTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGGCCCCAGGCTTCATCAGGAATCGCACTTTCAg                                           <  1:634532/150‑1 (MQ=255)
                                                                                                             gtggtAAAGGTGAAAGGAGAATGAGCAAAGACCGCCGATGTGCGGGCTTCAGGGGGGGGCAGAAGGACTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATCGCACTTTCAGGt                                         <  2:1001659/150‑1 (MQ=255)
                                                                                                                      gTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATCGCACTTTCAGGTGACTGCCtt                                >  2:1024516/1‑150 (MQ=255)
                                                                                                                       tGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATCGCACTTTCAGGTGACTGCCttt                               <  2:785161/150‑1 (MQ=255)
                                                                                                                                aaTGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGc                                                                                <  2:1060385/92‑1 (MQ=255)
                                                                                                                                aaTGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGc                                                                                >  1:1060385/1‑92 (MQ=255)
                                                                                                                                        aaGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTCTCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATCGCACTTTCAGGTGACTGCCTTTGGATAAGTTGTATCCAg              >  2:505975/1‑150 (MQ=255)
                                                                                                                                          gcccgcccATCGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATc                                                     >  2:208431/1‑109 (MQ=255)
                                                                                                                                          gcccgcccATCGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATc                                                     <  1:208431/109‑1 (MQ=255)
                                                                                                                                            ccgccgATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATCGCACTTTCAGGTGACTGCCTTTGGCTAAGTTGTATCCAGCCTg          >  2:1174071/1‑150 (MQ=255)
                                                                                                                                                    gggcgggcTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATCGCACTTTCAGGTGACTGCCTTTGGATAAGTTGTATCCAGCCTGTTCGTCTg  >  1:442464/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTGCGAATGGGCATTTTGGTCTATCCGGCGGAGAGTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATCGCACTTTCAGGTGACTGCCTTTGGATAAGTTGTATCCAGCCTGTTCGTCTG  >  NC_002947/1274945‑1275242

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: