Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,068,157 T→C intergenic (‑77/+16) yggE ← / ← argO uncharacterized protein YggE/L‑arginine exporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,068,1570TC98.0% 165.9 / ‑4.2 49intergenic (‑77/+16)yggE/argOuncharacterized protein YggE/L‑arginine exporter
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (26/22);  total (27/22)
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.96e-01

TGCGCTGCCATCCCGCTAATACCCATTAATGCCGCCAGGGCGATAACTTTGAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTTTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAACACATCCCACAACCAGATTGATAATGCGCTGTGCTTTTGCCGTGCGCAGACG  >  NC_000913/3068022‑3068292
                                                                                                                                       |                                                                                                                                       
tGCGCTGCCATCCCGCTAATACCCATTAATGCCGCCAGGGCGATAACTTTGAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCa                                                                                                                                     >  2:717705/1‑140 (MQ=255)
  cgcTGCCATCCCGCTAATACCCATTAATGCCGCCAGGGCGATAACTTTGAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATc                                                                                                                                   >  2:529698/1‑140 (MQ=255)
            ccGCTAATACCCATTAATGCCGCCAGGGCGATAACTTTGAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGcc                                                                                                                         <  2:616860/140‑1 (MQ=255)
               cTAATACCCATTAATGCCGCCAGGGCGATAACTTTGAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTaa                                                                                                                      <  1:54324/140‑1 (MQ=255)
                    aCCCATTAATGCCGCCAGGGCGATAACTTTGAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGaa                                                                                                                 >  2:105895/1‑140 (MQ=255)
                      ccATTAATGCCGCCAGGGCGATAACTTTGAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACa                                                                                                               >  2:153871/1‑140 (MQ=255)
                         ttAATGCCGCCAGGGCGATAACTTTGAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAgg                                                                                                            >  1:94623/1‑140 (MQ=255)
                          tAATGCCGCCAGGGCTATAACTTTGGACTGCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGc                                                                                                           <  1:399096/140‑1 (MQ=255)
                          tAATGCCGCCAGGGCGATAACTTTGAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTGTTATTTCCCCATCAAATAGCGACTAACTGAACAAGGc                                                                                                           >  1:272350/1‑140 (MQ=255)
                            aTGCCGCCAGGGCGATAACTTTGAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCtt                                                                                                         <  2:428405/140‑1 (MQ=255)
                             tGCCGCCAGGGCGATAACTTTGAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGACTtg                                                                                                        <  1:437321/140‑1 (MQ=255)
                               ccgccAGGGCGATAACTTTGAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGc                                                                                                     <  2:48002/140‑1 (MQ=255)
                                cgccAGGGCGATAGCTTTGAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGc                                                                                                     <  1:161679/140‑1 (MQ=255)
                                  ccAGGGCGATAACTTTGAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCAt                                                                                                   <  2:94623/140‑1 (MQ=255)
                                  ccAGGGCGATAACTTTGAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCAt                                                                                                   >  2:788841/1‑140 (MQ=255)
                                                ttGAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGACTAACTGGACAAGGGTTGTGCCTGAGCAATACCGtct                                                                                     >  1:394403/1‑140 (MQ=255)
                                                 tGAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGtctc                                                                                    >  2:148664/1‑140 (MQ=255)
                                                             tAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCa                                                                        <  1:153871/140‑1 (MQ=255)
                                                              aGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGGCCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCaa                                                                       <  1:394512/140‑1 (MQ=255)
                                                               gTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAg                                                                      <  1:148664/140‑1 (MQ=255)
                                                                   ctcCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCa                                                                        <  2:406779/134‑1 (MQ=255)
                                                                   ctcCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCa                                                                        >  1:406779/1‑134 (MQ=255)
                                                                     ccGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATa                                                                <  1:611480/140‑1 (MQ=255)
                                                                     ccGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGAATGAGCAATACGGTGTCTCGCGAGCTGCAAGGCAATa                                                                >  1:228459/1‑140 (MQ=255)
                                                                       gtgtTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATaaa                                                              <  1:529698/140‑1 (MQ=255)
                                                                          ttAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTGTCGCCAGCTACAAGGCAATAAACca                                                           >  1:236914/1‑140 (MQ=255)
                                                                            aaTTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATaaa                                                              >  1:645598/1‑135 (MQ=255)
                                                                            aaTTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATaaa                                                              <  2:645598/135‑1 (MQ=255)
                                                                            aaTTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACcaca                                                         >  1:704647/1‑140 (MQ=255)
                                                                              ttGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATATCGTGTCTCGCCAGCTGCGAGGCAATAAACCACATa                                                       >  2:254359/1‑140 (MQ=255)
                                                                              ttGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATa                                                       >  2:500810/1‑140 (MQ=255)
                                                                                      cccAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAACACATcc                                               >  1:526687/1‑140 (MQ=255)
                                                                                      cccAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAACACATcc                                               >  2:114616/1‑140 (MQ=255)
                                                                                      cccAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAACACATcc                                               >  1:763610/1‑140 (MQ=255)
                                                                                        cAAAATTCAGGACCTCTAGAGGTAAGCTTAGCGCCCCTGCTTTATTTCCCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAACACATCCca                                             <  2:228459/140‑1 (MQ=255)
                                                                                          aaaTTCAGGACCTCTGGAGGCACGCTTAGCGACTCGGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAACACATCCcaca                                           <  2:272350/140‑1 (MQ=255)
                                                                                                aGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAACACATCCca                                             >  1:25709/1‑132 (MQ=255)
                                                                                                aGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAACACATCCca                                             <  2:25709/132‑1 (MQ=255)
                                                                                                       ctGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAACACATCCCACAAc                                         >  2:67597/1‑129 (MQ=255)
                                                                                                       ctGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAACACATCCCACAAc                                         <  1:67597/129‑1 (MQ=255)
                                                                                                          gAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAACACATCCCACAACCAGATTGATAATgc                           >  1:822969/1‑140 (MQ=255)
                                                                                                             gCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAACACATCCCACAACCAGATTGATAATGCGCt                        >  2:27154/1‑140 (MQ=255)
                                                                                                               aaGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAACACATCCCACAACCAGATTGATAATgcgc                         <  2:486790/137‑1 (MQ=255)
                                                                                                               aaGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAACACATCCCACAACCAGATTGATAATgcgc                         >  1:486790/1‑137 (MQ=255)
                                                                                                               aaGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAACACATCCACAACCAGATTGATAATGCGCtgtg                     >  1:820696/1‑140 (MQ=255)
                                                                                                                     aGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAACACATCCCACAACCAGATTGATAATGCGCTGTGCtttt                >  2:725509/1‑140 (MQ=255)
                                                                                                                            ctGTTTTATTTTTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAACACATCCCACAACCAGATTGATAATGCGCTGTGCTTTTGCCGTgc         >  1:352162/1‑140 (MQ=255)
                                                                                                                                   aTTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAACACATCCCACAACCAGATTGATAATGCGCTGTGCTTTTGCCGTGCGCAGAcg  <  1:620557/140‑1 (MQ=255)
                                                                                                                                   aTTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAACACATCCCACAACCAGATTGATAATGCGCTGTGCTTTTGCCGTGCGCAGAcg  <  1:105895/140‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                       
TGCGCTGCCATCCCGCTAATACCCATTAATGCCGCCAGGGCGATAACTTTGAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTTTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAACACATCCCACAACCAGATTGATAATGCGCTGTGCTTTTGCCGTGCGCAGACG  >  NC_000913/3068022‑3068292

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