Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	590,771	C→T	D62N (GAC→AAC)	nfrA ←	bacteriophage N4 receptor, outer membrane subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	590,771	0	C	T	100.0%	166.9 / NA	43	D62N (GAC→AAC)	nfrA	bacteriophage N4 receptor, outer membrane subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/21);  total (22/21)

GACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTCAGCGGAATATTATCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCGTCGTTATTTTTCTGTGCCTTCAGCGCTTTATCGAGACGGGGATAAATAACAAAATGGCGATAATCGCTCAGCCCCAGCTCTTCTGCGCTGGTGCCGATATTGTCTGCGAGTGCGCTGGTACTCAATAAAGACGTCAGCAGTAAACC  >  NC_000913/590624‑590915                                                                                                                                                    |                                                                                                                                                 gACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTCAGCGGAATATTATCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCgttgtt                                                                                                                                               >  1:158045/1‑151 (MQ=255)  aCGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTCAGCGGAATATTATCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCgttgtt                                                                                                                                               >  1:623430/1‑150 (MQ=255)                   cAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTCAGCGGAATATTATCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCGTTGTTATTTTTCTGTGCCTTc                                                                                                                               <  2:183678/149‑1 (MQ=255)                        aaCAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGCATAAAGTCAGCGGAATATTATCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCGTTGTTATTTTTCTGTGCCTTCAGCGCtt                                                                                                                        >  2:679676/1‑151 (MQ=255)                         aCAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTCAGCGGAATATTATCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCGTTGTTATTTTTCTGTGCCTTCAGCGCttt                                                                                                                       <  2:211773/151‑1 (MQ=255)                            gCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTCAGCGGAATATTATCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCGTTGTTATTTTTCTGTGCCTTCAGCGCTTTATc                                                                                                                    >  2:448296/1‑151 (MQ=255)                               gcgcCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTCAGCGGAATATTATCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCGTTGTTATTTTTCTGTGCCTTCAGCGCTTTATCg                                                                                                                   >  2:158605/1‑149 (MQ=255)                                 gcCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTCAGCGGAATATTATCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCGTTGTTATTTTTCTGTGCCTTCAGGGCTTTATCGAGAc                                                                                                               >  2:31868/1‑151 (MQ=255)                                    cGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTCAGCGGAATATTATCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCGTTGTTATTTTTCTGTGCCTTCAGCGCTTTATCGAGACggg                                                                                                            >  2:574788/1‑151 (MQ=255)                                     ggTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTCAGCGGAATATTATCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCGTTGTTATTTTTCTGTGCCTTCAGCGCTTTATCGAGACgggg                                                                                                           >  1:27771/1‑151 (MQ=255)                                            atGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTCAGCGGAATATTATCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCGTTGTTATTTTTCTGTGCCTTCAGCGCTTTATCGAGACGGGGataaata                                                                                                    <  2:333596/151‑1 (MQ=255)                                            atGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTCAGCGGAATATTATCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCGTTGTTATTTTTCTGTGCCTTCAGCGCTTTATCGAGACGGGGataaat                                                                                                     <  2:213643/150‑1 (MQ=255)                                              gACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTCAGCGGAATATTATCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCGTTGTTATTTTTCTGTGCCTTCAGCGCTTTATCGAGACGGGGATAAATAAc                                                                                                  >  1:323924/1‑151 (MQ=255)                                                ccAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTCAGCGGAATATTATCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCGTTGTTATTTTTCTGTGCCTTCAGCGCTTTATCGAGACGGGGATAAATAACaa                                                                                                <  2:462140/151‑1 (MQ=255)                                                  aaaaTGGCGATAGGCTTCCGCAAGGTATAAAGTCAGCGGAATATTATCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCGTTGTTATTTTTCTGTGCCTTCAGCGCTTTATCGAGACGGGGATAAATAACaaaa                                                                                              <  1:674906/151‑1 (MQ=255)                                                       ggCGATAGGCTTCCGCAAGGTATAAAGTCAGCGGAATATTATCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCGTTGTTATTTTTCTGTGCCTTCAGCGCTTTATCGAGACGGGGATAAATAACAAAATGGc                                                                                          <  1:509472/150‑1 (MQ=255)                                                              ggCTTCCGCAAGGTATAAAGTCAGCGGAATATTATCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCGTTGTTATTTTTCTGTTCCTTCAGCGCTTTATCGAGACGGGGATAAATAACAAAATGGCGATAAt                                                                                    >  2:339294/1‑149 (MQ=255)                                                               gCTTCCGCAAGGTATAAAGTCAGCGGAATATTATCTGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCGTTGTTATTTTTCTGTGCCTTCAGCGCTTTATCGAGACGGGGATAAATAACAAAATGGCGATAATCGc                                                                                 <  2:281363/151‑1 (MQ=255)                                                               gCTTCCGCAAGGTATAAAGTCAGCGGAATATTATCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCGTTGTTATTTTTCTGTGCCTTCAGCGCTTTATCGAGACGGGGATAAATAACAAAATGGCGATAATCGc                                                                                 <  2:356373/151‑1 (MQ=255)                                                                     gCAAGGTATAAAGTCAGCGGAATATTATCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCGTTGTTATTTTTCTGTGCCTTCAGCGCTTTATCGAGACGGGGATAAATAACAAAATGGCGATAATCGCTCAGcc                                                                           <  1:53729/151‑1 (MQ=255)                                                                        aGGTATAAAGTCAGCGGAATATTATCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCGTTGTTATTTTTCTGTGCCTTCAGCGCTTTATCGAGACGGGGATAAATAACAAAATGGCGATAATCGCTCAGcccc                                                                         >  1:29024/1‑150 (MQ=255)                                                                              aaaGTCAGCGGAATATTATCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCGTTGTTATTTTTCTGTGCCTTCAGCGCTTTATCGAGACGGGGATAAATAACAAAATGGCGATAATCGCTCAGCCCCAGctct                                                                   <  1:205138/150‑1 (MQ=255)                                                                              aaaGTCAGCGGAATATTATCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCGTTGTTATTTTTCTGTGCCTTCAGCGCTTTATCGAGACGGGGATAAATAACAAAATGGCGATAATCGCTCAGCCCCAGctct                                                                   <  1:115001/150‑1 (MQ=255)                                                                                      cGGAATATTATCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCGTTGTTATTTTTCTGTGCCTTCAGCGCTTTATCGAGACGGGGATAAATAACAAAATGGCGATAATCGCTCAGCCCCAGCTCTTCTGCGCTg                                                          >  1:556969/1‑151 (MQ=255)                                                                                       ggAATATTATCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCGTTGTTATTTTTCTGTGCCTTCAGCGCTTTATCGAGACGGGGATAAATAACAAAATGGCGATAATCGCTCAGCCCCAGCTCTTCTGCGCTgg                                                         <  2:438119/151‑1 (MQ=255)                                                                                        gAATATTATCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCGTTGTTATTTTTCTGTGCCTTCAGCGCTTTATCGAGACGGGGATAAATAACAAAATGGCGATAATCGCTCAGCCCCAGCTCTTCTGCGCTGGt                                                        <  2:655569/151‑1 (MQ=255)                                                                                         aaTATTATCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCGTTGTTATTTTTCTGTGCCTTCAGCGCTTTATCGAGACGGGGATAAATAACAAAATGGCGATAATCGCTCAGCCCCAGCTCTTCTGCGCTGGTg                                                       <  2:283544/151‑1 (MQ=255)                                                                                          atatTATCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCGTTGTTATTTTTCTGTGCCTTCAGCGCTTTATCGAGACGGGGATAAATAACAAAATGGCGATAATCGCTCAGCCCCAGCTCTTCTGCGCTGGTGc                                                      >  1:14157/1‑151 (MQ=255)                                                                                            attatCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCGTTGTTATTTTTCTGTGCCTTCAGCGCTTTATCGAGACGGGGATAAATAACAAAATGGCGATAATCGCTCAGCCCCAGCTCTTCTGCGCTGGTGCCg                                                    <  2:170537/151‑1 (MQ=255)                                                                                             ttatCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCGTTGTTATTTTTCTGTGCCTTCAGCGCTTTATCGAGACGGGGATAAATAACAAAATGGCGATAATCGCTCAGCCCCAGCTCTTCTGCGCTGGTg                                                       >  2:43519/1‑147 (MQ=255)                                                                                             ttatCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCGTTGTTATTTTTCTGTGCCTTCAGCGCTTTATCGAGACGGGGATAAATAACAAAATGGCGATAATCGCTCAGCCCCAGCTCTTCTGCGCTGGTGCCg                                                    >  1:656530/1‑150 (MQ=255)                                                                                                   ggCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCGTTGTTATTTTTCTGTGCCTTCAGCGCTTTATCGAGACGGGGATAAATAACAAAATGGCGATAATCGCTCAGCCCCAGCTCTTCTGCGCTGGTGCCGATATTGt                                             <  2:147157/151‑1 (MQ=255)                                                                                                   ggCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCGTTGTTATTTTTCTGTGCCTTCAGCGCTTTATCGAGACGGGGATAAATAACAAAATGGCGATAATCGCTCAGCCCCAGCTCTTCTGCGCTGGTGCCGATATTGt                                             <  1:156041/151‑1 (MQ=255)                                                                                                       ccTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCGTTGTTATTTTTCTGTGCCTTCAGCGCTTTATCGAGACGGGGATAAATAACAAAATGGCGATAATCGCTCAGCCCCAGCTCTTCTGCGCTGGTGCCGATATTGTCTGc                                         <  2:519484/151‑1 (MQ=255)                                                                                                            tgGTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCGTTGTTATTTTTCTGTGCCTTCAGCGCTTTATCGAGACGGGGATAAATAACAAAATGGCGATAATCGCTCAGCCCCAGCTCTTCTGCGCTGGTGCCGATATTGTCTGCGAGTg                                    >  1:100234/1‑151 (MQ=255)                                                                                                             ggTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCGTTGTTATTTTTCTGTGCCTTCAGCGCTTTATCGAGACGGGGATAAATAACAAAATGGCGATAATCGCTCAGCCCCAGCTCTTCTGCGCTGGTGCCGATATTGTCTGCGAGTgc                                   >  1:504013/1‑151 (MQ=255)                                                                                                              gtgtATATATTCAAATTCGCGGATGGCGGTTGCTTCGTTGTTATTTTTCTGTGCCTTCAGCGCTTTATCGAGACGGGGATAAATAACAAAATGGCGATAATCGCTCAGCCCCAGCTCTTCTGCGCTGGTGCCGATATTGTCTGCGAGTgcg                                  >  1:516558/1‑151 (MQ=255)                                                                                                                  atatatTCAAATTCGCGGATGGCGGTTGCTTCGTTGTTATTTTTCTGTGCCTTCAGCGCTTTATCGAGACGGGGATAAATAACAAAATGGCGATAATCGCTCAGCCCCAGCTCTTCTGCGCTGGTGCCGATATTGTCTGCGAGTGCGCTg                               <  1:354692/150‑1 (MQ=255)                                                                                                                    atatTCAAATTCGCGGATGGCGGTTGCTTCGTTGTTATTTTTCTGTGCCTTCAGCGCTTTATCGAGACGGGGATAAATAACAAAATGGCGATAATCGCTCAGCCCCAGCTCTTCTGCGCTGGTGCCGATATTGTCTGCGAGTGCGCTGGt                             >  2:486190/1‑150 (MQ=255)                                                                                                                        tCAAATTCGCGGATGGCGGTTGCTTCGTTGTTATTTTTCTGTGCCTTCAGCGCTTTATCGAGACGGGGATAAATAACAAAATGGCGATAATCGCTCAGCCCCAGCTCTTCTGCGCTGGTGCCGATATTGTCTGCGAGTGCGCTGGTACTc                         >  2:244883/1‑150 (MQ=255)                                                                                                                           aaTTCGCGGATGGCGGTTGCTTCGTTGTTATTTTTCTGTGCCTTCAGCGCTTTATCGAGACGGGGATAAATAACAAAATGGCGATAATCGCTCAGCCCCAGCTCTTCTGCGCTGGTGCCGATATTGTCTGCGAGTGCGCTGGTACTCAAt                      >  1:449129/1‑150 (MQ=255)                                                                                                                                        cGGTTGCTTCGTTGTTATTTTTCTGTGCCTTCAGCGCTTTATCGAGACGGGGATAAATAACAAAATGGCGATAATCGCTCAGCCCCAGCTCTTCTGCGCTGGTGCCGATATTGTCTGCGAGTGCGCTGGTACTCAATAAAGACGTcagcag        <  2:529144/151‑1 (MQ=255)                                                                                                                                             gCTTCGTTGTTATTTTTCTGTGCCTTCAGCGCTTTATCGAGACGGGGATAAATAACAAAATGGCGATAATCGCTCAGCCCCAGCTCTTCTGCGCTGGTGCCGATATTGTCTGCGAGTGCGCTGGTACTCAATAAAGACGTCAGCAGTAAAc   >  1:500684/1‑151 (MQ=255)                                                                                                                                             gCTTCGTTGTTATTTTTCTGTGCCTTCAGCGCTTTATCGAGACGGGGATAAATAACAAAATGGCGATAATCGCTCAGCCCCAGCTCTTCTGCGCTGGTGCCGATATTGTCTGCGAGTGCGCTGGTACTCAATAAAGACGTCAGCAGTAAAc   >  2:93209/1‑151 (MQ=255)                                                                                                                                               ttCGTTGTTATTTTTCTGTGCCTTCAGCGCTTTATCGAGACGGGGATAAATAACAAAATGGCGATAATCGCTCAGCCCCAGCTCTTCTGCGCTGGTGCCGATATTGTCTGCGAGTGCGCTGGTACTCAATAAAGACGTCAGCAGTAAAcc  >  2:665240/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                 GACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTCAGCGGAATATTATCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCGTCGTTATTTTTCTGTGCCTTCAGCGCTTTATCGAGACGGGGATAAATAACAAAATGGCGATAATCGCTCAGCCCCAGCTCTTCTGCGCTGGTGCCGATATTGTCTGCGAGTGCGCTGGTACTCAATAAAGACGTCAGCAGTAAACC  >  NC_000913/590624‑590915

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑