Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,897,782 G→C 13.6% intergenic (‑182/+52) PP_5172 ← / ← PP_5173 hypothetical protein/RND family transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,897,7820GC13.6% 87.5 / 5.8 36intergenic (‑182/+52)PP_5172/PP_5173hypothetical protein/RND family transporter
Reads supporting (aligned to +/- strand):  ref base G (12/18);  new base C (3/2);  total (16/20)
Fisher's exact test for biased strand distribution p-value = 6.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.10e-01

AGGTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACGATGA  >  NC_002947/5897634‑5897922
                                                                                                                                                    |                                                                                                                                            
agGTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATg                                                                                                                                              <  2:146303/149‑1 (MQ=255)
  gtgtGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCa                                                                                                                                            <  5:86151/146‑1 (MQ=255)
  gTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCa                                                                                                                                            >  1:224766/1‑149 (MQ=255)
       aaGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCAt                                                                                                                                        <  6:144243/148‑1 (MQ=255)
                tGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGACCACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTcgc                                                                                                                               <  6:397261/148‑1 (MQ=255)
                  cgcCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCgg                                                                                                                             <  1:37461/148‑1 (MQ=255)
                  cgcCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCgg                                                                                                                             <  2:171683/148‑1 (MQ=255)
                  cgcCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCgg                                                                                                                             <  5:175455/148‑1 (MQ=255)
                     cTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTCGCgg                                                                                                                             >  4:139398/1‑145 (MQ=25)
                     cTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTCGCgg                                                                                                                             <  3:139398/145‑1 (MQ=25)
                     cTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGAGTCTGGTAGGAGCGGGTTTACACGCGAATACGGTGGTCGATGCAGAATCGCATACGCGGGtc                                                                                                                          >  5:8205/1‑147 (MQ=255)
                        tcTTGAGGGCAATTGTAGTTTGGCGGGATTATTCTTAAAAATACTGTTTAATAATTTTTATATACCCATTTAAGAGAACTCGCCTCTGGTAGGAGCGGGTTTACCGGCGAATACGGTGGTCGATGAACCATCGCATTCGCGGGTAAAcc                                                                                                                      <  6:8205/148‑1 (MQ=38)
                        gCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAAcc                                                                                                                      <  2:98536/149‑1 (MQ=255)
                        gCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAAc                                                                                                                       <  8:70244/148‑1 (MQ=255)
                        gCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAAc                                                                                                                       <  2:143889/148‑1 (MQ=255)
                         cTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAAccc                                                                                                                     <  8:166878/149‑1 (MQ=255)
                                  aaTAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCCCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTCGCGGGTAAAcc                                                                                                                      >  6:194467/1‑139 (MQ=21)
                                  aaTAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCCCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTCGCGGGTAAAcc                                                                                                                      <  5:194467/139‑1 (MQ=21)
                                                       ttCTTAAAAATACTTTTTAAGAATTTTAATATAACCATTAAAGAGAACACGCCTCTGGTAGGAGCGGGTTAACCCGCTAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAg                                                                                       <  2:224766/149‑1 (MQ=38)
                                                          ttAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTAACCGCGAATACGGTGGTCGATGCACCATAGCATTAGCGGGTAAAACCGCTCCTACCAGGGGCCATATTGGCCTTCAGGCg                                                                                    >  3:238294/1‑149 (MQ=255)
                                                                           gAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTGGATTCACCATCGCATTAGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTGAGGCGTGTTTCTGCGGCTcggg                                                                   >  1:216865/1‑147 (MQ=255)
                                                                                        aCCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCAATCGCGGGTAAACCCGCTCCTAACAGGGGCCATATTGGCCTTCAAGCGGTTTTCTTCGGCTCGCCGATCATGTAACaa                                                      >  5:6671/1‑148 (MQ=37)
                                                                                                atataaCAGGCGTTTGGTTGGAGGGTGTTTACCCGCTAATAGGGTGGTCGATGCACCATCGCATTCGCGGTTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTa                                              <  6:6671/145‑1 (MQ=255)
                                                                                                agagAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTa                                              >  3:46386/1‑149 (MQ=255)
                                                                                                 gagaACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTa                                              >  1:80452/1‑148 (MQ=255)
                                                                                                             ctGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAgg                                 >  1:236380/1‑149 (MQ=255)
                                                                                                                gTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGaa                               <  6:378118/148‑1 (MQ=255)
                                                                                                                   ggAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTCGCGGGTAAACCCGCACCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCAGGCAGGAAcag                            >  5:300331/1‑148 (MQ=21)
                                                                                                                      gCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACagcagc                        >  3:126484/1‑149 (MQ=255)
                                                                                                                          gTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAAcag                            <  4:395340/141‑1 (MQ=255)
                                                                                                                          gTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAAcag                            >  3:395340/1‑141 (MQ=255)
                                                                                                                          gTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGt                      >  4:61366/1‑147 (MQ=255)
                                                                                                                             tACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGc                  <  4:238294/148‑1 (MQ=255)
                                                                                                                                ccGCGAATACGGTGGTCGATGCACCATCGCATTTGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAgg              >  4:329386/1‑149 (MQ=255)
                                                                                                                                 cgcgAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGACACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGt             >  2:317082/1‑149 (MQ=255)
                                                                                                                                         cGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACga     <  3:61366/149‑1 (MQ=255)
                                                                                                                                            tggtCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACgatga  <  2:236380/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
AGGTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACGATGA  >  NC_002947/5897634‑5897922

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑