Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,984,734 G→A G112G (GGC→GGT araG ← L‑arabinose ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,984,7340GA96.2% 245.1 / 0.1 80G112G (GGC→GGTaraGL‑arabinose ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (43/32);  total (45/35)
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.83e-01

GAGATATTTCAGCGGCGTGTCCGGGTCAATATCCATACCAAGATGTTTAAGTTGTAAACCCGCCTCATAATTCAGCAATGAGCGATTCACAATGCCGCCTTTATGCGGCAGCTGGCCGAGATAGATGTTTTCCGCGACGGTCATTTCCGGCACGAGATGCAGTTCCTGGTAAATAATCGCCACGCCCGCG  >  NC_000913/1984638‑1984827
                                                                                                |                                                                                             
gagaTATTTCAGCGGCGTGTCCGGGTCAATATCCATACCAAGATGTTTAAGTTGTAAACCCGCCTCATAATTCAGCAATGAGCGATTCACAATGCCACCt                                                                                            <  2:634830/100‑1 (MQ=255)
gagaTATTTCAGCGGCGTGTCCGGGTAAATATCCATACCAAGATGTTTAAGTTGTAAACCCGCCTCATAATTCAGCAATGAACGATTCACAATGCCACCt                                                                                            >  1:1589702/1‑100 (MQ=255)
 agaTATTTCAGCTGCGTGTCCGGGTCAATATCCATACCAACATGTTTAAGTTGTAAACCCGCCTCATAATTCAGCAATGAGCGATTCACAATGCCACCtt                                                                                           <  1:321514/100‑1 (MQ=255)
  gaTATTTCAGCGGCGTGTCCGGGTCAATATCCATACCAAGATGTTTAAGTTGTAAACCCGCCTCTTAATTCAGCAGTGAGCGATTCACAATGGCACCtt                                                                                           <  1:348756/99‑1 (MQ=255)
  gaTATTTCAGCGGCGTGTCCGGGTCAATATCCATACCAAGATGTTTAAGTTGTAAACCCGCCTCATAATTCAGCAATGAGCGATTCACAATGCCACCttt                                                                                          >  1:452295/1‑100 (MQ=255)
   atctttCAGCGGCGTGTCCAAGTCAATATCCATACTAAGATGTTTAAGTAGTAAACCCGCCACATAATTCAGCAATGAGCGATTTACAAAGCCACCTTTa                                                                                         >  1:1243630/4‑100 (MQ=255)
   atatTTCAGCGGCGTGTCCGGGTCAATATCCATACCAAGATGTTTAAGTTGTAAACCCGCCTCATAATTCAGCAATGAGCGATTCACAATGCCACCTTTa                                                                                         <  2:1360196/100‑1 (MQ=255)
     atTTCAGCGGCGTGTCCGTGTCAATATCCATACCAATATTTTTTTTTTTTTATCCCGCCTCATAATTCAGCAATGAGCGATTCACAATGCCACCTTTATg                                                                                       <  2:8509/100‑1 (MQ=255)
     atTTCAGCGGCGTGTCCGGGTCAATATCCATACCAAGATGTTTAAGTTGTAAACCCGCCTCATAATTCAGCAATGAGCGATTCACAATGCCACCTTTATg                                                                                       <  2:955516/100‑1 (MQ=255)
     atTTCAGCGGCGTGTCCGGGTCAATATCCATACCAAGATGTTTAAGTTGTAAACCCGCCTCATAATTCAGCAATGAGCGATTCACAATGCCACCTTTATg                                                                                       >  2:477847/1‑100 (MQ=255)
     atTTCAGCGGCGTGTCCGGGTCAATATCCATACCAAGATGTTTAAGTTGTAAACCCGCCTCATAATTCAGCAATGAGCGATTCACAATGCCACCTTTATg                                                                                       >  1:1445015/1‑100 (MQ=255)
          aGCGGCGTGTCCGGGTCAATATCCATACCAAGATGTTTAAGTTGTAAACCCGCCTCATAATTCAGCAATGAGCGATTCACAATGCCACCTTTATGCGGc                                                                                   <  1:1205299/99‑1 (MQ=255)
             ggcgTGTCCGGGTCAATATCCATACCAAGATGTTTAAGTTGTAAACCCGCCTCATAATTCAGCAATGAGCGATTCACAATGCAACCTTTATGCGGAAGCt                                                                               >  2:1414548/1‑100 (MQ=255)
              gcgTGTCCGGGTCAATATACATACCAAGATGTTTAAGTTGTAAACCCGCCTCCTAATTCAGCAATTAGCGATTCACAATGCAACCTTTATGCGGAAg                                                                                 >  2:736091/1‑97 (MQ=255)
                gtgtCCGTGTCAATATCCATACCAAGATGTTTAAGTTGTAAACCCGCCTCATAATTCAGCAATGAGCGATTCACAATGCCACCTTTATGCGGCAg                                                                                 <  2:1017073/95‑1 (MQ=255)
                gtgtCCGGGTCAATATCCATACCAAGATGTTTAAGTTGTAAACCCGCCTCATAATTCAGCAATGAGCGATTCACAATGCCACCTTTATGCGGCAGCTGGc                                                                            <  2:543635/100‑1 (MQ=255)
                 tgtCCGGGTCAATATCCATACCAAGATGTTTAAGTTGTAAACCCGCCTCATAATTCAGCAATGAGCGATTCACAATGCCACCTTTATGCGGCAGCTGGcc                                                                           >  1:865223/1‑100 (MQ=255)
                  gtCCGGGTCAATATCCATACCAAGATGTTTAAGTTGTAAACCCGCCTCATAATTCAGCAATGAGCGATTCACAATGCCACCTTTATGCGGCAGCTGGCCg                                                                          <  1:981474/100‑1 (MQ=255)
                  gtCCGGGTCAATATCCATACCAAGATGTTTAAGTTGTAAACCCGCCTCATAATTCAGCAATGAGCGATTCACAATGCCACCTTTATGCGGCAGCTGGCCg                                                                          <  1:1168876/100‑1 (MQ=255)
                    ccGGGTCAATATCCATACCAAGATGTTTAAGTTGTAAACCCGCCTCATAATTCAGCCATGAGCGATTCACAATGCCACCTTTATGCGGCAGCTGGCCgag                                                                        >  2:1710095/1‑100 (MQ=255)
                    ccGGGTCAATATCCATACCAAGATGTTTAAGTTGTAAACCCGCCTCATAATTCAGCAATGAGCGATTCACAATGCCACCTTTATGCGGCAGCTGGCCgag                                                                        <  1:405272/100‑1 (MQ=255)
                    ccGGGTCAATATCCATACCAAGATGTTTAAGTTGTAAACCCGCCTCATAATTCAGCAATGAGCGATTCACAATGCCACCTTTATGCGGCAGCTGGCCga                                                                         >  2:355444/1‑99 (MQ=255)
                     cGGGTCAATATCCATACCAAGATGTTTAAGTTGTAAACCCGCCTCATAATTCAGCAATGAGCGATTCACAATGCCACCTTTATGCGGCAGTTGGCCgaga                                                                       <  2:1408328/100‑1 (MQ=255)
                     cGGGTCAATATCCATACCAAGATGTTTAAGTTGTAAACCCGCCTCATAATTCAGCAATGAGCGATTCACAATGCCACCTTTATGCGGCAGCTGGCCgag                                                                        >  1:1493144/1‑99 (MQ=255)
                       ggTCAATATCCATACCAAGATGTTTAAGTTGTAAACCCGCCTCATAATTCAGCAATGAGCGATTCACAATGCCACCTTTATGCGGCAGCTGGCCGagata                                                                     >  1:1585115/1‑100 (MQ=255)
                       ggTCAATATCCATACCAAGATGTTTAAGTTGTAAACCCGCCTCATAATTCAGCAATGAGCGATTCACAATGCCAACTTTATGCGGCAGATgg                                                                             >  1:565763/1‑92 (MQ=255)
                        gTCAATATCCATACCAAGATGTTTAAGTTGTAAACCCGCCTCATAATTCAGCAATGAGCGATTCACAATGCCACCTTTATGCGGCAGCTGGCCGagatag                                                                    <  2:765010/100‑1 (MQ=255)
                          cAATATCCATACCAAGATGTTTAAGTTGTAAACCCGCCTCATAATTCAGCAATGAGCGATTCACAATGCCACCTTTATGCGGCAGCTGGCCGagatagat                                                                  >  1:2037/1‑100 (MQ=255)
                            atatCCATACCAAGATGTTTAAGTTGTAAACCCGCCTCATAATTCAGCAATGAGCGATTCACAATGCCACCTTTATGCGGCAGCTGGCCGAGATAGATGt                                                                >  1:1031716/1‑100 (MQ=255)
                             tatCCATACCAAGATGTTTAAGTTGTAAACCCGCCTCATAATTCAGCAATGAGCGATTCACAATGCCACCTTTATGCGGCAGCTGGCCGAGATAGATGtt                                                               <  2:706580/100‑1 (MQ=255)
                             tatCCATACCAAGATGTTTAAGTTGTAAACCCGCCTCATAATTCAGCAATGAGCGATTCACAATGCCACCTTTATGCGGCAGCTGGCCGAGATAGATGtt                                                               <  2:465804/100‑1 (MQ=255)
                              atCCATACCAAGATGTTTAAGTTGTAAACCCGCCTCATAATTCAGCAATGAGCGATTCACAATGCCACCTTTATGCGGCAGCTGGCCGAGATAGATGttt                                                              >  1:476324/1‑100 (MQ=255)
                               tCCATACCAAGATGTTTAAGTTGTAAACCCGCCTCATAATTCAGCAATGAGCGATTCACACTGCCACCTTTATGCGGCAGCTGGCCGAGATAGATGtttt                                                             >  1:343116/1‑100 (MQ=255)
                                ccATACCAAGATGTTTAAGTTGTAATCCCGCCTCATAATTCAGCAATGAGCGATTCACAATGCCACCTTTATGCGGCAGCTGGCCGAGATAGATGTTTTc                                                            <  1:881731/100‑1 (MQ=255)
                                 cATACCAAGATGTTTAAGTTGTAAACCCGCCTCATAATTCAGCAATTAGCGATTCACAATGCCACCTTTATGCGGCAGCTGGCCGAGATAGATGTTTTcc                                                           >  1:607813/1‑100 (MQ=255)
                                 cATACCAAGATGTTTAAGTTGTAAACCCGCCTCATAATTCAGCAATGAGCGATTCACAATGCCACCTTTATGCGGCAGCTGGCCGAGATAGATGTTTTcc                                                           >  2:651359/1‑100 (MQ=255)
                                   tACCAAGATGTTTAAGTTGTAAACCCGCCTCATAATTCAGCAATGAGCGATTCACAATGCCACCTTTATGCGGCAGCTGGCCGAGATAGATGTTTTCcgc                                                         <  2:764381/100‑1 (MQ=255)
                                     ccAAGATGTTTAAGTTGTAAACCCGCCTCATAATTCAGCAATGAGCGATTCACAATGCCACCTTTATGCGGCAGCTGGCCGAGATAGATGTTTTCCGCGa                                                       >  1:221362/1‑100 (MQ=255)
                                      cAAGATGTTTAAGTTGTAAACCCGCCTCATAATTCAGCAATGAGCGATTCACAATGCCTCCTTTATGCGGCAGCTGGCCGAGATAGATGTGTTCCGCGAc                                                      <  1:1701012/100‑1 (MQ=255)
                                      cAAGATGTTTAAGTTGTAAACCCGCCTCATAATTCAGCAATGAGCGATTCACAATGCCACCTTTATGCGGCAGCTGGCCGAGATAGATGTTTTCCGCGAc                                                      >  2:407993/1‑100 (MQ=255)
                                      cAAGATGTTTAAGTTGTAAACCCGCCTCATAATTCAGCAATGAGCGATTCACAATGCCACCTTTATGCGGCAGCTGGCCGAGATAGATGTTTTCCGCGAc                                                      >  2:408100/1‑100 (MQ=255)
                                      cAAGATGTTTAAGTTGTAAACCCGCCTCATAATTCAGCAATGAGCGATTCACAATGCCACCTTTATGCGGCAGCTGGCCGAGATAGATGTTTTCCGCGAc                                                      >  1:932205/1‑100 (MQ=255)
                                           tGTTTAAGTTGTAAACCCGCCTCATAATTCAGCAATGAGCGATTCACAATGCCGCCTTTATGCGGCAGCTGGCCGAGATAGATGTTTTCCGCGACGGTCa                                                 >  2:1270309/1‑100 (MQ=255)
                                             tttAAGTTGTAAACCCGCCTCATAATTCAGCAATGAGCGATTCACAATGCCGCCTTTATGCGGCAGCTGGCCGAGATAGATGTTTTCCGCGACGGTCAtt                                               >  1:136658/1‑100 (MQ=255)
                                             tttAAGTTGTAAACCCGCCTCATAATTCAGCAATGAGCGATTCACAATGCCACCTTTATGCGGCAGCTGGCCGAGATAGATGTTTTCCGCGACGGTCAtt                                               >  2:43727/1‑100 (MQ=255)
                                               tAAGTTGTAAACCCGCCTCATAATTCAGCAATGAGCGATTCACAATGCCACCTTTATGCGGCAGCTGGCCGAGATAGATGTTTTCTGCGACGGTCATTTc                                             >  1:625646/1‑100 (MQ=255)
                                                aaGTTGTAAACCCGCCTCATAATTCAGCAATAAGCGATTCACAATGCCACCTTTATGCGGCAGCTGGCCGAGATAGATGTTTTCCGCGACGGTCATTTcc                                            >  1:145028/1‑100 (MQ=255)
                                                 aGTTGTAAACCCGCCTCATAATTCAGCAATGAGCGATTCACAATGCCACCTTTATGCGGCAGCTGGCCGAGATAGATGTTTTCCGCGACGGTCATTTCCg                                           >  2:1775673/1‑100 (MQ=255)
                                                  gTTGTAAACCCGCCTCATAATTCAGCAATGAGCGATTCACAATGCCACCTTTATGCGGCAGCTGGCCGAGATAGATGTTTTCCGCGACGGTCATTTCCgg                                          <  2:1461166/100‑1 (MQ=255)
                                                   ttGTAAACCCGCCTCATAATTCAGCAATGAGCGATTCACAATGCCACCTTTATGCGGCAGCTGGCCGAGATAGATGTTTTCCGCGACGGTCATTTCCGGc                                         <  1:475681/100‑1 (MQ=255)
                                                       aaaCCCGCCTTATAATTCATCAATGAGCGATTCACAATGCCACCTTTATGCGGCAGCTGGCCGAGATAGATGTTTTCCGCGACGGTCATTTCCGGCACg                                      >  1:1097614/1‑99 (MQ=255)
                                                       aaaCCCGCCTCATAATTCAGCAATGAGCGATTCACAATGCCGCCTTTATGCGGCAGCTGGCCGAGATAGATGTTTTCCGCGACGGTCATTTCCGGCACga                                     <  1:1270298/100‑1 (MQ=255)
                                                        aaCCCGCCTCATAATTCAGCAATGAGCGATTCACAATGCCACCTTTATGCGGCAGCTGGCCGAGATAGATGTTTTCCGCGACGGTCATTTCCGGCACgag                                    <  2:1031726/100‑1 (MQ=255)
                                                        aaCCCGCCTCATAATTCAGCAATGAGCGATTCACAATGCCACCTTTATGCGGCAGCTGGCCGAGATAGATGTTTTCCGCGACGGTCATTTCCGGCACgag                                    >  1:867549/1‑100 (MQ=255)
                                                            cGCCTCATAATTCAGCAATGAGCGATTCACAATGCCACCTTTATGCGGCAGCTGGCCGAGATAGATGTTTTCCGCGACGGTCATTTCCGGCACGAGATGc                                <  1:1069095/100‑1 (MQ=255)
                                                              ccTCATAATTCAGCAATGAGCGATTCACAATGCCACCTTTATGCGGCAGCTGGCCGAGATAGATGTTTTCCGCGACGGTCATTTCCGGCACgag                                    <  2:637755/94‑1 (MQ=255)
                                                              ccTCATAATTCAGCAATGAGCGATTCACAATGCCACCTTTATGCGGCAGCTGGCCGAGATAGATGTTTTCCGCGACGGTCATTTCCGGCACGAGATGCAg                              <  2:1296727/100‑1 (MQ=255)
                                                                tCATAATTCAGCAATGAGCGATTCACAATGCCACCTTTATGCGGCAGCTGGCCGAGATAGATGTTTTCCGCGACGGTCATTTCCGGCACGAGATGCAGtt                            <  2:867555/100‑1 (MQ=255)
                                                                 cATAATTCAGCAATGAGCGATTCACAATGCCACCTTTATGCGGCAGCTGGCCGAGATAGATGTTTTCCGCGACGGTCATTTCCGGCACGAGATGCAGTTc                           >  1:1737164/1‑100 (MQ=255)
                                                                 cATAATTCAGCAATGAGCGATTCACAATGCCACCTTTATGCGGCAGCTGGCCGAGATAGATGTTTTCCGCGACGGTCATTTCCGGCACGAGATGCAGTTc                           >  1:1263203/1‑100 (MQ=255)
                                                                 cATAATTCAGCAATGAGCGATTCACAATGCCACCTTTATGCGGCAGCTGGCCGAGATAGATGTTTTCCGCGACGGTCATTTCCGGCACGAGATGCAGTTc                           <  2:452297/100‑1 (MQ=255)
                                                                 cATAATTCAGCAATGAGCGATTCACAATGCCACCTTTATGCGGCAGCTGGCCGAGATAGATGTTTTCCGCGACGGTCATTTCCGGCACGAGATGCAGTTc                           >  1:554830/1‑100 (MQ=255)
                                                                   tAATTCAGCAATGAGCGATTCACAATGCCACCTTTATGCGGCAGCTGGCCGAGATAGATGTTTTCCGCGACGGTCATTTCCGGCACGAGATGCAGTTcc                          >  1:1311708/1‑99 (MQ=255)
                                                                   tAATTCAGCAATGAGCGATTCACAATGCCACCTTTATGCGGCAGCTGGCCGAGATAGATGTTTTCCGCGACGGTCATTTCCGGCACGAGATGCAGTTCCt                         >  2:53418/1‑100 (MQ=255)
                                                                    aaTTCAGCAATGAGCGATTCACAATGCCACCTTTATGCGGCAGCTGGCCGAGATAGATGTTTTCCGCGACGGTCATTTCCGGCACGAGATGCAGTTCCt                         >  2:806030/1‑99 (MQ=255)
                                                                    aaTTCAGCAATGAGCGATTCACAATGCCACCGTTATGCGGCAGCTGGCCGAGATAGATGTTTTCCGCGACGGTCATTTCCGGCACGAGATGCAGTTCCTg                        <  1:1224771/100‑1 (MQ=255)
                                                                     aTTCAGCAATGAGCGATTCACAATGCCACCTTTATGCGGCAGCTGGCCGAGATAGATGTTTTCCGCGACGGTCATTTCCGGCACGAGATGCAGTTCCTgg                       >  1:1514794/1‑100 (MQ=255)
                                                                      ttCAGCAATGAGCGATTCACAATGCCACCTTTATGCGGCAGCTGGCCGAGATAGATGTTTTCCGCGACGGTCATTTCCGGCACGATATGCAGTTCCTGGt                      <  1:553319/100‑1 (MQ=255)
                                                                       tCAGCAATGAGCGATTCACAATGCCACCTTTATGCGGCAGCTGGCCGAGATAGATGTTTTCCGCGACGGTCATTTCCGGCACGAGATGCAGTTCCTGGTa                     >  2:1123445/1‑100 (MQ=255)
                                                                       tCAGCAATGAGCGATTCACAATGCCACCTTTATGCGGCAGCTGGCCGAGATAGATGTTTTCCGCGACGGTCATTTCCGGCACGAGATGCAGTTCCTGGTa                     <  2:74130/100‑1 (MQ=255)
                                                                       tCAGCAATGAGCGATTCACAATGCCACCTTTATGCGGCAGCTGGCCGAGATAGATGTTTTCCGCGACGGTCATTTCCGGCACGAGATGCAGTTCCTGGTa                     >  1:199776/1‑100 (MQ=255)
                                                                         aGCAATGAGCGATTCACAATGCCACCTTTATGCGGCAGCTGGCCGAGATAGATGTTTTCCGCGACGGTCATTTCCGGCACGAGATGCAGTTCCTGGTaaa                   >  2:1597845/1‑100 (MQ=255)
                                                                          gCAATGAGCGATTCACAATGCCACCTTTATGCGGCAGCTGGCCGAGATAGATGTTTTCCGCGACGGTCATTTCCGGCACGAGATGCAGTTCCTGGTAaat                  >  2:1303052/1‑100 (MQ=255)
                                                                          gCAATGAGCGATTCACAATGCCACCTTTATGCGGCAGCTGGCCGAGATAGATGTTTTCCGCGACGGTCATTTCCGGCACGAGATGCAGTTCCTGGTAaat                  <  2:1247970/100‑1 (MQ=255)
                                                                              tGATCGATTTACAATGTCACCTTTATGCGGCAGCTGGCCGAGATAGATGTTTTCCGCGACGGTCATTTCCGGCACGAGATGCAGTTCCTGGTAaat                  <  2:1589713/96‑1 (MQ=255)
                                                                              tGAGCGATTAACAATGCCCCCTTTATGCGGCAGCTGGCCGAGATAGATGTTTTCTGCGACGGTCATTTTCGGCACGTGTTGCAGTTCCTGTTAaataa                <  1:1243612/98‑1 (MQ=255)
                                                                               gAGCGATTCACAATGCCACCTTTATGCGGCAGCTGGCCGAGATAGATGTTTTCCGCGACGGTCATTTCCGGCACGAGATGCAGTTCCTGGTAAATAATCg             >  2:208840/1‑100 (MQ=255)
                                                                                  cGATTCACAATGCCACCTTTATGCGGCAGCTGGCCGAGATAGATGTTTTCCGCGACGGTCATTTCCGGCACGAGATGCAGTTCCTGGTAAATAATCGCCa          >  1:357111/1‑100 (MQ=255)
                                                                                   gcttCACAATGCCACCTTTATGCGGAAGCTGGGCGAGATAGATGGTTGCCGCGACGGTCATTTCCGGCACTAGATGCAGTTCCTCGTAAATAATCTCCAc         >  1:148549/3‑100 (MQ=255)
                                                                                          aaTGCCACCTTTATGCGGCAGCTGGCCTAGATAGATGTTTTCCGCGACGGTCATTTCCGGCACGAGATGCAGTTCCTGGTAAATAATCGCCACGCCcgcg  <  2:428171/100‑1 (MQ=255)
                                                                                                |                                                                                             
GAGATATTTCAGCGGCGTGTCCGGGTCAATATCCATACCAAGATGTTTAAGTTGTAAACCCGCCTCATAATTCAGCAATGAGCGATTCACAATGCCGCCTTTATGCGGCAGCTGGCCGAGATAGATGTTTTCCGCGACGGTCATTTCCGGCACGAGATGCAGTTCCTGGTAAATAATCGCCACGCCCGCG  >  NC_000913/1984638‑1984827

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