Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,093,223 A→G V72V (GTT→GTC)  rhaM ← L‑rhamnose mutarotase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,093,2230AG98.8% 335.0 / ‑4.8 84V72V (GTT→GTC) rhaML‑rhamnose mutarotase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (27/56);  minor base T (1/0);  total (28/56)
Fisher's exact test for biased strand distribution p-value = 3.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

GGTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAATCT  >  NC_000913/4093177‑4093269
                                              |                                              
ggTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACa                                             <  1:3661659/50‑1 (MQ=255)
 gTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACAt                                            >  1:770532/1‑50 (MQ=255)
  ttAGCGGGCATAACATCGGTCATATATTTCCCCCCACTTCGGCCTACATc                                           >  1:5611926/1‑50 (MQ=19)
    aGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCa                                          <  1:1941667/49‑1 (MQ=255)
    aGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAg                                         >  1:4418896/1‑50 (MQ=255)
     gCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGt                                        >  1:3921286/1‑50 (MQ=255)
     gCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGt                                        <  1:893487/50‑1 (MQ=255)
       gggCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGc                                      >  1:5368852/1‑50 (MQ=255)
       gggCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGc                                      <  1:3264772/50‑1 (MQ=255)
         gCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTg                                    >  1:1896414/1‑50 (MQ=255)
         gCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTg                                    <  1:5725482/50‑1 (MQ=255)
          cATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTgg                                   <  1:4432227/50‑1 (MQ=255)
           aTAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGc                                  <  1:5826766/50‑1 (MQ=255)
           aTAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGc                                  <  1:2406141/50‑1 (MQ=255)
            tAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGc                                  >  1:810738/1‑49 (MQ=255)
            tAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCa                                 <  1:351578/50‑1 (MQ=255)
            tAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCa                                 <  1:5124014/50‑1 (MQ=255)
             aaCATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCaa                                >  1:5156955/1‑50 (MQ=255)
              aCATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAAc                               <  1:2959185/50‑1 (MQ=255)
              aCATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAAc                               <  1:3366616/50‑1 (MQ=255)
               cATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAAcc                              >  1:5497396/1‑50 (MQ=255)
               cATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAAcc                              >  1:264658/1‑50 (MQ=255)
               cATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAAcc                              <  1:5347791/50‑1 (MQ=255)
               cATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAAcc                              <  1:640020/50‑1 (MQ=255)
               cATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAAcc                              <  1:1646614/50‑1 (MQ=255)
                aTCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCg                             >  1:2019748/1‑50 (MQ=255)
                aTCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCg                             >  1:5379792/1‑50 (MQ=255)
                 tCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGc                            <  1:4382695/50‑1 (MQ=255)
                   ggTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAt                          <  1:1990121/50‑1 (MQ=255)
                   ggTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAt                          <  1:1975431/50‑1 (MQ=255)
                    gTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAtt                         <  1:3849651/50‑1 (MQ=255)
                    gTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAtt                         <  1:2196712/50‑1 (MQ=255)
                    gTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAtt                         >  1:932612/1‑50 (MQ=255)
                     tCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                        >  1:5069859/1‑50 (MQ=255)
                     tCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                        <  1:5617134/50‑1 (MQ=255)
                     tCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                        <  1:199294/50‑1 (MQ=255)
                      cATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                        >  1:5049854/1‑49 (MQ=255)
                       atatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTcc                       <  1:5117123/49‑1 (MQ=255)
                       atatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTcc                       <  1:671890/49‑1 (MQ=255)
                       atatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                        >  1:5754084/1‑48 (MQ=255)
                       atatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCa                      <  1:657080/50‑1 (MQ=255)
                       atatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCa                      <  1:3768591/50‑1 (MQ=255)
                       atatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCa                      <  1:3708391/50‑1 (MQ=255)
                       atatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCa                      <  1:2653708/50‑1 (MQ=255)
                        tatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAg                     <  1:4599960/50‑1 (MQ=255)
                        tatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAg                     <  1:4270001/50‑1 (MQ=255)
                         atatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGc                    >  1:1378312/1‑50 (MQ=255)
                         atatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGc                    >  1:1607617/1‑50 (MQ=255)
                         atatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGc                    >  1:1917243/1‑50 (MQ=255)
                          tatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGc                    <  1:635456/49‑1 (MQ=255)
                          tatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGc                    <  1:3081897/49‑1 (MQ=255)
                          tatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGc                    <  1:2426720/49‑1 (MQ=255)
                          tatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCg                   >  1:2322931/1‑50 (MQ=255)
                          tatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCg                   >  1:2485304/1‑50 (MQ=255)
                          tatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCg                   >  1:4001166/1‑50 (MQ=255)
                            tttCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGt                  <  1:5768996/49‑1 (MQ=255)
                              tCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttc                <  1:2660884/49‑1 (MQ=255)
                               ccaccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttct               <  1:4259417/49‑1 (MQ=255)
                                caccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttcttc             <  1:3465994/50‑1 (MQ=255)
                                caccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttcttc             <  1:353246/50‑1 (MQ=255)
                                caccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttcttc             <  1:3275561/50‑1 (MQ=255)
                                caccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt              <  1:5360991/49‑1 (MQ=255)
                                caccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt              <  1:5535799/49‑1 (MQ=255)
                                caccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt              <  1:735368/49‑1 (MQ=255)
                                caccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt              <  1:1649812/49‑1 (MQ=255)
                                caccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt              <  1:211915/49‑1 (MQ=255)
                                caccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt              >  1:312442/1‑49 (MQ=255)
                                caccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt              >  1:1687661/1‑49 (MQ=255)
                                 accaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCa            <  1:2710624/50‑1 (MQ=255)
                                 accaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCa            <  1:1537252/50‑1 (MQ=255)
                                  ccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt              >  1:1527572/1‑47 (MQ=255)
                                   caACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGa          <  1:3872241/50‑1 (MQ=255)
                                    aaCGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAt         >  1:2660157/1‑50 (MQ=255)
                                    aaCGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAt         <  1:5070449/50‑1 (MQ=255)
                                    aaCGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAt         <  1:4133401/50‑1 (MQ=255)
                                     aCGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAtt        <  1:2711949/50‑1 (MQ=255)
                                     aCGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAtt        <  1:5634669/50‑1 (MQ=255)
                                      cGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTc       <  1:1705334/50‑1 (MQ=255)
                                      cGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTc       <  1:1271150/50‑1 (MQ=255)
                                       gTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTc       >  1:2286451/1‑49 (MQ=255)
                                       gTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCa      <  1:1823035/50‑1 (MQ=255)
                                        ttGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCaa     <  1:508825/50‑1 (MQ=255)
                                           gCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAtct  >  1:4926283/1‑50 (MQ=255)
                                           gCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAtc   <  1:3056839/49‑1 (MQ=255)
                                              |                                              
GGTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAATCT  >  NC_000913/4093177‑4093269

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