Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 592,150 T→G K343N (AAA→AAC 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,1500TG100.0% 66.4 / NA 20K343N (AAA→AACnfrBbacteriophage N4 receptor, inner membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (10/10);  total (10/10)

GGCCACAAACTTTCATACGCCAGCAACAGCAAAAGCTGGATCATCACCAGCATCGCGAGGAAGCTGACAAAGTTACTGATTGCCCCTTTGCGGTCGCGCCAGAGAAAGTAGTTCAGCGTCAGGCTGGAGGTCCATTTATGGGTTTTAAAGCCTTGGAAAACAATGCCGATGATCCAGCGGGATTTTTGTCGAACCGCAGTCGAAAAGGTATCGGGGAAATATTCGCGCACGCAGATCATGTTTGATGTCCGCGCGTGCTGTAAAAA  >  NC_000913/592016‑592281
                                                                                                                                      |                                                                                                                                   
ggCCACAAACTTTCATACGCCAGCAACAGCAAAAGCTGGATCATCACCAGCATCGCGAGGAAGCTGACAAAGTTACTGATTGCCCCTTTGCGGTCGCGCCAGAGAAAGTAGTTCAGCGTCAGGCTGGAGGTCCAGTTATg                                                                                                                                <  1:736972/140‑1 (MQ=255)
                  gCCAGCAACAGCAAAAGGTGGACCATCACCAGCATGGCGAGGAAGCTGACAAAGTTACTGATTGCCCCTTTGCGGTCGCGCCAGAGAAAGTAGTTCAGCGTCAGGCTGGAGGTCCAGTTATGGGTTTTAAAGCCTTGGaa                                                                                                              <  2:332723/140‑1 (MQ=255)
                        aacagcaaAAGCTGGATCATCACCAGCATCGCGAGGAAGCTGACAAAGTTACTGATTGCCCCTTTGCGGTCGCGCCAGAGAAAGTAGTTCAGCGTCAGGCTGGAGGTCCAGTTATGGGTTTTAAAGCCTTGGAAAACAAt                                                                                                        <  1:539337/140‑1 (MQ=255)
                        aacagcaaAAGCTGGATCATCACCAGCATCGCGAGGAAGCTGACAAAGTTACTGATTGCCCCTTTGCGGTCGCGCCAGAGAAAGTAGTTCAGCGTCAGGCTGGAGGTCCAGTTATGGGTTTTAAAGCCTTGGAAAACAAt                                                                                                        <  1:786625/140‑1 (MQ=255)
                                            caCCAGCATCGCGAGGAAGCTGACAAAGTTACTGATTGCCCCTTTGCGGTCGCGCCAGAGAAAGTAGTTCAGCGTCAGGCTGGAGGTCCAGTTATGGGTTTTAAAGCTTGGAAAACAATGCCGATGATCCAGCGGGAttt                                                                                   >  1:551956/1‑140 (MQ=255)
                                                aGCATCGCGAGGAAGCTGACAAAGTTACTGATTGCCCCTTTGCGGTCGCGCCAGAGAAAGTAGTTCAGCGTCAGGCTGGAGGTCCAGTTATCGGTTTTAAAGCCATGGAAAACAATGCCGATCATCCAGCGGGATTTTTg                                                                                >  2:325583/1‑140 (MQ=255)
                                                           gAAGCTGACAAAGTTGCTGATTTCCTCTTTGCGGTCGTGCCAGAGAAAGTAGTTCAGCGTCAGGCTGGAGGTCCAGTTATGGGTTTTAAAGCCTTGGAAAACAATGCCGATGATCCAGCGGGATTTTTGTCGAACCGCAg                                                                     <  1:340147/140‑1 (MQ=255)
                                                                                 gCCCCTTTGCGGTCGTGCCAGACAAAGTATTTCAGCGACAGGCTGGAGGTCCAGTTATGGGTTTTAAAGCCTTGGAAAACAATGCCGATGATCCAGCGGGATTTTTGTCGAACCGCAGTCGAAAAGGTATCGGGGAAata                                               <  1:325583/140‑1 (MQ=255)
                                                                                  ccccTTTGCGGTCGCGCCAGAGAAAGTAGTTCAGCGTCAGGCTGGAGGTCCAGTTATGGGTTTTAAAGCCTTGGAAAACAATGCCGATGATCCAGCGGGATTTTTGTCGAACCGCAGTCGAAAAGGTATCGGGGAAatat                                              >  1:534035/1‑140 (MQ=255)
                                                                                  ccccTTTGCGGTCGCGCCAGAGAAAGTAGTTCAGCGTCAGGCTGGAGGTCCAGTTATGGGTTTTAAAGCCTTGGAAAACAATGCCGATGATCCAGCGGGATTTTTGTCGAACCGCAGTCGAAAAGGTATCGGGGAAatat                                              >  2:816563/1‑140 (MQ=255)
                                                                                   cccTTTGCGGTCGCGCCAGAGAAAGTAGTTCAGCGTCAGGCTGGAGGTCCAGTTATGGGTTTTAAAGCCTTGGAAAACAATGCCGATGATCCAGCGGGATTTTTGTCGAACCGCAGTCGAAAAGGTATCGGGGAAATAtt                                             >  2:855465/1‑140 (MQ=255)
                                                                                        tGCGGTCGCGCCAGAGAAAGTAGTTCAGCGTCAGGCTGGAGGTCCAGTTATGGGTTTTAAAGCCTTGGAAAACAATGCCGATGATCCAGCGGGATTTTTGTCGAACCGCAGTCGAAAAGGTAt                                                         <  2:546408/123‑1 (MQ=255)
                                                                                        tGCGGTCGCGCCAGAGAAAGTAGTTCAGCGTCAGGCTGGAGGTCCAGTTATGGGTTTTAAAGCCTTGGAAAACAATGCCGATGATCCAGCGGGATTTTTGTCGAACCGCAGTCGAAAAGGTAt                                                         >  1:546408/1‑123 (MQ=255)
                                                                                         gCGGTCGCGCCAGAGAAAGTAGTTCAGCGTCAGGCTGGAGGTCCAGTTATGGGTTTTAAAGCCTTGGAAAACAATGCCGATGATCCAGCGGGATTTTTGTCGAACCGCAGTCGAAAAGGTATCGGGGAAATATTCGcgca                                       <  1:855465/140‑1 (MQ=255)
                                                                                                cgcCAGAGAAAGTAGTTCAGCGTCAGGCTGGAGGTCCAGTTATGGGTTTTAAAGCCTTGGAAAACAATGCCGATGATCCAGCGGGATTTTTGTCGAACCGCAGTCGAAAAGGTATCGGGGAAATATTCGCGCACGCAGAt                                <  1:816563/140‑1 (MQ=255)
                                                                                                      agaAAGTAGTTCAGCGTCAGGCTGGAGGTCCAGTTATGGGTTTTAAAGCCTTGGAAAACAATGCCGATGATCCAGCGGGATTTTTGTCGAACCGCAGTCGAAAAGGTATCGGGGAAATATTCGCGCACGCAGATCATGtt                          <  2:534035/140‑1 (MQ=255)
                                                                                                               ttCAGCGTCAGGCTGGAGGTCCAGTTATGGGTTTTAAAGCCTTGGAAAACAATGCCGATGATCCAGCGGGATTTTTGTCGAACCGCAGTCGAAAAGGTATCGGGGAAATATTCGCGCACGCAGATCATGTTTGATGTCcg                 >  2:450730/1‑140 (MQ=255)
                                                                                                               ttCAGCGTCAGGCTGGAGGTCCAGTTATGGGTTTTAAAGCCTTGGAAAACAATGCCGATGATCCAGCGGGATTTTTGTCGAACCGCAGTCGAAAAGGTATCGGGGAAATATTCGCGCACGCAGATCATGTTTGATGTCcg                 >  1:123757/1‑140 (MQ=255)
                                                                                                                         ggCTGGAGGTCCAGTTATGGGTTTTAAAGCCTTGGAAAACAATGCCGATGATCCAGCGGGATTTTTGTCGAACCGCAGTCGAAAAGGTATCGGGGAAATATTCGCGCACGCAGATCATGTTTGATGTCCGCGCGTGCTGt       >  1:549225/1‑140 (MQ=255)
                                                                                                                              gAGGTCCAGTTATGGGTTTTAAAGCCTTGGAAAACAATGCCGATGATCCAGCGGGATTTTTGTCGAACCGCAGTCGAAAAGGTATCGGGGAAATATTCGCGCACGCAGATCATGTTTGATGTCCGCGCGTGCTGTaaaaa  >  2:200451/1‑140 (MQ=255)
                                                                                                                                      |                                                                                                                                   
GGCCACAAACTTTCATACGCCAGCAACAGCAAAAGCTGGATCATCACCAGCATCGCGAGGAAGCTGACAAAGTTACTGATTGCCCCTTTGCGGTCGCGCCAGAGAAAGTAGTTCAGCGTCAGGCTGGAGGTCCATTTATGGGTTTTAAAGCCTTGGAAAACAATGCCGATGATCCAGCGGGATTTTTGTCGAACCGCAGTCGAAAAGGTATCGGGGAAATATTCGCGCACGCAGATCATGTTTGATGTCCGCGCGTGCTGTAAAAA  >  NC_000913/592016‑592281

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