Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 588,045 A→G 27.6% D711D (GAT→GAC nfrA ← bacteriophage N4 receptor, outer membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb588,0450AG27.6% 40.7 / 17.7 29D711D (GAT→GACnfrAbacteriophage N4 receptor, outer membrane subunit
Reads supporting (aligned to +/- strand):  ref base A (15/6);  new base G (6/2);  total (21/8)
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

GCAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTATCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCATGTCATCCAGACGCTGGTTC  >  W3110S.gb/587978‑588111
                                                                   |                                                                  
gCAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTa                                                                    <  1:10309/68‑1 (MQ=255)
 cAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTGTc                                                                  >  1:1155489/1‑69 (MQ=255)
 cAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTATc                                                                  <  1:1728044/69‑1 (MQ=255)
             aTTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTGTCAATGTCatcaat                                                      >  1:1862278/1‑69 (MQ=255)
             aTTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTGTCAATGTCatcaat                                                      >  1:3678806/1‑69 (MQ=255)
                gAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTATCAATGTCATCAATCAc                                                   <  1:2127124/69‑1 (MQ=255)
                gAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTATCAATGTCATCAATCAc                                                   <  1:2223499/69‑1 (MQ=255)
                gAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTATCAATGTCATCAATCAc                                                   <  1:264905/69‑1 (MQ=255)
                          tgATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTATc                                                                  >  1:345981/1‑44 (MQ=255)
                           gATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTGTCAATGTCATCAATCACCAGCCGGGCg                                         <  1:1376160/68‑1 (MQ=255)
                                            aGTGGGGTTATCAGCGCCTGATTATCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGc                       >  1:1441229/1‑69 (MQ=255)
                                            aGTGGGGTTATCAGCGCCTGATTATCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGc                       >  1:1802353/1‑69 (MQ=255)
                                            aGTGGGGTTATCAGCGCCTGATTATCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGACGCAGGc                       >  1:2304106/1‑69 (MQ=255)
                                            aGTGGGGTTATCAGCGCCTGATTATCAATGTCATAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCa                      >  1:6198/1‑69 (MQ=255)
                                               ggggTTATCAGCGCCTGATTATCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCAt                     >  1:3238783/1‑68 (MQ=255)
                                                  gTTATCAGCGCCTGATTATCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCATGTCa                 >  1:1863300/1‑69 (MQ=255)
                                                  gTTATCAGCGCCTGATTATCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCATGTCa                 >  1:1865302/1‑69 (MQ=255)
                                                     aTCAGCGCCTGATTGTCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCATGTCa                 >  1:500626/1‑66 (MQ=255)
                                                     aTCAGCGCCTGATTATCAATGTCATCAATCACCCGCCGGGCGTAGTGCTGCGTCGCAGGCATGTCa                 >  1:956982/1‑66 (MQ=255)
                                                     aTCAGCGCCTGATTATCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCATGTCa                 >  1:1129664/1‑66 (MQ=255)
                                                       cAGCGCCTGATTATCAATGTCATCAATCACCAGCCGGGCGTa                                       >  1:2378550/1‑42 (MQ=255)
                                                        aGCGCCTGATTATCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCATGTCATc               >  1:3073303/1‑65 (MQ=255)
                                                        aGCGCCTGATTATCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCATGTCATc               >  1:2522450/1‑65 (MQ=255)
                                                        aGCGCCTGATTATCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCATGTCATc               >  1:3498265/1‑65 (MQ=255)
                                                            ccTGATTGTCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGc                           >  1:1762792/1‑49 (MQ=255)
                                                            ccTGATTGTCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGc                           <  1:1150093/49‑1 (MQ=255)
                                                              tGATTGTCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCATGTCATCCAGACGCTgg     >  1:1818664/1‑69 (MQ=255)
                                                                 ttATCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCATGTCATCCAGACGCTGGTTc  >  1:227277/1‑69 (MQ=255)
                                                                  tATCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGc                       <  1:2843829/47‑1 (MQ=255)
                                                                  tATCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGc                       <  1:3237341/47‑1 (MQ=255)
                                                                   |                                                                  
GCAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTATCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCATGTCATCCAGACGCTGGTTC  >  W3110S.gb/587978‑588111

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: