Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,985,112 C→G 17.9% intergenic (‑43/+27) araG ← / ← araF L‑arabinose ABC transporter ATPase/L‑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,1120CG17.9% 152.2 / 27.2 56intergenic (‑43/+27)araG/araFL‑arabinose ABC transporter ATPase/L‑arabinose ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (22/24);  new base G (5/5);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.56e-01

TGACCGGCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTACCC  >  NC_000913/1984969‑1985251
                                                                                                                                               |                                                                                                                                            
tGACCGGCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGc                                                                                                                                          >  1:48465/1‑148 (MQ=255)
   ccGGCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCaga                                                                                                                                       <  4:561716/148‑1 (MQ=255)
   ccGGCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCaga                                                                                                                                       <  3:1008346/148‑1 (MQ=255)
               tCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTc                                                                                                                          <  2:48465/148‑1 (MQ=255)
                            ccGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCaga                                                                                                                                       >  4:888713/1‑123 (MQ=255)
                            ccGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCaga                                                                                                                                       <  3:888713/123‑1 (MQ=255)
                                       ttAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCag                                                                                                                                        >  4:599730/1‑111 (MQ=255)
                                       ttAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCag                                                                                                                                        <  3:599730/111‑1 (MQ=255)
                                         aaCGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATGATGCAGAgg                                                                                                                                    <  3:1098052/112‑1 (MQ=255)
                                         aaCGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATGATGCAGAgg                                                                                                                                    >  4:1098052/1‑112 (MQ=255)
                                                  aaaCGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATGATGCAGAgg                                                                                                                                    >  3:1073024/1‑103 (MQ=255)
                                                  aaaCGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATGATGCAGAgg                                                                                                                                    <  4:1073024/103‑1 (MQ=255)
                                                                           gagagaTACGGGGTAGACTGTTGCATAACGTCTCCGTGTATCACTGGTCGTTCACACCCCCTCTGCATCATGc                                                                                                                                          <  1:858471/73‑1 (MQ=255)
                                                                           gagagaTACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGc                                                                                                                                          >  2:858471/1‑73 (MQ=255)
                                                                                 tACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTaa                                                                          >  2:29202/1‑130 (MQ=255)
                                                                                 tACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTaa                                                                          <  1:29202/130‑1 (MQ=255)
                                                                                        tAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATGATGCAGAggg                                                                                                                                   >  3:642089/1‑66 (MQ=255)
                                                                                        tAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATGATGCAGAggg                                                                                                                                   <  4:642089/66‑1 (MQ=255)
                                                                                                 gCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCAcg                                                                >  3:848748/1‑124 (MQ=255)
                                                                                                 gCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCAcg                                                                <  4:848748/124‑1 (MQ=255)
                                                                                                      aCGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCgg                                   >  1:331842/1‑148 (MQ=255)
                                                                                                          ctcCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTc                                                                                                                          <  4:1042003/57‑1 (MQ=255)
                                                                                                          ctcCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTc                                                                                                                          >  3:1042003/1‑57 (MQ=255)
                                                                                                               tGAATCACTGGTCGTTCACACCCCCTCTGCATCTTGCAGAGGGCAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACTTCGGTAACTTCGGTaaa                               <  4:83825/143‑1 (MQ=255)
                                                                                                               tGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTaaa                               <  2:398465/143‑1 (MQ=255)
                                                                                                               tGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTaaa                               >  1:398465/1‑143 (MQ=255)
                                                                                                               tGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTaaa                               >  3:83825/1‑143 (MQ=255)
                                                                                                                gAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGttcttc                                                                               <  4:397169/94‑1 (MQ=255)
                                                                                                                gAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGttcttc                                                                               >  3:397169/1‑94 (MQ=255)
                                                                                                                      cTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGGTc                   >  3:409423/1‑148 (MQ=255)
                                                                                                                       tGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGttcttc                                                                               >  3:338005/1‑87 (MQ=255)
                                                                                                                       tGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGttcttc                                                                               <  4:338005/87‑1 (MQ=255)
                                                                                                                        ggTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGa                                                           >  4:771479/1‑106 (MQ=255)
                                                                                                                        ggTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGa                                                           <  3:771479/106‑1 (MQ=255)
                                                                                                                        ggTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCaa                 <  2:331842/148‑1 (MQ=255)
                                                                                                                                  ttccccTCTGCATGATGCAGAGGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCttttt                                                                                            >  4:210024/3‑64 (MQ=25)
                                                                                                                                  ttccccTCTGCATGATGCAGAGGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCttttt                                                                                            <  3:210024/62‑1 (MQ=25)
                                                                                                                                   cccccTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGttctt                                                                                >  1:431044/1‑74 (MQ=255)
                                                                                                                                   cccccTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGttctt                                                                                <  2:840882/74‑1 (MQ=255)
                                                                                                                                   cccccTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGttctt                                                                                <  2:431044/74‑1 (MQ=255)
                                                                                                                                   cccccTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGttctt                                                                                >  1:840882/1‑74 (MQ=255)
                                                                                                                                      ccTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGt                                                      <  1:415596/97‑1 (MQ=255)
                                                                                                                                      ccTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGt                                                      >  2:415596/1‑97 (MQ=255)
                                                                                                                                       ctctGCATGATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGttct                                                                                 >  3:788513/1‑69 (MQ=255)
                                                                                                                                       ctctGCATGATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGttct                                                                                 <  4:788513/69‑1 (MQ=255)
                                                                                                                                       ctctGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTTCCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTAccc  <  4:409423/148‑1 (MQ=255)
                                                                                                                                        tctGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCgg                                   >  1:200810/1‑114 (MQ=255)
                                                                                                                                        tctGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCgg                                   <  2:200810/114‑1 (MQ=255)
                                                                                                                                        tctGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTaaa                               >  1:62646/1‑118 (MQ=255)
                                                                                                                                        tctGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTaaa                               <  2:62646/118‑1 (MQ=255)
                                                                                                                                        tctGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCttt          <  2:854865/139‑1 (MQ=255)
                                                                                                                                        tctGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCttt          >  1:854865/1‑139 (MQ=255)
                                                                                                                                         ctGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGt                     <  1:1005255/127‑1 (MQ=255)
                                                                                                                                         ctGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGt                     >  2:1005255/1‑127 (MQ=255)
                                                                                                                                           gCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCg                                                                                                         <  2:732214/41‑1 (MQ=255)
                                                                                                                                           gCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCg                                                                                                         >  1:732214/1‑41 (MQ=255)
                                                                                                                                              tcatGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTg                                                            <  4:145649/83‑1 (MQ=255)
                                                                                                                                              tcatGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTg                                                            >  3:145649/1‑83 (MQ=255)
                                                                                                                                               |                                                                                                                                            
TGACCGGCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTACCC  >  NC_000913/1984969‑1985251

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