Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,068,285 Δ1 bp intergenic (+31/+226) ymdF → / ← rutG KGG family protein/pyrimidine permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,068,2850T.97.9% 191.5 / ‑2.1 48intergenic (+31/+226)ymdF/rutGKGG family protein/pyrimidine permease
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base . (21/26);  total (22/26)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

GACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCA  >  NC_000913/1068157‑1068394
                                                                                                                                |                                                                                                             
gACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTc                                                                                                   >  1:1237248/1‑139 (MQ=255)
gACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTc                                                                                                   >  1:955361/1‑139 (MQ=255)
          cGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCAt                                                                                         >  1:872924/1‑139 (MQ=255)
             gggTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCAc                                                                                      <  2:1237248/139‑1 (MQ=255)
                            tGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTCTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCg                                                                         >  2:1471591/1‑139 (MQ=255)
                            tGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCggt                                                                       >  1:571047/1‑139 (MQ=255)
                              aCCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCaaa                                                                                                <  2:248071/112‑1 (MQ=255)
                              aCCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCaaa                                                                                                >  1:248071/1‑112 (MQ=255)
                                    aGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATa                                                               <  1:99705/139‑1 (MQ=255)
                                           aTCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACtgatga                                                        <  2:571047/139‑1 (MQ=255)
                                            tCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGAc                                                       <  2:457691/139‑1 (MQ=255)
                                            tCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGAc                                                       <  1:313008/139‑1 (MQ=255)
                                            tCTGAAGCAGGCAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGAc                                                       <  2:955361/139‑1 (MQ=255)
                                                     ggTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTAcc                                                                           >  2:799102/1‑110 (MQ=255)
                                                     ggTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTAcc                                                                           <  1:799102/110‑1 (MQ=255)
                                                                  gtAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCg                                                                         >  2:653820/1‑99 (MQ=255)
                                                                  gtAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCg                                                                         <  1:653820/99‑1 (MQ=255)
                                                                  gtAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATgcg                                 <  2:1008221/139‑1 (MQ=255)
                                                                    aaGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGcgcg                               >  2:1234503/1‑139 (MQ=255)
                                                                     agagCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTAcc                                                                           >  1:1257073/1‑94 (MQ=255)
                                                                     agagCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTAcc                                                                           <  2:1257073/94‑1 (MQ=255)
                                                                         cAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCggtggt                                                                    <  1:1433283/97‑1 (MQ=255)
                                                                         cAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCggtggt                                                                    >  2:1433283/1‑97 (MQ=255)
                                                                         cAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACa                                                      >  2:800874/1‑111 (MQ=255)
                                                                         cAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACa                                                      <  1:800874/111‑1 (MQ=255)
                                                                            tCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCgg                                                                                <  2:361492/82‑1 (MQ=255)
                                                                            tCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCgg                                                                                >  1:361492/1‑82 (MQ=255)
                                                                            tCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATg                            >  2:509787/1‑134 (MQ=255)
                                                                            tCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATg                            <  1:509787/134‑1 (MQ=255)
                                                                                ggCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAgg                   <  1:1234503/139‑1 (MQ=255)
                                                                                         gACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGAc                                                       <  2:557639/94‑1 (MQ=255)
                                                                                         gACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGAc                                                       >  1:557639/1‑94 (MQ=255)
                                                                                                 accGGGCTAATCAATGTCGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATTAGGCCTACCAGAAGATCGCa  <  2:997852/138‑1 (MQ=255)
                                                                                                 accGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCa  <  2:533930/138‑1 (MQ=255)
                                                                                                         aaTCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTAcc                                                                           >  1:826691/1‑58 (MQ=255)
                                                                                                         aaTCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTAcc                                                                           <  2:826691/58‑1 (MQ=255)
                                                                                                         aaTCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCa  >  2:513110/1‑131 (MQ=255)
                                                                                                         aaTCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCa  <  1:513110/131‑1 (MQ=255)
                                                                                                                gACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGt                                                                               <  1:719455/47‑1 (MQ=255)
                                                                                                                gACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGt                                                                               >  2:719455/1‑47 (MQ=255)
                                                                                                                     aatcATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACt                                                             <  2:1361408/58‑1 (MQ=255)
                                                                                                                     aatcATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACt                                                             >  1:1361408/4‑61 (MQ=255)
                                                                                                                     aatcATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATg                            <  2:1394423/91‑1 (MQ=255)
                                                                                                                     aatcATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATg                            >  1:1394423/4‑94 (MQ=255)
                                                                                                                     aatcATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAgg                   >  2:550146/4‑103 (MQ=255)
                                                                                                                     aatcATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAgg                   <  1:550146/100‑1 (MQ=255)
                                                                                                                      atcATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCggtggt                                                                    <  2:1170641/51‑1 (MQ=255)
                                                                                                                      atcATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCggtggt                                                                    >  1:1170641/3‑53 (MQ=255)
                                                                                                                                |                                                                                                             
GACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCA  >  NC_000913/1068157‑1068394

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