Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 592,390 A→G S263S (AGT→AGC)  nfrB ← bacteriophage N4 receptor, inner membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913592,3900AG96.8% 102.8 / NA 31S263S (AGT→AGC) nfrBbacteriophage N4 receptor, inner membrane subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (16/14);  minor base T (0/1);  total (16/15)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTGATGTCCGCGCGTGCTGTAAAAATTTACGCTGCTCGCGTTCTTTGGCTTCGTCCACCACCGGAAAACGGACAAAAATTTCCGTCATACCTTTTTCTTTCAGGCGGAAGCCAATGTCGTAATCTTCAGTAAGACTCTGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCTGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTCACGCACCGGAACATCTTTGCCATGCA  >  NC_000913/592256‑592526
                                                                                                                                      |                                                                                                                                        
tttGATGTCCGCGCGTGCTGTAAAAATTTACGCTGCTCGCGTTCTTTGGCTTCGTCCACCACCGGAAAACGGACAAAAATTTCCGTCATACCTTTTTCTTTCAGGCGGAAGCCAATGTCGTAATCTTCAGTAAGGCTCTg                                                                                                                                     >  1:218176/1‑140 (MQ=255)
    aTGTCCGCGCGTGCTGTAAAAATTTACGCTGCTCGCGTTCTTTGGCTTCGTCCACCACCGGAAAACGGACAAAAATTTCCGTCATACCTTTTTCTTTCAGGCGGAAGCCAATGTCGTAATCTTCAGTAAGGCTCTGCACg                                                                                                                                 >  2:442704/1‑140 (MQ=255)
    aTGTCCGCGCGTGCTGTAAAAATTTACGCTGCTCGCGTTCTTTGGCTTCGTCCACCACCGGAAAACGGACAAAAATTTCCGTCATACCTTTTTATTTCAGGCGGAAGCCAATGTCGTAATCTTCAGTAAGGCTCTGCACg                                                                                                                                 >  1:625074/1‑140 (MQ=255)
        ccGCGCGTGCTGTAAAAATTTACGCTGCTCGCGTTCTTTGGCTTCGTCCACCACCGGAAAACGGACAAAAATTTCCGTCATACCTTTTTCTTTCAGGCGGAAGCCAATGTCGTAATCTTCAGTAAGGCTCTGCACGTCGa                                                                                                                             <  2:39907/140‑1 (MQ=255)
         cgcgcgTGCTGTAAAAATTTACGCTGCTCGCGTTCTTTGGCTTCGTCCACCACCGGAAAACGGACAAAAATTTCCGTCATACCTTTTTCTTTCAGGCGGAAGCCAATGTCGTAATCTTCAGTAAGGCTCTGCACGTCGaa                                                                                                                            >  2:1652/1‑140 (MQ=255)
         cgcgcgTGCTGTAAAAATTTACGCTGCTCGCGTTCTTTGGCTTCGTCCACCACCGGAAAACGGACAAAAATTTCCGTCATACCTTTTTCTCTCAGGCGGAAGCCAATGTCGTAATCTTCAGTAAGGCTCTg                                                                                                                                     <  1:361312/131‑1 (MQ=255)
         cgcgcgTGCTGTAAAAATTTACGCTGCTCGCGTTCTTTGGCTTCGTCCACCACCGGAAAACGGACAAAAATTTCCGTCATACCTTTTTCTCTCAGGCGGAAGCCAATGTCGTAATCTTCAGTAAGGCTCTg                                                                                                                                     >  2:361312/1‑131 (MQ=255)
              gTGCTGTAAAAATTTACGCTGCTCGCGTTCTTTGGCTTCGTCCACCACCGGAAAACGGACAAAAATTTCCGTCATACCTTTTTCTTTCAGGCGGAAGCCAATGTCGTAATCTTCAGTAAGGCTCTGCACGTCGAAAGCaa                                                                                                                       <  1:538059/140‑1 (MQ=255)
                    tAAAAATTTACGCTGCTCGCGTTCTTTGGCTTCGTCCACCACCGGAAAACGGACAAAAATTTCCGTCATACCTTTTTATTTCAGGCGGAAGCCAATGTCGTAATCTTCAGTAAGGCTCTGCACGTCGAAAGCAATaccgt                                                                                                                 <  2:625074/140‑1 (MQ=255)
                              cGCTGCTCGCGTTCTTTGGCTTCGTCCACCACCGGAAAACGGACAAAAATTTCCGTCATACCTTTTTCTTTCAGGCGGAAGCCAATGTCGTAATCTTCAGTAAGGCTCTGCACGTCGAAAGCAATACCGTCACCGTCAGc                                                                                                       >  2:75970/1‑140 (MQ=255)
                              cGCTGCTCGCGTTCTTTGGCTTCGTCCACCACCGGAAAACGGACAAAAATTTCCGTCATACCTTTTTCTTTCAGGCGGAAGCCAATGTCGTAATCTTCAGTAAGGCTCTGCACGTCGAAAGCAATACCGTCACCGTCAGc                                                                                                       >  1:420299/1‑140 (MQ=255)
                              cGCTGCTCGCGTTCTTTGGCTTCGTCCACCACCGGAAAACGGACAAAAATTTCCGTCATACCTTTTTCTTTCAGGCGGAAGCCAATGTCGTAATCTTCAGTAAGGCTCTGCACGTCGAAAGCAATACCGTCACCGTCAGc                                                                                                       >  1:543784/1‑140 (MQ=255)
                              cGCTGCTCGCGTTCTTTGGCTTCGTCCACCACCGGAAAACGGACAAAAATTTCCGTCATACCTTTTTCTTTCAGGCGGAAGCCAATGTCGTAATCTTCAGTAAGCCTCTGCACGTCGGAAGCGATACCGTCGCCGTCAGc                                                                                                       >  2:303740/1‑140 (MQ=255)
                                           cTTTGGCTTCGTCCACCACCGGAAAACGGACAAAAATTTCCGTCATACCTTTTTCTTTCAGGCGGAAGCCAATGTCGTAATCTTCAGTAAGGCTCTGCACGTCGAAAGCAATACAGTCACCGTCAGCTAACAGTGCGGTc                                                                                          >  2:640631/1‑140 (MQ=255)
                                                       ccaccaCCGGAAAACGGACAAAAATTTCCGTCATACCTTTTTCTTTCAGGCGGAAGCCAATGTCGTAATCTTCAGTAAGGCTCTGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGcggcgg                                                                              <  1:75970/140‑1 (MQ=255)
                                                              cGGAAAACGGACAAAAATTTCCGTCATACCTTTTTCTTTCAGGCGGAAGCCAATGTCGTAATCTTCAGTAAGGCTCTGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCTGAAAc                                                                       >  2:40463/1‑140 (MQ=255)
                                                               ggAAAACGGACAAAAATTTCCGTCATACCTTTTTCTTTCAGGCGGAAGCCAATGTCGTAATCTTCAGTAAGGCTCTGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCTGAAACa                                                                      >  1:631411/1‑140 (MQ=255)
                                                                gAAAACGGACAAAAATTTCCGTCATACCTTTTTCTTTCAGGCGGAAGCCAATGTCGTAATCTTCAGTAAGGCTCTGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCTGAAACAg                                                                     >  1:77388/1‑140 (MQ=255)
                                                                     cGGACAAAAATTTCCGTCATACCTTTTTCTTTCAGGCGGAAGCCAATGTCGTAATCTTCAGTAAGGCTCTGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCTGAAACAGGTGcc                                                                <  2:77388/140‑1 (MQ=255)
                                                                        aCAAAAATTTCCGTCATACCTTTTTCTTTCAGGCGGAAGCCAATGTCGTAATCTTCAGTAAGGCTCTGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCTGAAACAGGCGCCGAc                                                             >  1:412374/1‑140 (MQ=255)
                                                                                 tCCGTCATACCTTTTTCTTTCAGGCGGAAGCCAATGTCGTAATCTTCAGTAAGGCTCTGCACGTCGAAAGCAATACAGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCTGAAACAGGTGCCGACGCCTGCGCt                                                    <  1:640631/140‑1 (MQ=255)
                                                                                            tttttCTTTCAGGCGGAAGCCAATGTCGTAATCTTCAGTAAGGCTCTGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCTGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTc                                         >  1:373875/1‑140 (MQ=255)
                                                                                                     cAGGCGGAAGCCAATGTCGTAATCTTCAGTAAGGCTCTGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCTGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAggg                                >  1:605027/1‑140 (MQ=255)
                                                                                                     cAGGCGGAAGCCAATGTCGTAATCTTCAGTAAGGCTCTGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCTGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAggc                                >  1:331443/1‑139 (MQ=255)
                                                                                                            aaGCCAATGTCGTAATCTTCAGTAAGGCTCTGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCTGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTCACg                         <  1:653428/140‑1 (MQ=255)
                                                                                                            aaGCCAATGTCGTAATCTTCAGTAAGGCTCTGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCTGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTCACg                         <  1:131016/140‑1 (MQ=255)
                                                                                                            aaGCCAATGTCGTAATCTTCAGTAAGGCTCTGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCTGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTCACg                         <  1:442704/140‑1 (MQ=255)
                                                                                                             aGCCAATGTCGTAATGTTCAGTTAGTCTCGGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCTGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTCACGc                        <  1:303740/140‑1 (MQ=255)
                                                                                                                    gTCGTAATCTTCAGTAAGGCTCTGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCTGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTCACGCACCGGaa                 <  2:373875/140‑1 (MQ=255)
                                                                                                                    gTCGTAATCTTCAGTAAGGCTCAGCACGTCGAAAGCAATGCCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCTGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTCACGCACCGGaa                 <  2:218176/140‑1 (MQ=255)
                                                                                                                            cTTCAGTAAGGCTCTGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCTGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTCACGCACCGGAACATCTTTg         <  1:1652/140‑1 (MQ=255)
                                                                                                                                   aaGGCTCTGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCTGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTCACGCACCGGAACATCTTTGCCATGCa  <  2:605027/140‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                        
TTTGATGTCCGCGCGTGCTGTAAAAATTTACGCTGCTCGCGTTCTTTGGCTTCGTCCACCACCGGAAAACGGACAAAAATTTCCGTCATACCTTTTTCTTTCAGGCGGAAGCCAATGTCGTAATCTTCAGTAAGACTCTGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCTGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTCACGCACCGGAACATCTTTGCCATGCA  >  NC_000913/592256‑592526

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑