Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 747,140 C→A P210P (CCG→CCT abrB ← putative regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913747,1400CA96.4% 89.6 / ‑2.8 28P210P (CCG→CCTabrBputative regulator
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (18/9);  total (19/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.07e-01

TCCGCCAGCCAATTGCCATATACGCCATCGCCAGCAGCCATTCCGGTAGTTCGATGGTGATGAGCTGACCAGACTGGAGCACTGCGCCCGCCAGCATCGGGATGAGCATCGTGCCGGAAGGAAGACGCAGCAGGCATCCCACTGTACCGGCAACGACCGCGAGCAAAATAGTAAGCAGGAGATTAATGCTCACTGGCGGGAACCAGACGATATGCTGGTTAA  >  NC_000913/747043‑747264
                                                                                                 |                                                                                                                            
tCCGCCAGCCAATTGCCATATACGCCATCGCCAGCAGCCATTCCGGTAGTTCGATGGTGATGAGCTGACCAGACTGGAGCACTGCGCCCGCCAGCATAGGGATGAGCATCGTGCCGGAAGGAAGACGCAGCAGGCATccc                                                                                    >  2:694095/1‑140 (MQ=255)
      agccaATTGCCATATACGCCATCGCCAGCAGCCATTCCGGTAGTTCGATGGTGATGAGCTGACCAGACTGGAGCACTGCGCCCGCCAGCATAGGGATGAGCATCGTGCCGGAAGGAAGACGCAGCAGGCATCCCACTGTa                                                                              >  2:676203/1‑140 (MQ=255)
               ccATATACACCATCGCCAGCAGCCATTCCGGTAGTTCGATGGTGATGAGCTGACCAGACTGGAGCACTGCGCCCGCCAGCATAGGGATGAGCATCGTGCCGGAAGGAAGACGCAGCAGGCATCCCACTGTACCGGCAACg                                                                     >  1:112281/1‑140 (MQ=255)
                  tataCGCCATCGCCAGCAGCCATTCCGGTAGTTCGATGGTGATGAGCTGACCAGACTGGAGCACTGCGCCCGCCAGCATAGGGATGAGCATCGTGCCGGAAGGAAGACGCAGCAGGCATCCCACTGTACCGGCAACGAcc                                                                  >  1:522084/1‑140 (MQ=255)
                    taCGCCATCGCCAGCAGCCATTCCGGTAGTTCGATGGTGATGAGCTGACCAGACTGGAGCACTGCGCCCGCCAGCATAGGGATGAGCATCGTGCCGGAAGGAAGACGCAGCAGGCATCCCACTGTACCGGCAACGACcgc                                                                >  1:560949/1‑140 (MQ=255)
                       gCCATCGCCAGCAGCCATTCCGGTAGTTCGATGGTGATGAGCTGACCAGACTGGAGCACTGCGCCCGCCAGCATAGGGATGAGCATCGTGCCGGAAGGAAGACGCAGCAGGCATCCCACTGTACCGGCAACGACCGCGAg                                                             >  1:8888/1‑140 (MQ=255)
                         cATCGCCAGCAGCCATTCCGGTAGTTCGATGGTGATGAGCTGACCAGACTGGAGCACTGCGCCCGCCAGCATCGGGATGAGCATCGTGCCGGAAGGAAGACGCAGCAGGCATCCCACTGTACCGGCAACGACCGCGAGCa                                                           >  2:568497/1‑140 (MQ=255)
                         cATCGCCAGCAGCCATTCCGGTAGTTCGATGGTGATGAGCTGACCAGACTGGAGCACTGCGCCCGCCAGCATAGGGATGAGCATCGTGCCGGAAGGAAGACGCAGCAGGCATCCCACTGTACCGGCAACGACCGCGAGCa                                                           >  1:476759/1‑140 (MQ=255)
                         cATCGCCAGCAGCCATTCCGGTAGTTCGATGGTGATGAGCTGACCAGACTGGAGCACTGCGCCCGCCAGCATAGGGATGAGCATCGTGCCGGAAGGAAGACGCAGCAGGCATCCCACTGTACCGGCAACGACCGCGAGCa                                                           >  2:22517/1‑140 (MQ=255)
                         cATCGCCAGCAGCCATTCCGGTAGTTCGATGGTGATGAGCTGACCAGACTGGAGCACTGCGCCCGCCAGCATAGGGATGAGCATCGTGCCGGAAGGAAGACGCAGCAGGCATCCCACTGTACCGGCAACGACCGCGAGCa                                                           >  2:26739/1‑140 (MQ=255)
                         cATCGCCAGCAGCCATTCCGGTAGTTCGATGTGATGAGCTGACCAGACTGGAGCACTGCGCCCGCCAGCATAGGGATGAGCATCGTGCCGGAAGGAAGACGCAGCAGGCATCCCACTGTACCGGCAACGACCGCGAGCaa                                                          >  1:337103/1‑140 (MQ=255)
                            cGCCAGCAGCCATTCCGGTAGTTCGATGGTGATGAGCTGACCAGACTGGAGCACTGCGCCCGCCAGCATAGGGATGAGCATCGTGCCGGAAGGAAGACGCAGCAGGCATCCCACTGTACCGGCAACGACCGCGAGCaaaa                                                        >  1:545922/1‑140 (MQ=255)
                                     ccATTCCGGTAGTTCGATGGTGATGAGCTGACCAGACTGGAGCACTGCGCCCGCCAGCATAGGGATGAGCATCGTGCCGGAAGGAAGACGCAGCAGGCATCCCACTGTACCGGCAACGACCGCGAGCAAAATAGTAAGCa                                               >  2:615099/1‑140 (MQ=255)
                                        ttCCGGTAGTTCGATGGTGATGAGCTGACCAGACTGGAGCACTGCGCCCGCCAGCATAGGGATGAGCATCGTGCCGGAAGGAAGACGCAGCAGGCATCCCACTGTACCGGCAACGACCGCGAGCAAAATAGTAAGCAGga                                            <  1:22517/140‑1 (MQ=255)
                                             gTAGTTCGATGGTGATGAGCTGACCAGACTGGAGCACTGCGCCCGCCAGCATAGGGATGAGCATCGTGCCGGAAGGAAGACGCAGCAGGCATCCCACTGTACCGGCAACGACCGCGAGCAAAATAGTAAGCAGGAGATTa                                       >  1:295593/1‑140 (MQ=255)
                                               aGTTCGATGGTGATGAGCTGACCAGACTGGAGCACTGCGCCCGCCAGCATAGGGATGAGCATCGTGCCGGAAGGAAGACGCAGCAGGCATCCCACTGTACCGGCAACGACCGCGAGCAAAATAGTAAGCAGGAGATTAAt                                     >  1:374947/1‑140 (MQ=255)
                                                gTTCGATGGTGATGAGCTGACCAGACTGGAGCACTGCGCCCGCCAGCATAGGGATGAGCATCGTGCCGGAAGGAAGACGCAGCAGGCATCCCACTGTACCGGCAACGACCGCGAGCAAAATAGTAAGCAGGAGATTAATg                                    <  1:26739/140‑1 (MQ=255)
                                                    gATGGTGATGAGCTGACCAGACTGGAGCACTGCGCCCGCCAGCATAGGGATGAGCATCGTGCCGGAAGGAAGACGCAGCAGGCATCCCACTGTACCGGCAACGACCGCGAGCAAAATAGTAAGCAGGAGATTAATGCTCa                                <  2:522084/140‑1 (MQ=255)
                                                     aTGGTGATGAGCTGACCAGACTGGAGCACTGCGCCCGCCAGCATAGGGATGAGCATCGTGCCGGAAGGAAGACGCAGCAGGCATCCCACTGTACCGGCAACGACCGCGAGCAAAATAGTAAGCAGGAGATTAATg                                    >  1:480721/1‑135 (MQ=255)
                                                     aTGGTGATGAGCTGACCAGACTGGAGCACTGCGCCCGCCAGCATAGGGATGAGCATCGTGCCGGAAGGAAGACGCAGCAGGCATCCCACTGTACCGGCAACGACCGCGAGCAAAATAGTAAGCAGGAGATTAATg                                    <  2:480721/135‑1 (MQ=255)
                                                        gTGATGAGCTGACCAGACTGGAGCACTGCGCCCGCCAGCATAGGGATGAGCATCGTGCCGGAAGGAAGACGCAGCAGGCATCCCACTGTACCGGCAACGACCGCGAGCAAAATAGTAAGCAGGAGATTAATGCTCACTgg                            >  2:393579/1‑140 (MQ=255)
                                                            tgaGCTGACCAGACTGGAGCACTGCGCCCGCCAGCATAGGGATGAGCATCGTGCCGGAAGGAAGACGCAGCAGGCATCCCACTGTACCGGCAACGACCGCGAGCAAAATAGTAAGCAGGAGATTAATGCTCACTGGCggg                        <  2:295593/140‑1 (MQ=255)
                                                               gCTGACCAGACTGGAGCACTGCGCCCGCCAGCATAGGGATGAGCATCGTGCCGGAAGGAAGACGCAGCAGGCATCCCACTGTACCGGCAACGACCGCGAGCAAAATAGTAAGCAGGAGATTAATGCTCACTGGCGGGAAc                     <  2:337103/140‑1 (MQ=255)
                                                                cTGACCGGGCTGGAGCACTGCGCCCGCCAGCATAGGGATGAGCATCGTGCCGGAAGGAAGACGCAGCAGGCATCCCACTGTACCGGCAACGACCGCGAGCAAAATAGTAAGCAGGAGATTAATGCTCACTGGCGGGAAcc                    <  1:694095/140‑1 (MQ=255)
                                                                   aCCAGACTGGAGCACTGCGCCCGCCAGCATAGGGATGAGCATCGTGCCGGAAGGAAGACGCAGCAGGCATCCCACTGTACCGGCAACGACCGCGAGCAAAATAGTAAGCAGGAGATTAATGCTCACTgg                            <  2:416289/129‑1 (MQ=255)
                                                                   aCCAGACTGGAGCACTGCGCCCGCCAGCATAGGGATGAGCATCGTGCCGGAAGGAAGACGCAGCAGGCATCCCACTGTACCGGCAACGACCGCGAGCAAAATAGTAAGCAGGAGATTAATGCTCACTgg                            >  1:416289/1‑129 (MQ=255)
                                                                   aCCAGACTGGAGCACTGCGCCCGCCAGCATAGGGATGAGCATCGTGCCGGAAGGAAGACGCAGCAGGCATCCCACTGTACCGGCAACGACCGCGAGCAAAATAGTAAGCAGGAGATTAATGCTCACTGGCGGGAACCAGa                 <  2:476759/140‑1 (MQ=255)
                                                                        aCTGGAGCACTGCGCCCGCCAGCATAGGGATGAGCATCGTGCCGGAAGGAAGACGCAGCAGGCATCCCACTGTACCGGCAACGACCGCGAGCAAAATAGTAAGCAGGAGATTAATGCTCACTGGCGGGAACCAGACGata            >  1:442806/1‑140 (MQ=255)
                                                                                  tGCGCCCGCCAGCATAGGGATGAGCATCGTGCCGGAAGGAAGACGCAGCAGGCATCCCACTGTACCGGCAACGACCGCGAGCAAAATAGTAAGCAGGAGATTAATGCTCACTGGCGGGAACCAGACGATATGCTGGTTaa  >  2:557152/1‑140 (MQ=255)
                                                                                                 |                                                                                                                            
TCCGCCAGCCAATTGCCATATACGCCATCGCCAGCAGCCATTCCGGTAGTTCGATGGTGATGAGCTGACCAGACTGGAGCACTGCGCCCGCCAGCATCGGGATGAGCATCGTGCCGGAAGGAAGACGCAGCAGGCATCCCACTGTACCGGCAACGACCGCGAGCAAAATAGTAAGCAGGAGATTAATGCTCACTGGCGGGAACCAGACGATATGCTGGTTAA  >  NC_000913/747043‑747264

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