Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,985,112 C→G 13.2% intergenic (‑43/+27) araG ← / ← araF arabinose ABC transporter ATP binding subunit/arabinose ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,985,1120CG13.2% 131.4 / 8.2 49intergenic (‑43/+27)araG/araFarabinose ABC transporter ATP binding subunit/arabinose ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (22/20);  new base G (5/2);  total (27/22)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

GGCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTACCCAGT  >  NC_000913/1984974‑1985254
                                                                                                                                          |                                                                                                                                              
ggCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGACGCCGCGAAATGACGCATACAGGGTAAACTGTAGCATACGGGCTCCAGGAATCAAGCGTCTTTCACACCCCCCCTCcagcacac                                                                                                                                            >  2:277549/1‑140 (MQ=255)
 gCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACGGTGGCATAACGTCTCCGCGAATCACCGGGCGTTCACACCCCCTCTGCATCATGCa                                                                                                                                           >  2:1093801/1‑143 (MQ=255)
     aGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGCGAATCACTGGTCGTTCACACGCCCTCTGCCTCCTGCACAgg                                                                                                                                       >  2:779467/1‑143 (MQ=255)
      gCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATGATGCAGAggg                                                                                                                                      <  1:311651/143‑1 (MQ=255)
       cAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATGATGCAgggggg                                                                                                                                     >  1:932643/1‑143 (MQ=255)
            aaaaCTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATGATGCAGAgggggt                                                                                                                                   <  2:883865/140‑3 (MQ=255)
            aaaaCTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATGATGCAGAgggggt                                                                                                                                   >  1:883865/1‑138 (MQ=255)
               cgtAATAGCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGATTTGTTc                                                                                                                             <  2:1056241/141‑1 (MQ=255)
               aCTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTGCACACCCCATCTGCATCGTGCAGAGGGGAAATTTTc                                                                                                                             >  2:354699/1‑143 (MQ=255)
                cTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTcc                                                                                                                            <  1:910984/143‑1 (MQ=255)
                  aaTATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCGCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCgg                                                                                                                          <  2:1047068/143‑1 (MQ=255)
                       ccGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGGCTCCGTGAATCACTGGTCGTTCACACCCCCCCTGCATCATGCAGAGGGGAATTTTTCCGGCAAAt                                                                                                                     <  1:47944/143‑1 (MQ=255)
                             gcgcCTTAACGCCGGGAAACGTTTTACCGATGCCGCGACATGAGAGGTACGGAGTAGACTGTTGCACAACGTCTCCGTGAATCACTGGTCGTTCACACCGACTCTGCATCATGCAGCAGGGAAGTTTTCCCGCAGAttaccta                                                                                                               >  2:244113/1‑143 (MQ=255)
                              cgcCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAAttacttac                                                                                                              >  1:538236/1‑143 (MQ=255)
                                    aaCGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAGTCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTa                                                                                                        <  1:193706/143‑1 (MQ=255)
                                          gggAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCACCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAAccat                                                                                                  >  2:281311/1‑141 (MQ=255)
                                           ggAAACGTCTTACCGACGCCGCGAAACGAGCGACACGGCGTAGACCGTTGCACAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCttt                                                                                                 <  1:191264/143‑1 (MQ=255)
                                             aaaCGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCCTCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAACCCttttt                                                                                               >  1:1102637/1‑143 (MQ=255)
                                                         gATGCCTCAAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGttctt                                                                                   <  2:410700/143‑1 (MQ=255)
                                                              cgcgAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTa                                                                              >  2:425401/1‑143 (MQ=255)
                                                               gcgAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTaa                                                                             <  1:62146/143‑1 (MQ=255)
                                                                cgAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCTGCAAATTACTTCCCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTaaa                                                                            >  1:945326/1‑143 (MQ=255)
                                                                      gtgagaTACGGAGTAGGCTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCTCCTCTGCATCATGCAGAGGTGAATTTTTCCGGCAAATTACTTACCTCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTc                                                                      <  2:151705/141‑1 (MQ=255)
                                                                             aCGGGGTAGACTGTTGCATAACGTCTTCGTGAATCGCTGGTCGTGCACCCCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTg                                                               <  2:870831/143‑1 (MQ=255)
                                                                             aCGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCCTCATGCAGAGTGGGATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAATGTCGTCAcgcgcg                                                               >  2:834993/1‑141 (MQ=255)
                                                                                   tAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGt                                                         <  2:40891/143‑1 (MQ=255)
                                                                                    agacgcttGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTa                                                        <  2:932806/137‑1 (MQ=255)
                                                                                     gACTGTTGCATAACGTCTCCGTGAATCACTGGTCGATCACAACGCCTCTGCATGGTGCAGAGGGGTATTTTTCCGGCAAATTACATAACGCCTAAACCTTTGGTCTCCAGATCATCGTTAAAGACGTCAGGCGTGGTCAggca                                                       >  2:368224/1‑140 (MQ=255)
                                                                                        tGTTACAGAACGCCTCCGTGAATCACTGGTCGTCCACACCCCCTCTGCATCATGCAGAGGAGAATTTATCCGGCAAAGTACTAACCGCCGAAACATTTTTTCTCCATCTCCTTCTTAAATTTGTCACCCGTGATCAGTACCAc                                                    >  2:409191/1‑143 (MQ=255)
                                                                                            gCAGAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCg                                                <  2:387295/143‑1 (MQ=255)
                                                                                            ccATAACGTCTCCGTGAATCACTGGTCGTTAACAACCCCTCTGCCTCGTGGTGTGGGGCGTTTTTTCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTATTTACTAAAGTTGTCACGCGTGATCAGTACCACGTCg                                                <  2:288892/142‑1 (MQ=255)
                                                                                               tAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTa                                             >  2:153441/1‑143 (MQ=255)
                                                                                                   gTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACtt                                         <  2:889818/143‑1 (MQ=255)
                                                                                                    tctcCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTc                                        >  2:877212/1‑143 (MQ=255)
                                                                                                         gTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTATAGTTGTCACGCGTGATCAGTACCACGTAGGTAACTTCGGTaa                                   >  1:1142/1‑143 (MQ=255)
                                                                                                         gTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTaa                                   >  2:792929/1‑143 (MQ=255)
                                                                                                               ctcTGGTCGCTCACACCCCCTCTGCATCATGCAGAGGGGGATTTCTCCGGCAAATTACTTACCGCCTAAGCCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGAGCAGTACCACGTCGGTAACTTCGGTAAATTTTg                             <  2:152300/141‑1 (MQ=255)
                                                                                                                  tGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTCCGGTCAATCTTGGCg                          >  2:49439/1‑143 (MQ=255)
                                                                                                                    gTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCGCGCGTGATCACTACTACGTCGGTAACTTCGGTCAATTTTGGCGGt                        >  2:229597/1‑143 (MQ=255)
                                                                                                                       gTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCa                     >  2:1019457/1‑143 (MQ=255)
                                                                                                                        ttCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCaa                    >  1:1044899/1‑143 (MQ=255)
                                                                                                                         tCACTCCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAAc                   <  2:6977/143‑1 (MQ=255)
                                                                                                                              tccccTCTGCATGATGCAGAGGGGGATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCTAGTTCTTCTTTAAAGTTGTCACGCGTGTTCAGTCCCACGTCGGTAACTTCGGCAGAGTTTGGCGGTTCACGGTCtt              >  2:428812/2‑143 (MQ=255)
                                                                                                                              tccccTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTACATTTTGGCGGTTCAACGTCtt              >  2:53059/2‑143 (MQ=255)
                                                                                                                              cccccTCTGCATCATGCAGAGGCGAATTTTTCCGGCAAATTACTTACCGCCTAAACCCTTTTCCTCCAGTTCTTCTTTAAAGTTGTCCCGCCTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCtt              <  1:552508/143‑1 (MQ=255)
                                                                                                                               ccccTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCttt             <  1:507029/143‑1 (MQ=255)
                                                                                                                                 ccTCTGCATGATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCCCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTCACTTCGGTAAATTTTGGCGCTTCACCGTCcttcg           >  1:1069659/1‑141 (MQ=255)
                                                                                                                                  ctctGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTCCCACGTCGGTAACTTCGGTAAATTTTGGCGCTTCAACGTCTTTTGc          >  2:1119449/1‑143 (MQ=255)
                                                                                                                                    ctGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTGTTTCTCCAGTTCTTCTTTAAGGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTAAGGCTGTTTATCGTCTTTTGCTa        >  2:841508/1‑143 (MQ=255)
                                                                                                                                      gCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTAcc      <  1:792929/143‑1 (MQ=255)
                                                                                                                                         tcatGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTACCCAg   <  1:425401/143‑1 (MQ=255)
                                                                                                                                          catGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTACCCAGt  >  1:509957/1‑143 (MQ=255)
                                                                                                                                          |                                                                                                                                              
GGCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTACCCAGT  >  NC_000913/1984974‑1985254

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