Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,897,018 C→T 9.1% A285T (GCC→ACC)  ygcS ← putative transporter YgcS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,897,0180CT 9.1% 114.3 / 3.2 44A285T (GCC→ACC) ygcSputative transporter YgcS
Reads supporting (aligned to +/- strand):  ref base C (21/19);  new base T (2/2);  total (23/21)
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.48e-01

GAAGGCAAACAGGCCATGACTACCAGCGTTGCCGCCAGCAGCAAAAAACTTCCCAGCAAAAATTTGCGATGTGCCAGCAGGTGCGTCAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCATCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGTGCGCCGCCAGTAACGCGAAGAGAACAAGGTTT  >  NC_000913/2896876‑2897158
                                                                                                                                              |                                                                                                                                            
gAAGGCAAACAGGCCATGACTACCAGCGTTGCCGCCAGCAGCAAAAAACTTCCCAGCAAAAATTTGCGATGTGCCAGCAGGTGCGTCAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCGTCAGGCTGGc                                                                                                                                              >  2:459057/1‑143 (MQ=255)
    gCAAACAGGCCATGACTACCAGCGTTGCCGCCAGCAGCAAAAAACTTCCCAGCAAAAATTTGCGATGTGCCAGCAGGTGCGGCAGAACTAATCCCAGCAGCGCGCCCACCATTAACAACGCATTAAGCATCAGGCTGGCAGTc                                                                                                                                          >  1:613249/1‑143 (MQ=255)
    acAAACAGCCCATGACTACCAGCGTTGCCGCCAGCAGCAAAAAACTTCCCAGCAAAAATTTGCGATGTGCCAGCAGGTGCGTCAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTc                                                                                                                                          <  2:841517/142‑1 (MQ=255)
       aaCAGGCCATGACTACCAGCGTTGCCGCCAGCAGCAAAAAACTTCCCAGCAAAAATTTGCGATGTGCCAGCAGGTGCGTCAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAgc                                                                                                                                       >  1:320722/1‑143 (MQ=255)
               aTGACTACCAGCGTTGCCGCCAGCAGCAAAAAACTTCCCAGCAAACATTTGCGATGTGCCAGCAGGTGCGTCAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCATCTTc                                                                                                                               <  2:1004015/143‑1 (MQ=255)
                 gACTACCAGCGTTGCCGCCAGCAGCAAAAAACTTCCCAGCAAAAATTTGCGATGTGCCAGCAGGTGCGTCAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCATCTTCCa                                                                                                                             <  2:952801/143‑1 (MQ=255)
                  tctatgaGCTTTGCCGCCAGGAGCAAAAGACTTCGCAGCAAAAGTTTGCGATGTGCCAGCAGGTGCGTCAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCATCTTCCAg                                                                                                                            <  1:322806/137‑1 (MQ=255)
                       cAGCGTTGCCGCCAGCAGCAAAAAACTTCCCAGCAAAAATTTGCGATGTGCCAGCAGGTGCGTCAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCATCTTCCAGACCaa                                                                                                                       >  2:526242/1‑143 (MQ=255)
                            ttGCCGCCAGCAGCAAAAAACTTCCCAGCAAAAATTTGCGATGTGCCAGCAGGTGCGTCAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGCCAGCACATCTTCCAGAACAATAgcc                                                                                                                  >  2:956101/1‑141 (MQ=255)
                                cgccAGCAGCAAAAAACTTCCCAGCAAAAATTTGCGATGTGCCAGCAGGTGCGTCAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCATCTTCCAGACCAATAGTCTggg                                                                                                              <  2:801834/143‑1 (MQ=255)
                                 gccAGCAGCAAAAAACTTCCCAGCAAAAATTTGCGATGTGCCAGCAGGTGCGTCAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCATCTTCCAGACCAATAGTCTGGGc                                                                                                             >  1:9924/1‑143 (MQ=255)
                                      cagcaAAAAACTTCCCAGCAAAAATTTGCGATGTGCCAGCAGGTGCGTCAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCATCTTCCAGACCAATAGCCTGGACGATAg                                                                                                        >  1:179795/1‑143 (MQ=255)
                                         caAAAAACTTCCCAGCAAACATTTGCGATGTGCCCGCAGGTGCGTCAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCATCTTCCAGACCAATAGTCTGGGCGATAGTTg                                                                                                     <  2:264386/143‑1 (MQ=255)
                                          aaaaaaCTTCCCAGCAAAAATTTGCGATGTGCCAGCAGGTGCGTCAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGTAGTCAGCGCATCTTCCAGACCAATAGTCTGGGCGATAGTTgg                                                                                                    >  2:836799/1‑143 (MQ=255)
                                          aaaaaaCTTCCCAGCAAAAATTTGCGATGTGCCAGCAGGTGCGTCAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGTAGTCAGCGCATCTTCCAGACCAATAGTCTGGGCGATAGTTgg                                                                                                    >  2:833499/1‑143 (MQ=255)
                                           aaaaaCTTCCCAGCAAAAATTTGCGATGTGCCAGCAGGTGCGTCAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCATCTTCCAGACCAATAGTCTGGGCGATAGTTGGc                                                                                                   >  1:782214/1‑143 (MQ=255)
                                            aaaaCTTCCCAGCAAAAATTTGCGATGTGCCAGCAGATGCGTCAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCATCTTCCAGCACAATAGTCTGGGCGATCGTTGGCa                                                                                                  >  1:323956/1‑143 (MQ=255)
                                                cTTCCCAGCAAAAATTTGCGATGTGCCAGCAGGTGCGTCAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCATCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCa                                                                                              <  1:776409/143‑1 (MQ=255)
                                                 ttCCCAGCAAAAATTTGCGATGTGCCAGCAGGTGCGTCAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCATCTTCCAGACCAATAGTCTGGGCGATAGGTGGCAGCCAg                                                                                             >  2:703129/1‑143 (MQ=255)
                                                    ccAGCAAAAATTTGCGATGTGCCAGCAGGTGCGTCAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCATCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGta                                                                                          >  1:883845/1‑143 (MQ=255)
                                                         aaaaaTTTGCGATGTGCCAGCAGGTGCGTCAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCATCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAAt                                                                                     >  2:42094/1‑143 (MQ=255)
                                                               ttGCGACGTGCCAGCACGTGCGTCAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCATCTTCCAGACCAATAGGCTGGTCGATAGTTGGCAGCCAGGTATAAATCACaaa                                                                               <  1:402961/143‑1 (MQ=255)
                                                                                gTGCGTCAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCATCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGCTTACgccgc                                                              >  1:995481/1‑139 (MQ=255)
                                                                                 tGCGTCAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGTAGTCAGCGCATCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCa                                                             <  1:833499/143‑1 (MQ=255)
                                                                                 tGCGTCAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGTAGTCAGCGCATCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCa                                                             <  1:836799/143‑1 (MQ=255)
                                                                                   cGTCAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCATCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACACACCATGGGATTACGAGGCAGa                                                           >  1:432310/1‑143 (MQ=255)
                                                                                    gTCAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCATCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGAc                                                          <  1:526242/143‑1 (MQ=255)
                                                                                     tCAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCATCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGCATACATCACAAACCATGGGATTACGAGGCAGACa                                                         >  2:561254/1‑143 (MQ=255)
                                                                                        gAACTAATCCCAGCCGCGCGCCCACAATGAACAACGCATTAAGCATCGGGCTGGCAGTCAGCGCATCTTCCAGACCAATAGTCTGGGCGCTAGTATGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAaag                                                      <  2:453525/143‑1 (MQ=255)
                                                                                          aCTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCATCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAaagaa                                                    <  1:803988/143‑1 (MQ=255)
                                                                                              aTCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCATCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAagaagaag                                                >  2:325290/1‑141 (MQ=255)
                                                                                               tCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCATGCTGGCAGTCAGCGCATCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGc                                               >  1:671230/1‑143 (MQ=255)
                                                                                                      agcGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCATCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTaaa                                        <  1:683370/143‑1 (MQ=255)
                                                                                                       gcgcgcCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCATCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGGCAAAGAAGACGCTGTGAAAc                                       >  1:59766/1‑143 (MQ=255)
                                                                                                            cccACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCATCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGATGTTAAACGCCGt                                  >  2:777896/1‑143 (MQ=255)
                                                                                                               acaATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCATCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGTgcg                               <  1:785708/143‑1 (MQ=255)
                                                                                                                   ttAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCATCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGTGcgccgc                           <  2:782214/143‑1 (MQ=255)
                                                                                                                   ttAACAACGCATTAAGCAGGAGGCTGGCAGTCAGCGCATCTTCCAGACAAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGTGcgccgc                           <  2:792091/143‑1 (MQ=255)
                                                                                                                    tAACAACTCATTAAGCATCAGGCTGGCAGTCAGCGCATCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGTGCgccgcc                          <  1:680815/143‑1 (MQ=255)
                                                                                                                     cacaacGCATTGAGCATCAGGCTGGCAGTGAGCGCATCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAGATCTCCAACCATGGGATTACGAGGCAGACAAAGGAGACGCTTTTAAACGCCGTGCGCCGCCa                         <  2:221103/142‑1 (MQ=255)
                                                                                                                     aacaacGCATTAAGCATCAGGCTGGCAGTCAGCGCATCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGTGCGCCGCCa                         <  2:320722/143‑1 (MQ=255)
                                                                                                                       caacGCATTAAGCATCAGGCTGGCAGTCAGCGCATCTTCCAGACCAATAGTCTGTGCGATAGTTGGCAGCCAGGGATAAATCACCAACCATGGGATTACGAGGCAGACAAAGAAGACGCGGCTACACGCCGTGCGCCGCCact                       >  1:381678/1‑141 (MQ=255)
                                                                                                                        aacGCATTAAGCATCAGGCTGGCAGTCAGCGCATCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGTGCGCCGCCAGTa                      >  1:443114/1‑143 (MQ=255)
                                                                                                                                       aGGCTGGCAGTCAGCGCATCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGTGCGCCGCCAGTAACGCGAAGAGAACaa       <  1:820209/143‑1 (MQ=255)
                                                                                                                                            ggCAGTCAGCGCATCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGTGCGACGCCAGTACCGCGAGGAGAACAAGGttt  >  1:474586/1‑143 (MQ=255)
                                                                                                                                              |                                                                                                                                            
GAAGGCAAACAGGCCATGACTACCAGCGTTGCCGCCAGCAGCAAAAAACTTCCCAGCAAAAATTTGCGATGTGCCAGCAGGTGCGTCAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCATCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGTGCGCCGCCAGTAACGCGAAGAGAACAAGGTTT  >  NC_000913/2896876‑2897158

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