Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,273,861 C→T 9.0% L349L (CTG→CTA)  garP ← galactarate/glucarate/glycerate transporter GarP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,273,8610CT 9.0% 104.9 / 3.5 45L349L (CTG→CTA) garPgalactarate/glucarate/glycerate transporter GarP
Reads supporting (aligned to +/- strand):  ref base C (17/23);  new base T (2/2);  total (20/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.12e-01

CAATGGAGGCAACATTGCCAAAGACGTTAAAGACGCCGCCGCAGAGGCCAACAATCTCTTTCGGCGCGGTGTCAGAAATCACCGGCCAGCCCAGCGCACCAAATCCTTTGCCAAAGAAAGCCAGCGCCATCAGCATGACCACCAGCGTGGTGTTGTTGGTGTAGTTACATAAGATGATGGTGGAAGCCAGCAACATTCCCAGCACAATCGGTAGCTTACGTGCCAGGGTCAGGGATAAACCGCGTTTGATCAGATAATCCGAGAAGACACCTCCCAGCA  >  NC_000913/3273719‑3273997
                                                                                                                                              |                                                                                                                                        
cAATGGAGGCAACATTGCCAAAGACGTTAAAGACGCCGCCGCAGAGGCCAACAATCTCTTTCGGCGCGGTGTCAGAAATCACCGGCCAGCCCAGCGCACCAAATCCTTTGCCAAAGAAAGCCAGCGCCACCAGCATGaccacc                                                                                                                                          <  1:286246/143‑1 (MQ=255)
   tGGAGGCAACATTGCCAAAGACGTTAAAGACGCCGCCCAAACCGCCAACAATCTCTTTCGGCGCGGTGTCAGAAATCACCGGCCAGCCCAGCGCACCAAATCCTTTGCCAAAGAAAGCCAGCGCCATCAGCATGACCACCAGc                                                                                                                                       <  2:983570/143‑1 (MQ=255)
        gCAACATTGCCAAAGACGTTAAAGACGCCGCCGCAGAGGCCAACAATCTCTTTCGGCGCGGTGTCAGAAATCACCGGCCAGCCCAGCGCACCAAATCCTTTGCCAAAGAAAGCCAGCGCCATCAGCATGACCACCAGCgtggt                                                                                                                                  >  2:667596/1‑143 (MQ=255)
        gCAACATTGCCAAAGACGTTAAAGACGCCGCCGCAGAGGCCAACAATCTCTTTCGGCGCGGTGTCAGAAATCACCGGCCAGCCCAGCGCACCAAATCCTTTGCCAAAGAAAGCCAGCGCCATCAGCATGACCACCAGCgtggt                                                                                                                                  <  1:854536/143‑1 (MQ=255)
        gCAACATTGCCAAAGACGTTAAAGACGCCGCCGCAGAGGCCAACAATCTCTTTCGGCGCGGAGTCAGAAATCACCGGCCAGCCCAGCGCACCAAATCCTTTGCCAAAGAAAGCCAGCGCCATCAGCATGACCACCAGCgtggt                                                                                                                                  >  1:202393/1‑143 (MQ=255)
          aaCATTGCCAAAGACGTTAAAGACGCCGCCGCAGAGGCCAACAATCTCTTTCGGCGCGGTGTCAGAAATCACCGGCCAGCCCAGCGCACCAAATCCTTTGCCAAAGAAAGCCAGCGCCATCAGCATGACCACCAGCGTGgtgt                                                                                                                                <  1:683828/143‑1 (MQ=255)
          aaCATGGCCAAAGACGTTAAAGACGCCGCCGCAGAGGCCAACAAGCTCTTTCGGCGCGGTGTCAGAAATCACCGGCCAGCCCAGCGCACCAAATCCTTTGCCAAAGAAAGCCAGCGCCATCAGCATGACCACCAGCGTGgtgt                                                                                                                                <  1:288510/143‑1 (MQ=255)
             aTTGCCAAAGACGTTAAAGACGCCGCCGCAGAGGCCAACAATCTCTTTCGGCGCGGTGTCAGAAATCACCGGCCAGCCCAGCGCACCAAATCCTGTGCCAAAGAAAGCCAGCGCCAGCAGCATGACCACAAGCccggtgcggt                                                                                                                             >  2:670225/1‑139 (MQ=255)
              ttGCCAATGCCCTTAAAGACGCCGCCGCAGAGGCCAACAATCTCTTTCGGCGCGTTGTCAGAAATCACCGGCCAGCCAAGCGCACCAAGTCCTTTGCCAAAGAAAGCCAGCGCCATCAGCATGACCACCAGCGTGGTgttgtt                                                                                                                            <  1:575183/143‑1 (MQ=255)
               tgcgagaGACGTTGGAGACGCCGCCGCAGAGGCCAGCAATCTCTCTCGGCGCGGTGTCAGAAATCACCGGCCAGCCGAGCGCACCAAAGCCTTGGCCAAGGAAAGCCAGCGCCATCAGCATGACCACCAGCGTGGTGttgttg                                                                                                                           <  2:523095/137‑1 (MQ=255)
                           cagaGACGCCGCCGCAGCGCCCAACATTCGCTTTCGGCGCGGTGTCAGAAATCACCGGCCAGCCCAGCGCACCAAATCCTTTGCCAAAGAATGCCAGCGCCATCAGCATGACCACCAGCGTGGTGTTGTTGGTGTAGTTACAt                                                                                                               <  2:929346/140‑1 (MQ=255)
                            aaaGACGCCGTCGCAGAGGCCAACAATCTCTTTCGGCGCGGTGTCAGAAATCACCGGCCAGCCCAGCGCACCAAATCCTTTGCCAAAGAAAGCCAGCGCCATCAGCATGACCACCAGCGTGGTGTTGTTGGTGTAGTTACATa                                                                                                              <  2:1261/143‑1 (MQ=255)
                              aGACGCCGCCGCAGAGGCCAACAATCTCCCTCGCCCGTGTGTCAGAAATCACCGGCCAGCCCAGCGCACCAAATCCTTTGCCAAAGAAAGCCAGCGCCATCAGCATGACCACCAGCGTGGTGTTGTTGGTGTAGTTACATAAg                                                                                                            <  1:273170/143‑1 (MQ=255)
                                aCGCCGCCGCAGATGCTAACAATCTCTTTCGGCGCGGTGTCAGACATCACCGGCCAGCCCAGCGCACCAAATCCTTTGCCGAAGAAAGCCAGCGCCATCAGCATGACCACCAGCGTGGTGTTGTTGGTGTAGTTACATAAgat                                                                                                          <  2:294733/143‑1 (MQ=255)
                                    cgccgcAGAGGCCAACAATCTCTTTCGGCGCGGTGTCAGAAATCACCGGCCAGCCCAGCGCACCAAATCCTTTGCCAAAGAAAGCCAGCGCCATCAGCATGACCACCAGCGTGGTGTTGTTGGTGTAGTTACATAAgatgacg                                                                                                      >  2:694730/1‑141 (MQ=255)
                                      ccgcAGAGGCCAACAATCTCTTTCGGCGCGGTGTCAGAAATCACCGGCCAGCCCAGCGCACCAAATCCTTTGCCAAAGAAAGCCAGCGCCATCAGCATGACCACCAGCGTGGTGTTGTTGGTGTAGTTACATAAGATGAtggt                                                                                                    <  1:883341/143‑1 (MQ=255)
                                         cgggggCCAACAACCCCTTTCGGCGCGTGGCCGGAAATCACCGGCCAGCCCAGCGCACCAATTCCTTTGCCAAAGAAAGCCGGCGCCATCAGCATGACCACCGGCGTGGTGTTGTGGGTGTAGTTACATAAGATGATGGTGGa                                                                                                 <  2:427500/139‑1 (MQ=255)
                                           gagGCCAACAATCTCCTTCGGCGCGGTGTCAGAAAGCACCGGCCAGCCCAGCGCACCAAAGCCTTGGCCAAAGAAAGCCAGCGCCATCAGCATGACCACCAGCGTGGTGTTGTTGGTGTAGTTACATAAGATGATGGTGGAAg                                                                                               <  1:174579/143‑1 (MQ=255)
                                                    aaTCTCTTTCGGCGCGGTGTCAGAAATCACCGGCCAGCCCAGCGCACCAACTCCTTTGCCAATGAAAGCCAGCGCCATCAGCATGACCACCAGCGTGGTGTTGTTGGTGTAGTTACATAAGATGATGGTGGAAGCCAGCAACa                                                                                      <  2:879554/143‑1 (MQ=255)
                                                          tttCGGCGCGGTGTCAGAAATCACCGGCCAGCCCAGCGCACCAAATCCTTTGCCAAAGAAAGCCAGCGCCATCAGCATGACCACCAGCGTGGTGTTGTTGGTGTAGTTACATAAGATGATGGTGGAAGCCAGCAACATTCCCa                                                                                >  1:581670/1‑143 (MQ=255)
                                                          tttCGGCGCGGTGTCAGAAATCACCGGCCAGACCAGCGGAACAAATCCTTTGACAAAGAAAGACAGCGACATCAGCATGACCAGCAGCGCGGAGATGATGGTGGAGTTACATAACATGATGCTGGAAGCAAGCAACATAccga                                                                                >  2:89237/1‑141 (MQ=255)
                                                           ttCGGCGCGGTGTCAGAAATCACCGGCCAGCCCAGCGCACCAAATCCTTTGCCAAAGAAAGCCAGCGCCATCAGCATGACCACCAGCGTGGTGTTGTTGGTGTAGTTACATAAGATGATGGTGGAAGCCAGCAACATTCCCAg                                                                               <  1:667596/143‑1 (MQ=255)
                                                               gcgcGGTGTCAGAAATCACCGGCCAGCCCAGCGCACCAAATCCTTTGCCAAAGAAAGCCAGCGCCATCAGCATGACCACCAGCGTGGTGTTGTTGGTGTAGTTACATAAGACGATGGTGGAAGCCAGCCACATCCCCTGccca                                                                           >  1:396140/1‑143 (MQ=255)
                                                                 gcgGTGACAGAAATCACCGGCCAGCCCAGCGCACCAAATCCTTTGCCAAAGAAAGCCAGCGCCATCAGCATGACCACCAGCGTGGTGTTGTTGGTGTAGTTACATAAGATGATGGTGGAAGCCAGCAACATTCCCAGCACAAt                                                                         <  1:380910/143‑1 (MQ=255)
                                                                         aGAAATCACCGGCCAGCCCAGCGCACCAAATCCTTTGCCAAAGAAAGCCAGCGCCATCAGCATGACCACCAGCGTGGTGTTGTTGGTGTAGTTACATAAGATGATGGTGGAAGCCAGCAACATTCCCAGCACAATCGGTAGCt                                                                 >  1:666406/1‑143 (MQ=255)
                                                                                cccGGCCAGCCCAGCGCACCAAATCCTTTGCCAAAGAAAGCCAGCGCCATCAGCATGACCACCAGCGTGGTGTTGTTGGTGTAGTTACATAAGATGATGGTGGAAGCCAGCAACATTCCCAGCACAATCGGTAGCTTACGTGc                                                          <  1:105977/142‑1 (MQ=255)
                                                                                      cagcccagcGCACCAAATCCGTTGGCAAAGAAAGCCAGGGCCATCAGCATGACCACCAGCGTGGTGTTGTTGGTGGAGTTACATAAGATGATGGTGGAAGCCAGCAACATTCCCAGCACAATCGGTAGCTTACGTGCCAGGGt                                                    <  1:184449/143‑1 (MQ=255)
                                                                                                  ccAAATCCTTTGCCAAAGAAAGCCAGCGCCATCAGCATGACCACCAGCGTGGTGTTGTTGGTGTAGTTACATAAGATGATGGTGGAAGCCAGCAACATGCCCAGCACAATCGGTAGCTTACGTGCCAGGGTCAGGGATAAAcc                                        <  2:849053/143‑1 (MQ=255)
                                                                                                  ccAAATCCTTTGCCAAAGAAAGCCAGCGCCATCAGCATGACCACCAGCGTGGTGTTGTTGGTGTAGTTACATAAGATGATGGTGGAAGCCAGCAACATGCCCAGCACAATCGGTAGCTCACGTGCCCGGGTCCGGGATAGAcc                                        >  2:329855/1‑143 (MQ=255)
                                                                                                     aaTCCTTTGCCAAAGAAAGCCAGCGCCATCAGCATGACCACCAGCGTGGTGTTGTTGGTGTAGTTACATAAGATGATGGTGGAAGCCAGCACCATGCCCAGCACAATCGGTAGCGTACGTGCCCGGGCCAGCGATAAACcgcg                                     >  2:221332/1‑143 (MQ=255)
                                                                                                     aaTCCTTTGCCAAAGAAAGCCAGCGCCATCAGCATGACCACCAGCGGGGTGTTGTTGGTGTAGTTACATAAGATGATGGTGGAAGCCAGCAACATTCCCAGCACAATCGGTAGCTTACGTGCCAGGGTCAGGGATAAACcgcg                                     <  1:930888/143‑1 (MQ=255)
                                                                                                          tttGCCAAAGAAAGCCAGCGCCATCAGCATGACCACCAGCGTGGTGTTGTTGGTGTAGTTACATAAGCTGATGGTGGAAGCCAGCAACATTCCCAGCACAATCGGTAGCTTACGTGCCAGGGTCAGGGATAAACCGCGTTTGa                                >  1:130914/1‑143 (MQ=255)
                                                                                                           tggCCAAAGGAAGCCAGCGCCAGCAGCATGACCGCCAGCGTGGTGTTGTTGGTGTAGCTACATAAGATGATGGTGGAAGCCAGCAACATTCCCAGCACAATCGGTAGCTTACGTGCCAGGGTCAGGGATAAACCGCGTTTGAt                               <  2:781491/141‑1 (MQ=255)
                                                                                                             gCCAAAGAAAGCCAGCGCCATCAGCATGACCACCAGCGTGGTGTTGTTGGTGTAGTTACATCAGATGCTGGTGGAAGCCAGCAACATTCCCAGCACAATCGGTAGCTCACGTGCACGGGTCAGGGGTAATCGGCGTTTGagcc                             >  2:118815/1‑140 (MQ=255)
                                                                                                              ccAAAGAAAGCCAGCGCCATCAGCATGACCACCAGCGTGGTGTTGTTGGTGTAGTTACATAAGATGATGGTGGAAGCCAGCAACATTCCCAGCACAATCGGTAGCTTACGTGCCAGGGGCAGGGATAAACCGCGTTTGGTCAg                            >  1:33982/1‑143 (MQ=255)
                                                                                                               cAAAGAAGGGCAGCGCTATCAGCATGACCACCAGCGTGGTGTTGTTAGTGTGGTTACGTAAGATGATGGTGGAAGCCAGCAACATTCCCAGCACAATCGGTAGCTTACGTGCCAGGGTCAGGGATAAACCGCGTTTGATCAGa                           <  2:472645/143‑1 (MQ=255)
                                                                                                                aaagaaagCCAGCGCCATCAGCATGACCACCAGCGTGGTGTTGTTGGTGTAGTTACATAAGATGATGGTGGAAGCCAGCAACATTCCCAGCACAATCGGTAGCTTACGTGCCAGGGTCAGGGATAAACGGCGTTGGATCAGag                          >  2:266077/1‑142 (MQ=255)
                                                                                                                aaagaaagCCAGCGCCATCAGCATGACCACCAGCGTGGTGTTGTTGGTGTAGTTACATAAGATGATGGTGGAAGCCAGCAACATTCCCAGCACAATCGATAGCTTACGTGCCAGGGTCAGGCATAAACCGCGTTTGATCAGAt                          >  2:572719/1‑143 (MQ=255)
                                                                                                                    aaagCCAGCGCCATCAGCATGACCACTAGCGTGGTGTTGTTGGTGTAGTTACATa                                                                                                              >  2:719014/1‑55 (MQ=255)
                                                                                                                    aaagCCAGCGCCATCAGCATGACCACTAGCGTGGTGTTGTTGGTGTAGTTACATa                                                                                                              <  1:979314/55‑1 (MQ=255)
                                                                                                                    aaagCCAGCGCCATCAGCATGACCACTAGCGTGGTGTTGTTGGTGTAGTTACATa                                                                                                              <  1:719014/55‑1 (MQ=255)
                                                                                                                    aaagCCAGCGCCATCAGCATGACCACTAGCGTGGTGTTGTTGGTGTAGTTACATa                                                                                                              >  2:979314/1‑55 (MQ=255)
                                                                                                                        ccAGCGCCATCAGCATGACCACCAGCGTGGTGTTGTTGGTGTAGTTACATAAGATGATGGTGGAAGCCAGCAACATTCCCAGCACAATCGGTAGCTTACGTGCCAGGGTCAGGGATAAACCGCGTTTGATCAGATAATCCgag                  >  1:285442/1‑143 (MQ=255)
                                                                                                                         cAGCGCCATCAGCATGACCACCAGCGTGGTGTTGTTGGTGTAGTTACATAAGATGATGGTGGAAGCCAGCAACATTCTCAGCACAATCGGTAGCTTAAGTGCCAGGGTCAGGGATAAACCGCGTTTGATCAGATAATCCgaga                 >  2:83225/1‑143 (MQ=255)
                                                                                                                           gcgcCAGCAGCATGACCACCAGCGGGGGGTTGTTGGGGTAGTTACATAAGATGATGGTGGAAGCCAGCAACATTCCCAGCACAATCGGTAGCTTACGTGCCAGGGTCAGGGATAAACCGCGTTTGATCAGATAATCCGagaag               <  1:699730/143‑1 (MQ=255)
                                                                                                                                      acggctcccaGCGTGGTGTTGTTGGTATGCTGACAGAAGCTGATGGTGGAATCTAGCAACATTCCCAGCACACTCGGTAGCTTACGTGCCAGGGTCAGGGATAAACCGCGTTTGATCAGATAATCCGAGAAGACACCTCCCAg    <  2:9184/136‑1 (MQ=255)
                                                                                                                                        gACCACCAGCGTGGTGTTGTTGGTGTAGTTACATAAGATCACGGTGGAAGCCAGCAACATTCCCAGCACAATCGGTAGCTTACGTGCCAGGCTCAGCGCTAACCCGCGTCCGATGAGCTAGTCCGAGAGGACACCACCCCGcc  >  2:294103/1‑142 (MQ=255)
                                                                                                                                              |                                                                                                                                        
CAATGGAGGCAACATTGCCAAAGACGTTAAAGACGCCGCCGCAGAGGCCAACAATCTCTTTCGGCGCGGTGTCAGAAATCACCGGCCAGCCCAGCGCACCAAATCCTTTGCCAAAGAAAGCCAGCGCCATCAGCATGACCACCAGCGTGGTGTTGTTGGTGTAGTTACATAAGATGATGGTGGAAGCCAGCAACATTCCCAGCACAATCGGTAGCTTACGTGCCAGGGTCAGGGATAAACCGCGTTTGATCAGATAATCCGAGAAGACACCTCCCAGCA  >  NC_000913/3273719‑3273997

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: ‑