Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 588,641 A→G Y772H (TAT→CAT)  nfrA ← bacteriophage N4 receptor, outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913588,6410AG97.4% 122.9 / ‑4.1 38Y772H (TAT→CAT) nfrAbacteriophage N4 receptor, outer membrane protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (20/17);  minor base T (1/0);  total (21/17)
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.50e-01

TTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGACCTCCTCATGCAAACGGCGGAAATTGAA  >  NC_000913/588511‑588773
                                                                                                                                  |                                                                                                                                    
gagacaggCATCATCACCCCGTGTTCTCCGGTATCCGCAGAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATGGCTTTTccc                                                                                                                             <  2:289476/134‑1 (MQ=255)
        catcatCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATGGCTTTTCCCTGGCGCTg                                                                                                                     >  1:813132/1‑140 (MQ=255)
              caCCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGTTACGATGGCTTTTCCCTGGCGCTGcgccgc                                                                                                               <  1:830414/140‑1 (MQ=255)
              caCCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATGGCTTTTCCCTGGCGCTGcgccgc                                                                                                               <  1:666644/140‑1 (MQ=255)
                 cccGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATTGCTTTTCCCTGGCGCTGCGCCGCCGa                                                                                                            >  2:642435/1‑140 (MQ=255)
                 cccGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATGGCTTTTCCCTGGCGCTGCGCCGCCGa                                                                                                            >  1:616258/1‑140 (MQ=255)
                  ccGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGGAGCTACGGTGGTTTGTCCCTGGCGCTGcgcctgtctc                                                                                                           >  1:289476/1‑134 (MQ=255)
                    gTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATGGCTTTTCCCTGGCGCTGCGCCGCCGACat                                                                                                         >  1:578735/1‑140 (MQ=255)
                     ttttCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAACATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATGGCTTTTCCCTGGCGCTGCGCCGCCGACatt                                                                                                        >  1:187490/1‑140 (MQ=255)
                          tcCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATGGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGt                                                                                                   <  1:517159/140‑1 (MQ=255)
                          tcCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATGGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGt                                                                                                   <  2:616258/140‑1 (MQ=255)
                          tcCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATGGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGt                                                                                                   <  2:813132/140‑1 (MQ=255)
                                 tCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATGGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGt                                                                                            <  1:752498/140‑1 (MQ=255)
                                        aGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATGGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCAtt                                                                                     >  1:475894/1‑140 (MQ=255)
                                                   taAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATGGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACg                                                                          <  2:475894/140‑1 (MQ=255)
                                                       cTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATGGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAg                                                                      >  1:464904/1‑140 (MQ=255)
                                                          agagCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATGGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCg                                                                   >  2:205916/1‑140 (MQ=255)
                                                          agagCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATGGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCg                                                                   >  2:577900/1‑140 (MQ=255)
                                                          agagCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATGGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCg                                                                   >  1:106801/1‑140 (MQ=255)
                                                           gagCAGGTCGCCTTCCAGCAGCCTATTGCGTGCGAAGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATGGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGa                                                                  <  2:347592/140‑1 (MQ=255)
                                                                 gTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATGGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGgaag                                                            >  2:718644/1‑140 (MQ=255)
                                                                      cTTCCAGCAGCATATTGCGTCCGATGCGGTACTCAGCTTCCAGTTGTCCGTAGCTACGATGGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCg                                                       >  1:420/1‑140 (MQ=255)
                                                                                  tatTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATGGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGTAAGCCGAGGGAAAAATCGAAACTGAACGTc                                           >  2:380846/1‑140 (MQ=255)
                                                                                  tatTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATGGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTc                                           >  2:513176/1‑140 (MQ=255)
                                                                                          ccGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATGGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAgcggcg                                   >  1:182343/1‑140 (MQ=255)
                                                                                           cGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATGGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGcga                                  <  2:182343/140‑1 (MQ=255)
                                                                                           cGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATGGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGcga                                  <  2:187490/140‑1 (MQ=255)
                                                                                           cGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATGGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGcga                                  <  1:718644/140‑1 (MQ=255)
                                                                                              tGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATGGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGcgaccg                               >  1:740637/1‑140 (MQ=255)
                                                                                              tGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATGGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCCGAACGCAAGCCGATGGCAGAATCGACACTGAACGTCCAGCGGcgaccg                               >  2:257900/1‑140 (MQ=255)
                                                                                               gCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATGGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGcgaccga                              <  2:740637/140‑1 (MQ=255)
                                                                                                 ggTACTCGGCTTCCAGTTGTCCGTATCTATGATGGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCgaccgacc                            <  1:380846/140‑1 (MQ=255)
                                                                                                          cTTCCAGTTGTCCGTAGCTACGATGGCTTTTCCCTGGCGCTGCGCCGCCGGCATTATTGTTAGCGGTACTCCTTGCGCCGGAACGCAAGCCGAAGGGAGAATCGAAAATGAACGACCACCGGCGACCGACCTCCTCactg                   >  2:743806/1‑137 (MQ=255)
                                                                                                          cTTCCAGTTGTCCGTAGCTACGATGGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGACCTCCTCATGc                   >  2:824769/1‑140 (MQ=255)
                                                                                                           ttCCAGTTGTCCGTAGCTACGACGGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGACCTCCTCATGCa                  <  1:287771/140‑1 (MQ=255)
                                                                                                             ccAGTTGTCCGTAGCTACGATGGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGACCTCCTCATGCaaa                >  2:589467/1‑140 (MQ=255)
                                                                                                                   gTCCGTAGCTACGATGGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTAAACCTCCAGCGGCGACAGACCCCCTCATGCAAAcggcgg          >  2:287771/1‑140 (MQ=255)
                                                                                                                           cTACGATGGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGACCTCCTCATGCAAACGGCGGAAATTGaa  <  2:420/140‑1 (MQ=255)
                                                                                                                           cTACGATGGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGACCTCCTCATGCAAACGGCGGAAATTGaa  <  1:792135/140‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                    
TTCACCGGCATCATCACCCCGTTTTCTCCGGTATCGGCAAAGACGCGGCTATAAACTGAGAGCAGGTCGCCTTCCAGCAGCATATTGCGTCCGATGCGGTACTCGGCTTCCAGTTGTCCGTAGCTACGATAGCTTTTCCCTGGCGCTGCGCCGCCGACATTATTGTTAGCGGTACTCATTGCGCCGGAACGCAAGCCGATGGAAGAATCGAAACTGAACGTCCAGCGGCGACCGACCTCCTCATGCAAACGGCGGAAATTGAA  >  NC_000913/588511‑588773

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