Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,670,581 G→A G295R (GGA→AGA)  ynfM → putative arabinose efflux transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,670,5810GA96.8% 190.4 / ‑2.0 62G295R (GGA→AGA) ynfMputative arabinose efflux transporter
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (30/30);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GTCGTTCGTCACGCTGTTTAATTACATCGGCTATCGGTTGATGCTCTCCCCCTGGCATGTCAGTCAGGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCGGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGTTTTCGACGGGGGTTATGCTGTTTGGTTTACTGATGACCTTATTCAGCTCGCTGTGGCTGATCTTTGCCGGAATGTTACTCTTCT  >  NC_000913/1670454‑1670707
                                                                                                                               |                                                                                                                              
gTCGTTCGTCACGCTGTTTAATTACATCGGCTATCGGTTGATGCTCTCCCCCTGGCATGTCAGTCAGGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGacc                                                                                                                     >  2:360221/1‑139 (MQ=255)
  cgttcgtCACGCTGTTTAATTACATCGGCTATCGGTTGATGCTCTCCCCCTGGCATGTCAGTCAGGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGaccac                                                                                                                   <  2:735669/139‑1 (MQ=255)
      cgtCACGCTGTTTAATTACATCGGCTATCGGTTGATGCTCTCCCCCTGGCATGTCAGTCAGGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGc                                                                                                               >  1:839461/1‑139 (MQ=255)
        tCACGCTGTTTAATTACATCGGCTATCGGTTGATGCTCTCCCCCTGGCATGTCAGTCAGGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTa                                                                                                             <  1:1081365/139‑1 (MQ=255)
             cTGTTTAATTACATCGGCTATCGGTTGATGCTCTCCCCCTGGCATGTCAGTCAGGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGc                                                                                                        <  2:839461/139‑1 (MQ=255)
              tGTTTAATTACATCGGCTATCGGTTGATGCTCTCCCCCTGGCATGTCAGTCAGGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCg                                                                                                       <  2:483693/139‑1 (MQ=255)
              tGTTTAATTACATCGGCTATCGGTTGATGCTCTCCCCCTGGCATGTCAGTCAGGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCg                                                                                                       <  2:239060/139‑1 (MQ=255)
              tGTTTAATTACATCGGCTATCGGTTGATGCTCTCCCCCTGGCATGTCAGTCAGGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCg                                                                                                       <  2:420004/139‑1 (MQ=255)
              tGTTTAATTACATCGGCTATCGGTTGATGCTCTCCCCCTGGCATGTCAGTCAGGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCg                                                                                                       <  1:297335/139‑1 (MQ=255)
                          tCGGCTATCGGTTGATGCTCTCCCCCTGGCATGTCAGTCAGGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGaccacc                                                                                                                  <  2:427673/116‑1 (MQ=255)
                          tCGGCTATCGGTTGATGCTCTCCCCCTGGCATGTCAGTCAGGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGaccacc                                                                                                                  >  1:427673/1‑116 (MQ=255)
                          tCGGCTATCGGTTGATGCTCTCCCCCTGGCATGTCAGTCAGGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGacca                                                                                                                    >  2:591806/1‑114 (MQ=255)
                          tCGGCTATCGGTTGATGCTCTCCCCCTGGCATGTCAGTCAGGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGacca                                                                                                                    <  1:591806/114‑1 (MQ=255)
                          tCGGCTATCGGTTGATGCTCTCCCCCTGGCATGTCAGTCAGGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCg                                                                                                       >  1:942649/1‑127 (MQ=255)
                          tCGGCTATCGGTTGATGCTCTCCCCCTGGCATGTCAGTCAGGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCg                                                                                                       <  2:942649/127‑1 (MQ=255)
                          tCGGCTATCGGTTGATGCTCTCCCCCTGGCATGTCAGTCAGGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGAt                                                                                           >  1:281989/1‑139 (MQ=255)
                          tCGGCTATCGGTTGATGCTCTCCCCCTGGCATGTCAGTCAGGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGAt                                                                                           >  2:53459/1‑139 (MQ=255)
                                        aTGCTCTCCCCCTGGCATGTCAGTCAGGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCg                                                                                                       >  2:1069854/1‑113 (MQ=255)
                                        aTGCTCTCCCCCTGGCATGTCAGTCAGGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCg                                                                                                       <  1:1069854/113‑1 (MQ=255)
                                               cccccTGGCATGTCAGTCAGGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGTTTTCGACGGGGGTTAt                                                                      >  2:520981/1‑139 (MQ=255)
                                                   cTGGCATGTCAGTCAGGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGaccacc                                                                                                                  <  1:177285/91‑1 (MQ=255)
                                                   cTGGCATGTCAGTCAGGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGaccacc                                                                                                                  >  2:177285/1‑91 (MQ=255)
                                                     ggCATGTCAGTCAGGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGTTTTCGACGGGGGTTATGCt                                                                   <  1:154347/136‑1 (MQ=255)
                                                     ggCATGTCAGTCAGGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGTTTTCGACGGGGGTTATGCt                                                                   >  2:154347/1‑136 (MQ=255)
                                                          gtcagtcaGGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCg                                                                                                       <  1:744573/95‑1 (MQ=255)
                                                          gtcagtcaGGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCg                                                                                                       >  2:744573/1‑95 (MQ=255)
                                                          gtcagtcaGGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGTTTTCGAc                                                                               >  1:212952/1‑119 (MQ=255)
                                                          gtcagtcaGGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGTTTTCGAc                                                                               <  2:990612/119‑1 (MQ=255)
                                                          gtcagtcaGGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGTTTTCGAc                                                                               >  1:990612/1‑119 (MQ=255)
                                                          gtcagtcaGGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGTTTTCGAc                                                                               <  2:212952/119‑1 (MQ=255)
                                                                cagGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGTTTTCGACGGGGGTTATGCTGTTTGGTTTACtga                                                     <  1:520981/139‑1 (MQ=255)
                                                                    ccGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGTTGTTGTTTTCGACGGGGGTTATGCTGTTTGGTTTACTGATGAc                                                 >  1:322783/1‑139 (MQ=255)
                                                                     cGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGttt                                                                                    <  1:535450/103‑1 (MQ=255)
                                                                     cGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGttt                                                                                    >  2:535450/1‑103 (MQ=255)
                                                                             gCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGTTTTCGACGGGGGTTATGCTGTTTGGTTTACTGATGACCTTATTCAg                                        >  1:301508/1‑139 (MQ=255)
                                                                                       gCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGTTTTCGACGGGGGTTATGCTGTTTGGTTTACTGATGACCTTATTCAGCTCGCTGTgg                              >  1:353815/1‑139 (MQ=255)
                                                                                                    aCCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGTTTTCGACGGGGGTTATGCTGTTTGGTTTACTGATGACCTTATTCAGCTCGCTGTGGCTGATCTTTGCCg                 <  2:281989/139‑1 (MQ=255)
                                                                                                         tACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGTTTTCGACGGGGGTTATGCTGTTTGGTTTACTGATGACCTTATTCAGCTCGCTGTGGCTGATCTTTGCCg                 >  1:402134/1‑134 (MQ=255)
                                                                                                         tACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGTTTTCGACGGGGGTTATGCTGTTTGGTTTACTGATGACCTTATTCAGCTCGCTGTGGCTGATCTTTGCCg                 <  2:402134/134‑1 (MQ=255)
                                                                                                              ggAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTc                                                                                                  <  2:19727/48‑1 (MQ=255)
                                                                                                              ggAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTc                                                                                                  >  1:19727/1‑48 (MQ=255)
                                                                                                                aGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGTTTTCGACGGGGGTTATGCTGTTTGGtt                                                           >  1:880308/1‑85 (MQ=255)
                                                                                                                aGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGTTTTCGACGGGGGTTATGCTGTTTGGtt                                                           <  2:880308/85‑1 (MQ=255)
                                                                                                                  cTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCg                                                                                                       >  1:434534/1‑39 (MQ=255)
                                                                                                                  cTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCg                                                                                                       <  2:434534/39‑1 (MQ=255)
                                                                                                                  cTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGTTTTCGACGGGGGTTATGCTGTTTGGt                                                            >  2:603170/1‑82 (MQ=255)
                                                                                                                  cTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGTTTTCGACGGGGGTTATGCTGTTTGGt                                                            <  1:603170/82‑1 (MQ=255)
                                                                                                                  cTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGTTTTCGACGGGGGTTATGCTGTTTGGTTTACTGATGAc                                                 <  1:522432/93‑1 (MQ=255)
                                                                                                                  cTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGTTTTCGACGGGGGTTATGCTGTTTGGTTTACTGATGAc                                                 >  2:522432/1‑93 (MQ=255)
                                                                                                                  cTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGTTTTCGACGGGGGTTATGCTGTTTGGTTTACTGATGACCTTATTCAGCTCGCTGTGGCTGATCTTTGc                   >  2:363056/1‑123 (MQ=255)
                                                                                                                  cTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGTTTTCGACGGGGGTTATGCTGTTTGGTTTACTGATGACCTTATTCAGCTCGCTGTGGCTGATCTTTGc                   <  1:363056/123‑1 (MQ=255)
                                                                                                                  cTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGTTTTCGACGGGGGTTATGCTGTTTGGTTTACTGATGACCTTATTCAGCTCGCTGTGGCTGATCTTTGCCg                 >  1:901273/1‑125 (MQ=255)
                                                                                                                  cTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGTTTTCGACGGGGGTTATGCTGTTTGGTTTACTGATGACCTTATTCAGCTCGCTGTGGCTGATCTTTGCCg                 <  2:901273/125‑1 (MQ=255)
                                                                                                                  cTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGTTTTCGACGGGGGTTATGCTGTTTGGTTTACTGATGACCTTATTCAGCTCGCTGTGGCTGATCTTTGCCGGAATGTTAct       <  2:158186/135‑1 (MQ=255)
                                                                                                                  cTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGTTTTCGACGGGGGTTATGCTGTTTGGTTTACTGATGACCTTATTCAGCTCGCTGTGGCTGATCTTTGCCGGAATGTTAct       >  1:158186/1‑135 (MQ=255)
                                                                                                                   tCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCg                                                                                                       >  2:623338/1‑38 (MQ=255)
                                                                                                                   tCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCg                                                                                                       <  1:623338/38‑1 (MQ=255)
                                                                                                                   tCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGTTTTCGACGGGGGTTATGCTGTTTGGTTTACTGATGACCTTATTCAGCTCGCTGTGGCTGATCTTTGCCGGAATGTTACtcttct  >  1:441021/1‑139 (MQ=255)
                                                                                                                     aCCCAAAGCCGGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGTTTTCGACGGGGGTTATGCTGTTTGGtt                                                           >  2:860913/1‑80 (MQ=255)
                                                                                                                     aCCCAAAGCCGGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGTTTTCGACGGGGGTTATGCTGTTTGGtt                                                           <  1:860913/80‑1 (MQ=255)
                                                                                                                          aaGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATgttgtt                                                                                     >  2:773377/1‑49 (MQ=255)
                                                                                                                          aaGCCAGAACCATGACCACACGCTATGGGCGTGGTCCAGTGATgttgtt                                                                                     <  1:773377/49‑1 (MQ=255)
                                                                                                                               |                                                                                                                              
GTCGTTCGTCACGCTGTTTAATTACATCGGCTATCGGTTGATGCTCTCCCCCTGGCATGTCAGTCAGGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCGGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGTTTTCGACGGGGGTTATGCTGTTTGGTTTACTGATGACCTTATTCAGCTCGCTGTGGCTGATCTTTGCCGGAATGTTACTCTTCT  >  NC_000913/1670454‑1670707

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