Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	590,706	C→T	L83L (CTG→CTA)	nfrA ←	bacteriophage N4 receptor, outer membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	590,706	0	C	T	100.0%	141.6 / NA	42	L83L (CTG→CTA)	nfrA	bacteriophage N4 receptor, outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (20/22);  total (20/22)

TCACTTCAACCGGAATAGCCGCCAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTCAGCGGAATATTATCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCGTCGTTATTTTTCTGTGCCTTCAGCGCTTTATCGAGACGGGGATAAATAACAAAATGGCGATAATCGCTCAGC  >  NC_000913/590570‑590841                                                                                                                                         |                                                                                                                                        tCACTTCAACCGGAATAGCCGCCAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTTAGc                                                                                                                                      <  1:464640/140‑1 (MQ=255)       cAACCGGAATAGCCGCCAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTTAGCGGAata                                                                                                                                <  1:1117269/140‑1 (MQ=255)       cAACCGGAATAGCCGCCAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTTAGCGGAata                                                                                                                                >  2:79738/1‑140 (MQ=255)       cAACCGGAATAGCCGCCAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTTAGCGGAata                                                                                                                                <  2:341084/140‑1 (MQ=255)                  gccgccAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTTAGCGGAATATTATCCGGCAc                                                                                                                     <  2:6331/140‑1 (MQ=255)                   ccgccAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTTAGCGGAATATTATCCGGCAcc                                                                                                                    >  1:771941/1‑140 (MQ=255)                    cgccAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTTAGCGGAATATTATCCGGCAcc                                                                                                                    >  2:1358296/1‑139 (MQ=255)                    cgccAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTTAGCGGAATATTATCCGGCAcc                                                                                                                    <  1:1358296/139‑1 (MQ=255)                     gccAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTTAGCGGAATATTATCCGGCACctg                                                                                                                  <  1:59223/140‑1 (MQ=255)                      ccAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTTAGCGGAATATTATCCGGCACctgc                                                                                                                 <  2:771941/140‑1 (MQ=255)                        aGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTTAGCGGAATATTATCCGGCACctgctg                                                                                                               >  1:1292723/1‑140 (MQ=255)                              cgcTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTTAGCGGAATATTATCCGGc                                                                                                                       <  2:236960/126‑1 (MQ=255)                              cgcTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTTAGCGGAATATTATCCGGc                                                                                                                       >  1:236960/1‑126 (MQ=255)                              cgcTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTTAGCGGAATATTATCCGGCACCTGCTGGTGtat                                                                                                         >  1:1179092/1‑140 (MQ=255)                                        gggCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTTAGCGGAATATTATCCGGCACctgctg                                                                                                               >  2:1103995/1‑124 (MQ=255)                                        gggCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTTAGCGGAATATTATCCGGCACctgctg                                                                                                               >  2:1241875/1‑124 (MQ=255)                                        gggCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTTAGCGGAATATTATCCGGCACctgctg                                                                                                               <  1:1241875/124‑1 (MQ=255)                                        gggCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTTAGCGGAATATTATCCGGCACctgctg                                                                                                               <  1:1103995/124‑1 (MQ=255)                                         ggCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTTAGCGGAATATTATCCGGCACCTGCTGGTGTATATATTCAAAtt                                                                                              >  1:261372/1‑140 (MQ=255)                                           cATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTTAGCGGAATATTATCCGGCACCTGCTGGTGTATATATTCAAATTcg                                                                                            <  2:1292723/140‑1 (MQ=255)                                               tcCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTTAGCGGAATATTATCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGa                                                                                        >  1:1331912/1‑140 (MQ=255)                                                    ggTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTTAGCGGAATATTATCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCg                                                                                   >  1:1360970/1‑140 (MQ=255)                                                         cGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTTAGCGGAATATTATCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGc                                                                              >  2:442629/1‑140 (MQ=255)                                                          gTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTTAGCGGAATATTATCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCt                                                                             <  1:79738/140‑1 (MQ=255)                                                           tttCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTTAGCGGAATATTATCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCtt                                                                            >  1:1151734/1‑140 (MQ=255)                                                              cAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTTAGCGGAATATTATCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTtcgt                                                                         <  2:1344830/140‑1 (MQ=255)                                                              cAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTTAGCGGAATATTATCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTtcgt                                                                         <  2:1360970/140‑1 (MQ=255)                                                                    aTCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTTAGCGGAATATTATCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCGTCGTTAt                                                                   >  1:1348577/1‑140 (MQ=255)                                                                                 aGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTTAGCGGAATATTATCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCGTCGTTATTTTTCTGTGCCtt                                                      >  2:1330425/1‑140 (MQ=255)                                                                                      cgcCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTTAGCGGAATATTATCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCGTCGTTATTTTTCTGTGCCTTCAgcg                                                 >  1:716962/1‑140 (MQ=255)                                                                                      cgcCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTTAGCGGAATATTATCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCGTCGTTATTTTTCTGTGCCTTCAgcg                                                 <  1:1330425/140‑1 (MQ=255)                                                                                            gTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTTAGCGGAATATTATCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCGTCGTTATTTTTCTGTGCCTTCAGCGCTTTAt                                           >  2:1155035/1‑140 (MQ=255)                                                                                            gTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTTAGCGGAATATTATCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCGTCGTTATTTTTCTGTGCCTTCAGCGCTTTAt                                           >  2:1137978/1‑140 (MQ=255)                                                                                            gTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTTAGCGGAATATTATCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCGTCGTTATTTTTCTGTGCCTTCAGCGCTTTAt                                           <  2:273668/140‑1 (MQ=255)                                                                                                     aCCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTTAGCGGAATATTATCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCGTCGTTATTTTTCTGTGCCTTCAGCGCTTTATCGAGACggg                                  >  2:579358/1‑140 (MQ=255)                                                                                                      ccAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTTAGCGGAATATTATCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCGTCGTTATTTTTCTGTGCCTTCAGCGCTTTATCGAGACgggg                                 <  2:1331912/140‑1 (MQ=255)                                                                                                           aTGGCGATAGGCTTCCGCAAGGTATAAAGTTAGCGGAATATTATCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCGTCGTTATTTTTCTGTGCCTTCAGCGCTTTATCGAGACGGGGataaa                            <  2:1160210/140‑1 (MQ=255)                                                                                                           aTGGCGATAGGCTTCCGCAAGGTATAAAGTTAGCGGAATATTATCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCGTCGTTATTTTTCTGTGCCTTCAGCGCTTTATCGAGACGGGGataaa                            <  2:1179092/140‑1 (MQ=255)                                                                                                              gCGATAGGCTTCCGCAAGGTATAAAGTTAGCGGAATATTATCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCGTCGTTATTTTTCTGTGCCTTCAGCGCTTTATCGAGACgg                                   <  2:1283401/130‑1 (MQ=255)                                                                                                              gCGATAGGCTTCCGCAAGGTATAAAGTTAGCGGAATATTATCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCGTCGTTATTTTTCTGTGCCTTCAGCGCTTTATCGAGACgg                                   >  1:1283401/1‑130 (MQ=255)                                                                                                                  tAGGCTTCCGCAAGGTATAAAGTTAGCGGAATATTATCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCGTCGTTATTTTTCTGTGCCTTCAGCGCTTTATCGAGACGGGGATAAATAACaaa                     <  2:1151734/140‑1 (MQ=255)                                                                                                                                     aaGTTAGCGGAATATTATCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCGTCGTTATTTTTCTGTGCCTTCAGCGCTTTATCGAGACGGGGATAAATAACAAAATGGCGATAATCGCTCAGc  <  2:1348577/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                        TCACTTCAACCGGAATAGCCGCCAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTCAGCGGAATATTATCCGGCACCTGCTGGTGTATATATTCAAATTCGCGGATGGCGGTTGCTTCGTCGTTATTTTTCTGTGCCTTCAGCGCTTTATCGAGACGGGGATAAATAACAAAATGGCGATAATCGCTCAGC  >  NC_000913/590570‑590841

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