Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 588,114 G→A 14.8% T947T (ACC→ACT nfrA ← bacteriophage N4 receptor, outer membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913588,1140GA14.8% 70.4 / 5.2 27T947T (ACC→ACTnfrAbacteriophage N4 receptor, outer membrane subunit
Reads supporting (aligned to +/- strand):  ref base G (19/4);  new base A (2/2);  total (21/6)
Fisher's exact test for biased strand distribution p-value = 2.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

AATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGCTCTCCGTTACGTTGATTAATCGCCTTAAAGGTATGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGACTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGCTACCTTCTGATGCCAGCTGACGCGATAATCTGC  >  NC_000913/587967‑588260
                                                                                                                                                   |                                                                                                                                                  
aaTGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGCTCTCCGTTACGTTGATTAATCGCCTTAAAGGTATGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGt                                                                                                                                                  >  2:229126/1‑150 (MQ=255)
  tGAACTTACGCATTTACCAGTGCACTCCAAGGGTGAGAAACGCGTTGTTGCGCTCTCCGTTACGTTGATTAATCGCCTTAAAGGTATGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTcc                                                                                                                                                <  1:211789/150‑1 (MQ=255)
                           tCCAATGGTGAGAAACGCGTTGTTGCGCTCTCCGTTACGTTGATTAATCGCCTTAAAGGTATGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGACTcc                                                                                                                       >  1:216121/1‑150 (MQ=255)
                           tCCAATGGTGAGAAACGCGTTGTTGCGCTCTCCGTTACGTTGATTAATCGCCTTAAAGGTATGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGACTcc                                                                                                                       >  1:41254/1‑150 (MQ=255)
                              aaTGGTGAGAAACGCGTTGTTGCGCTCTCCGTTACGTTGATTAATCGCCTTAAAGGTATGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGACTCCGcc                                                                                                                    >  2:198095/1‑150 (MQ=255)
                                   tGAGAAACGCGTTGTTGCGCTCTCCGTTACGTTGATTAATCGCCTTAAAGGTATGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGACTCCGCCAAGCt                                                                                                               >  2:65635/1‑150 (MQ=255)
                                                 ttgCGCTCTCCGTTACGTTGATTAATCGCCTTAAAGGTATGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCAGTCCAGATATTCCAGCGGACCCCGACTCCGCCAAGCTGCGCGCCCTGAGTg                                                                                                 >  2:24145/1‑150 (MQ=255)
                                                 ttgCGCTCTCCGTTACGTTGATTAATCGCCTTAAAGGTATGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCAGTCCAGATATTCCAGCGGACCCCGACTCCGCCAAGCTGCGCGCCCTGAGTg                                                                                                 >  2:24141/1‑150 (MQ=255)
                                                                 gTTGATTAATCGCCTTAAAGGTATGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGACTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGc                                                                                 >  1:176935/1‑150 (MQ=255)
                                                                      ttAATCGCCTTAAAGGTATGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGACTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTg                                                                            >  2:1878/1‑150 (MQ=255)
                                                                            gCCTTAAAGGTATGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGACTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGa                                                                      >  2:86460/1‑150 (MQ=255)
                                                                            gCCTTAAAGGTATGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGACTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGa                                                                      >  2:151403/1‑150 (MQ=255)
                                                                            gCCTTAAAGGTATGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGACTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGa                                                                      >  2:151402/1‑150 (MQ=255)
                                                                             ccTTAAAGGTATGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGACTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGaa                                                                     >  1:45659/1‑150 (MQ=255)
                                                                                       aTGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGACTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTa                                                           <  2:160225/150‑1 (MQ=255)
                                                                                          ttGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGACTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGACCGTGAGCGTAAgg                                                        >  1:54828/1‑150 (MQ=255)
                                                                                                                       gCGGCCAGGCGTCGTAGTGCGTCTCGCCAGTCCAGATATTCCAGCGGACCCCGACTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGCTACCTTCt                           <  1:24145/150‑1 (MQ=255)
                                                                                                                       gCGGCCAGGCGTCGTAGTGCGTCTCGCCAGTCCAGATATTCCAGCGGACCCCGACTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGCTACCTTCt                           <  1:24141/150‑1 (MQ=255)
                                                                                                                         ggCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGACTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGCTACCTTCTGa                         <  2:176935/150‑1 (MQ=255)
                                                                                                                          gCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGACTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGCTACCTTCTGAt                        >  1:227918/1‑150 (MQ=255)
                                                                                                                          gCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGACTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGCTACCTTCTGAt                        >  1:194912/1‑150 (MQ=255)
                                                                                                                          gCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGACTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGCTACCTTCTGAt                        >  1:193712/1‑150 (MQ=255)
                                                                                                                          gCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGACTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGCTACCTTCTGAt                        >  1:193688/1‑150 (MQ=255)
                                                                                                                                        tGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGACTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGCTACCTTCTGATGCCAGCTGACGCGa          >  1:129067/1‑150 (MQ=255)
                                                                                                                                            tctcGCCGGTCCAGATATTCCAGCGGACCCCGACTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGCTACCTTCTGATGCCAGCTGACGCGATAAt      >  1:45337/1‑150 (MQ=255)
                                                                                                                                              tcGCCGGTCCAGATATTCCAGCGGACCCCGACTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGCTACCTTCTGATGCCAGCTGACGCGATAATCt    >  2:172316/1‑150 (MQ=255)
                                                                                                                                                gCCGGTCCAGATATTCCAGCGGACCCCGACTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGCTACCTTCTGATGCCAGCTGACGCGATAATCTGc  <  1:65635/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGCTCTCCGTTACGTTGATTAATCGCCTTAAAGGTATGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGACTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGCTACCTTCTGATGCCAGCTGACGCGATAATCTGC  >  NC_000913/587967‑588260

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