Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,405,345 G→A 27.6% intergenic (+67/‑9) hfq → / → hflX HF‑I, host factor for RNA phage Q beta replication/predicted GTPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,405,3450GA27.6% 43.4 / 17.4 29intergenic (+67/‑9)hfq/hflXHF‑I, host factor for RNA phage Q beta replication/predicted GTPase
Reads supporting (aligned to +/- strand):  ref base G (13/8);  new base A (5/3);  total (18/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

TAAGGTTTCGGGCTGTTTTTTTACACGGGGAGCCAGCGATCCTGCGTTCCCCGCTGATCTATTTAGAGGGTTATACGCTTGTTTGACCGTTATGATGCTGGTGAGCAGGCGGTACTGGTACACATCTATTTTAC  >  W3110S.gb/4405276‑4405409
                                                                     |                                                                
tAAGGTTTCGGGCTGTTTTTTTACACGGGGAGCCAGCGATCCTGCGTTCCCCGCTGATCTATTTAGAGGGt                                                                 >  1:1727305/1‑71 (MQ=255)
  aGGTTTCGGGCTGTTTTTTTACACGGGGAGCCAGCGATCCTGCGTTCCCCGCTGATCTATTTAGAGGGTta                                                               <  1:561954/71‑1 (MQ=255)
      ttCGGGCTGTTTTTTTACACGGGGAGCCAGCGATCCTGCGTTCCCCGCTGATCTATTTAGAGGGTTATACg                                                           >  1:3784107/1‑71 (MQ=255)
                 tttttACACGGGGAGCCAGCGATCCTGCGTTCCCCGCTGATCTATTTAGAGGGTTATACGCTTGTTTGa                                                  >  1:2629187/1‑69 (MQ=255)
                 tttttACACGGGGAGCCAGCGATCCTGCGTTCCCCGCTGATCTATTTAGAGGGTTATACGCTTGTTTGAcc                                                >  1:1670664/1‑71 (MQ=255)
                                ccAGCGATCCTGCGTTCCCCGCTGATCTATTTAGAGGGTTATACGCTTGTTTGACCGTTATGATGCTGGt                                  >  1:1549671/1‑70 (MQ=255)
                                ccAGCGATCCTGCGTTCCCCGCTGATCTATTTAGAGGGTTATACGCTTGTTTGACCGTTATGATGCTGGt                                  >  1:2056693/1‑70 (MQ=255)
                                   gCGATCCTGCGTTACCCGCTGATCTATTTAGAGGATTATACg                                                           >  1:1762902/1‑42 (MQ=37)
                                   ccGATCCTGCGTTCCCCGCTGATCTATTTAGAGGGTta                                                               <  1:2654785/37‑1 (MQ=255)
                                    cGATCCTGCGTTCCCCGCTGATCTATTTAGAGGGTTATACGCTTGTTTGACCGTTATGATGCTGGTGAGCa                             <  1:3277295/71‑1 (MQ=255)
                                     gATCCTGCGTTCCCCGCTGATCTATTTAGAGGATTATACGCTTGTTTGACCGTTATGATGCTGGt                                  >  1:3833096/1‑65 (MQ=255)
                                     gATCCTGCGTTCCCCGCTGATCTATTTAGAGGATTATACGCTTGTTTGACCGTTATGATGCTGGt                                  >  1:2068188/1‑65 (MQ=255)
                                      aTCCTGCGTTCCCCGCTGATCTATTTAGAGGGTTATACGCTTGTTTGACCGt                                              <  1:852725/52‑1 (MQ=255)
                                        ccTGCGTTCCCCGCTGATCTATTTAGAGGGTTATACGCTTGTTTGACCGTTat                                           >  1:3540830/1‑53 (MQ=255)
                                           gCGTTCCCCGCTGATCTATTTAGAGGGTTATACGCTTGTTTGACCGTTatga                                         >  1:347618/1‑52 (MQ=255)
                                           gCGTTCCCCGCTGATCTATTTAGAGGGTTATACGCTTGTTTGACCGTTatga                                         >  1:1083124/1‑52 (MQ=255)
                                            cGTTCCCCGCTGATCTATTTAGAGGGTTATACGCTTGTTTGACCGTTATGATGCTGGTGAGCAGGCGGTAc                     <  1:15278/71‑1 (MQ=255)
                                              ttCCCCGCTGATCTATTTAGAGGATTATACGCTTGTTTGACCGTTATGATGCTGGTGAGCAGGCGGTACt                    <  1:48067/70‑1 (MQ=255)
                                                ccccGCTGATCTATTTAGAGGATTATACGCTTGTTTGACCGTTATGATGCTGGTGAGCAGGCg                         >  1:4064915/1‑63 (MQ=255)
                                                ccccGCTGATCTATTTAGAGGATTATACGCTTGTTTGACCGTTATGATGCTGGTGAGCAGGCg                         >  1:4069565/1‑63 (MQ=255)
                                                 cccGCTGATCTATTTAGAGGGTTATACGCTTGTTTGACCGTTATGATGCTGGTGAGCAGGCGGTACTGGTa                <  1:2030964/71‑1 (MQ=255)
                                                 cccGCTGATCTATTTAGAGGGTTATACGCTTGTTTGACCGTTATGATGCTGGTGAGCAGGCGGTACTGGTa                <  1:2756057/71‑1 (MQ=255)
                                                    gCTGATCTATTTAGAGGGTTATACGCTTGTTTGACCGt                                              <  1:1723407/38‑1 (MQ=255)
                                                      tGATCTATTTAGAGGATTATACGCTTGTTTGACCGTTATGATGCTGGTGAg                               <  1:1893449/51‑1 (MQ=255)
                                                        aTCTATTTAGAGGATTATACGCTTGTTTGACCGTTATGATGCTGGTGAGCAGGCGGTACTgg                  <  1:3271305/62‑1 (MQ=255)
                                                         tCTATTTAGAGGGTTATACGCTTGTTTGACCGTTATGATGCTGGTGAGCAGGCGGTACTGGTACACATCTa        >  1:2455391/1‑71 (MQ=255)
                                                          cTATTTAGAGGGTTATACGCTTGTTTGACCGTTATGATGCTGGTGAGCAGGCGGTACTGGTACACATCTAt       >  1:2092452/1‑71 (MQ=255)
                                                          cTATTTAGAGGGTTATACGCTTGTTTGACCGTTATGATGCTGGTGAGCAGGCGGTACTGGTACACATCTAt       >  1:1880485/1‑71 (MQ=255)
                                                               tAGAGGGTTATACGCTTGTTTGACCGTTATGATGCTGGTGAGCAGGCGGTACTGGTACACATCTATTTTAc  >  1:2714213/1‑71 (MQ=255)
                                                                     |                                                                
TAAGGTTTCGGGCTGTTTTTTTACACGGGGAGCCAGCGATCCTGCGTTCCCCGCTGATCTATTTAGAGGGTTATACGCTTGTTTGACCGTTATGATGCTGGTGAGCAGGCGGTACTGGTACACATCTATTTTAC  >  W3110S.gb/4405276‑4405409

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