Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,022,351 G→C 13.8% intergenic (‑2/+666) PP_3547 ← / ← emrB NAD(P)‑binding oxidoreductase/multidrug transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,022,3510GC13.8% 73.2 / 6.3 29intergenic (‑2/+666)PP_3547/emrBNAD(P)‑binding oxidoreductase/multidrug transporter permease
Reads supporting (aligned to +/- strand):  ref base G (13/12);  new base C (2/2);  total (15/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.82e-01

TGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACA  >  NC_002947/4022207‑4022499
                                                                                                                                                |                                                                                                                                                    
tGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGcc                                                                                                                                                 >  2:228516/1‑150 (MQ=255)
 gCGTTGGAATGACGAACATCAGCACACCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATTGGGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGccc                                                                                                                                                <  1:336758/150‑1 (MQ=255)
        aaTGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAgc                                                                                                                                         >  2:724589/1‑150 (MQ=255)
                     aGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGTCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAg                                                                                                                            <  1:659484/150‑1 (MQ=255)
                           gCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATggg                                                                                                                      <  2:15279/150‑1 (MQ=255)
                                  ctctGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCtca                                                                                                                <  2:698423/149‑3 (MQ=255)
                                  ctctGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCtca                                                                                                                >  1:698423/1‑147 (MQ=255)
                                   tctGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAg                                                                                                                                   >  2:269138/1‑129 (MQ=255)
                                   tctGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAg                                                                                                                                   <  1:269138/129‑1 (MQ=255)
                                             atctggcttttgatctggcttttgATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTg                                                                                                    >  2:376717/1‑150 (MQ=255)
                                                          tctggcttttgatctTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTg                                                                                                    >  1:288911/1‑137 (MQ=37)
                                                          tctggcttttgatctTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTg                                                                                                    <  2:288911/137‑1 (MQ=37)
                                                              gcttttgatctTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGa                                                                                              <  2:723395/139‑1 (MQ=33)
                                                              gcttttgatctTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGa                                                                                              >  1:723395/1‑139 (MQ=33)
                                                              gcttttgatctTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGa                                                                                   <  1:590711/150‑1 (MQ=32)
                                                                            tcttgatcttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCg                                                                     >  2:553037/1‑150 (MQ=32)
                                                                                    ttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGagcgcag                                                             >  1:411201/1‑150 (MQ=32)
                                                                                                                cGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGc                                 <  1:724589/150‑1 (MQ=32)
                                                                                                                        ttGCGGAGGCAGGTGACGGGCATGGATGCCCATCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTg                         >  1:435019/1‑150 (MQ=17)
                                                                                                                          gCGGAGGCAGATGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGgatga                                             >  2:92392/1‑128 (MQ=18)
                                                                                                                          gCGGAGGCAGATGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGgatga                                             <  1:92392/128‑1 (MQ=18)
                                                                                                                                  aGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACttttt               <  2:403928/150‑1 (MQ=32)
                                                                                                                                    ttGACGGGCATCCATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACttttt               >  2:616868/2‑148 (MQ=12)
                                                                                                                                    ttGACGGGCATCCATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACttttt               >  2:616880/2‑148 (MQ=12)
                                                                                                                                    ttGACGGGCATCCATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACttttt               <  1:616880/147‑1 (MQ=12)
                                                                                                                                    ttGACGGGCATCCATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACttttt               <  1:616868/147‑1 (MQ=12)
                                                                                                                                    gTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCgg                                                 <  2:377910/114‑1 (MQ=33)
                                                                                                                                    gTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCgg                                                 >  1:377910/1‑114 (MQ=33)
                                                                                                                                         gggCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGt        >  2:404294/1‑150 (MQ=32)
                                                                                                                                               ggATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACa  >  2:219649/1‑150 (MQ=32)
                                                                                                                                                |                                                                                                                                                    
TGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACA  >  NC_002947/4022207‑4022499

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