Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,632,837 A→G intergenic (‑552/‑53) ydfJ ← / → ynfT putative transporter YdfJ/Qin prophage; protein YnfT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,632,8370AG91.7% 39.1 / ‑2.1 24intergenic (‑552/‑53)ydfJ/ynfTputative transporter YdfJ/Qin prophage; protein YnfT
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (8/14);  total (9/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAAGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCC  >  NC_000913/1632779‑1632949
                                                          |                                                                                                                
ttAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACAtctct                                 <  1:1636319/140‑1 (MQ=17)
ttAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACAtctct                                 <  1:1891900/140‑1 (MQ=17)
  aGCATTCCTGTATACCTTAAACCACAATGTTGATCAACGGACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTcc                               <  2:631619/140‑1 (MQ=11)
  aGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTcc                               <  2:2254366/140‑1 (MQ=17)
  aGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTcc                               <  1:2398735/140‑1 (MQ=17)
  aGCATTCCTGTATACCTGAAACAACAATGTTTATCTAAGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTcc                               <  1:993220/140‑1 (MQ=255)
     aTTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTa                            >  2:1729758/1‑140 (MQ=14)
     aTTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTa                            >  2:1354661/1‑140 (MQ=17)
     aTTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTa                            <  1:45923/140‑1 (MQ=18)
       tCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATc                          >  1:558867/1‑140 (MQ=14)
        ccTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAAGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTAAATGAACAACATACTGATAATAACATCTCTCCTTATCa                         >  1:1130575/1‑140 (MQ=255)
          tGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATCATATATAACACATTACTAATCCGAGGCTAAATGGAGAACATACTGATAATCCCATCTCTCCTTATCata                       >  2:1084414/1‑140 (MQ=11)
            tataCCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATAtt                     >  1:233183/1‑140 (MQ=14)
             ataCCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTc                    <  1:1467066/140‑1 (MQ=18)
                   gAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCAtt              >  2:31299/1‑140 (MQ=14)
                   gAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCAtt              >  2:351557/1‑140 (MQ=14)
                   gAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCAtt              <  1:1296620/140‑1 (MQ=18)
                     aaCGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCAtttt            <  1:351557/140‑1 (MQ=18)
                     aaCGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCAtttt            <  1:1840425/140‑1 (MQ=18)
                             tGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATg    <  1:2006640/140‑1 (MQ=17)
                             tGTTTATCTACGAACTACAAGAACACCCAAGATCAGAATAGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATg    <  2:1130575/140‑1 (MQ=11)
                               tttATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCAtcct  <  2:1577681/140‑4 (MQ=17)
                                         aaCTTTAAGAACACCCAGGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGc   >  2:164806/1‑129 (MQ=14)
                                         aaCTTTAAGAACACCCAGGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGc   <  1:164806/129‑1 (MQ=17)
                                                          |                                                                                                                
TTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAAGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCC  >  NC_000913/1632779‑1632949

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: