Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	590,608	C→A	5.9%	R116L (CGA→CTA)	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,608	0	C	A	5.9%	241.2 / 2.6	68	R116L (CGA→CTA)	nfrA	bacteriophage N4 receptor, outer membrane subunit
Reads supporting (aligned to +/- strand):  ref base C (26/38);  new base A (2/2);  total (28/40)
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.36e-01

CAGGGCATTCTGCCCGACTTCACTGCGACAACGCAGGGTCGGCGCAGCATCGCACGCTTTTTGCTGGGCAAGCAGTTCTTCAACTGTCGTAACGCTTTTCACTTCAACCGGAATAGCCGCCAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTCAGCGGAATATTATCCGGCACCTGCTGGTGTATATATT  >  NC_000913/590472‑590743                                                                                                                                         |                                                                                                                                        cAGGGCATTCTGCCCGACTTCACTGAGACAACGCAGGGTCGGCGCAGCATCGCACGCTTTTTGCTGGGCAAGCAGTTCTTCCACTGTCGTAACGCTTTTCACTTCAAACGGAATAGCCGCCAGACTGCGCTCAAGTCggg                                                                                                                                      <  2:628129/140‑1 (MQ=255)  aGGGCATTCTGCCCGACTTCACTGCGACAACGCAGGGTCGGCGCAGCATCGCACGCTTTTTGCTGGGCAAGCAGTTCTTCAACTGTCGTAACGCTTTTCACTTCAACCGGAATAGCCGCCAGACTGCGCTCAAGTCGGGc                                                                                                                                     <  2:169575/140‑1 (MQ=255)  aGGGCATTCTGCCCGACTTCACTGCGACAACGCAGGGTCGGCGCAGCATCGCACGCTTTTTGCTGGGCAAGCAGTTCTTCAACTGTCGTAACGCTTTTCACTTCAACCGGAATAGCCGCCAGACTGCGCTCAAGTCGGGc                                                                                                                                     <  1:653475/140‑1 (MQ=255)  aGGGCATTCTGCCCGACTTCACTGCGACAACGCAGGGTCGGCGCAGCATCGCACGCTTTTTGCTGGGCAAGCAGTTCTTCAACTGTCGTAACGCTTTTCACTTCAACCGGAATAGCCGCCAGACTGCGCTCAAGTCGGGc                                                                                                                                     <  1:653523/140‑1 (MQ=255)   gggCATTCTGCCCGACTTCACTGCGACAACGCAGGGTCGGCGCAGCATCGCACGCTTTTTGCTGGGCAAGCAGTTCTTCAACTGTCGTAACGCTTTTCACTTCAACCGGAATAGCCGCCAGACTGCGCTCAAGTCGGGCa                                                                                                                                    <  1:504565/140‑1 (MQ=255)             cccGACTTCACTGCGACAACGCAGGGTCGGCGCAGCATCGCACGCTTTTTGCTGGGCAAGCAGTTCTTCAACTGTCGTAACGCTTTTCACTTCAACCGGAATAGCCGCCAGACTGCGCTCAAGTCGGGCATCTCCTGGGt                                                                                                                          <  1:290177/140‑1 (MQ=255)               cGACTTCACTGCGACAACGCAGGGTCGGCGCAGCATCGCACGCTTTTTGCTGGGCAAGCAGTTCTTCAACTGTCGTAACGCTTTTCACTTCAACCGGAATAGCCGCCAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGa                                                                                                                        <  2:1014491/140‑1 (MQ=255)                     cACTGCGACAACGCAGGGTCGGCGCAGCATCGCACGCTTTTTGCTGGGCAAGCAGTTCTTCAACTGTCGTAACGCTTTTCACTTCAACCGGAATAGCCGCCAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTc                                                                                                                  >  2:11190/1‑140 (MQ=255)                       cTGCGACAACGCAGGGTCGGCGCAGCATCGCACGCGTTTTGCTGGGCAAGCAGTTCTTCTACTGTCGTAACGCTTTTCACTTCAACCGGAATTGCCGCCAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAg                                                                                                                >  2:384459/1‑140 (MQ=255)                            aCAACGCAGGGTCGGCGCAGCATCGCACGCTTTTTGCTGGGCAAGCCGTTCTTCAACTGTCGTAACTCTTTTCACTTCAACCGGAATAGCCGCCAGTCTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGAt                                                                                                           >  2:63477/1‑140 (MQ=255)                            aCAACGCAGGGTCGGCGCAGCATCGCACGCTTTTTGCTGGGCAAGCAGTTCTTCAACTGTCGTAACGCTTTTCACTTCCACCGGAATAGCCGCCCGACTGCGCTCAAGGCGGGCATCTCCTGGGTGACGTGTCAGTGGAt                                                                                                           >  2:1121092/1‑140 (MQ=255)                            aCAACGCAGGGTCGGCGCAGCATCGCACGCTTTTTGCTGGGCAAGCAGTTCTTCAACTGTCGTAACGCTTTTCACTTCAACCGGAATAGCCGCCAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGAt                                                                                                           >  1:798123/1‑140 (MQ=255)                            aCAACGCAGGGTCGGCGCAGCATCGCACGCTTTTTGCTGGGCAAGCAGTTCTTCAACTGTCGTAACGCTTTTCACTTCAACCGGAATAGCCGCCAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGAt                                                                                                           >  1:724428/1‑140 (MQ=255)                            aCAACGCAGGGTCGGCGCAGCATCGCACGCTTGTTGCTGGGCAAGCAGAACTTCAACTGTCGTAACGCTTTTCACTTCAACCGGAATAGCCGCCAGACTGCGAGCAAGTCGGGCATCTCCTGGGTCACGTTTCAggtcat                                                                                                           >  2:1121392/1‑135 (MQ=255)                            aCAACGCAGCGACGGCGCAGCATCGCACGCTTTTTGCTGGGCAAGCAGATCTTCAACTGTCGTAACGCTTTTCACTTCAACCGGAATAGCCGCCCGACTGCGCTCAAGTCGGGCATCACCTGGGTCACGTTTCAGTTGag                                                                                                           >  2:405267/1‑139 (MQ=255)                               aCGCAGGGTCGGCGCAGCATCGCACGCTTTTTGCTGGGCAAGCAGTTCTTCAACTGTCGTAACGCTTTTCACTTCAACCGGAATAGCCGCCAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCt                                                                                                        <  1:999274/140‑1 (MQ=255)                                     ggTCGGCGCAGCATCGCACGCTTTTTGCTGGGCAAGCAGTTCTTCAACTGTCGTAACGCTTTTCACTTCAACCGGAATAGCCGCCAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCa                                                                                                  <  2:161929/140‑1 (MQ=255)                                     ggTCGGCGCAGCATCGCACGCTTTTTGCTGGGCAAGCAGTTCTTCAACTGTCGTAACGCTTTTCACTTCAACCGGAATAGCCGCCAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCa                                                                                                  <  2:162010/140‑1 (MQ=255)                                          gcgcAGCATCGCACGCTTTTTGCTGGGCAAGCAGTTCTTCAACTGTCGTAACGCTTTTCACTTCAACCGGAATAGCCGCCAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGc                                                                                             >  2:902669/1‑140 (MQ=255)                                           cgcAGCATCGCACGCTTTTTGCTGGGCAAGCAGTTCTTCAACTGTCGTAACGCTTTTCACTTCAACCGGAATAGCCGCCAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGcc                                                                                            >  1:278790/1‑140 (MQ=255)                                              agcaTCGCACGCTTTTTGCTGGGCAAGCAGTTCTTCAACTGTCGTAACGCTTTTCACTTCAACCGGAATAGCCGCCAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCcgcg                                                                                         <  2:208030/140‑1 (MQ=255)                                                caTCGCACGCTTTTTGCTGGGCAAGCAGTTCTTCAACTGTCGTAACGCTTTTCACTTCAACCGGAATAGCCGCCAGACTGCGCTCTAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGACGCGcc                                                                                       >  2:461732/1‑140 (MQ=255)                                                caTCGCACGCTTTTTGCTGGGCAAGCAGTTCTTCAACTGTCGTAACGCTTTTCACTTCAACCGGAATAGCCGCCAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGcc                                                                                       <  2:1045867/140‑1 (MQ=255)                                                caTCGCACGCTTTTTGCTGGGCAAGCAGTTCTTCAACTGTCGTAACGCTTTTCACTTCAACCGGAATAGCCGCCAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGcc                                                                                       >  1:737394/1‑140 (MQ=255)                                                caTCGCACGCTTTTTGCTGGGCAAGCAGTTCTTCAACTGTCGTAACGCTTTTCACTTCAACCGGAATAGCCGCCAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGcc                                                                                       >  1:205135/1‑140 (MQ=255)                                                       cGCTTTTTGCTGGGCAAGCAGTTCTTCAACTGTCGTAACGCTTTTCACTTCAACCGGAATAGCCGCCAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGtcat                                                                                <  2:693878/140‑1 (MQ=255)                                                          tttttGCTGGGCAAGCAGTTCTTCAACTGTCGTAACGCTTTTCACTTCAACCGGAATAGCCGCCAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTcatcat                                                                             >  1:431469/1‑140 (MQ=255)                                                          tttttGCTGGGCAAGCAGTTCTTCAACTGTCGTAACGCTTTTCACTTCAACCGGAATAGCCGCCAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTcatcat                                                                             >  1:17264/1‑140 (MQ=255)                                                          tttttGCTGGGCAAGCAGTTCTTCAACTGTCGTAACGCTTTTCACTTCAACCGGAATAGCCGCCAGACTGCGCTCAAGTAGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTcatcat                                                                             >  1:82823/1‑140 (MQ=255)                                                          tttttGCTGGGCAAGCAGTTCTTCAACTGTCGTAACGCTTTTCACTTCAACCGGAATAGCCGCCAGACTGCGCTCAAGTAGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTcatcat                                                                             >  1:83938/1‑140 (MQ=255)                                                                   ggCAAGCAGTTCTTCAACTGTCGTAACGCTTTTCACTTCAACCGGAATAGCCGCCAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAAt                                                                    >  1:451212/1‑140 (MQ=255)                                                                   ggCAAGCAGTTCTTCAACTGTCGTAACGCTTTTCACTTCAACCGGAATAGCCGCCAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAAt                                                                    >  1:822883/1‑140 (MQ=255)                                                                     cAAGCAGTTCTTCAACTGTCGTAACGCTTTTCACTTCAACCGGAATAGCCGCCAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATgg                                                                  <  1:491273/140‑1 (MQ=255)                                                                        gtaGTTCTTCAACTGTCGTAACGCTTTTCACTTCAACCGGAATAGCCGCCAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGa                                                               >  1:200133/3‑140 (MQ=255)                                                                        gtaGTTCTTCAACTGTCGTAACGCTTTTCACTTCAACCGGAATAGCCGCCAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGa                                                               >  1:200110/3‑140 (MQ=255)                                                                        gtaGTTCTTCAACTGTCGTAACGCTTTTCACTTCAACCGGAATAGCCGCCAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGa                                                               >  1:200333/3‑140 (MQ=255)                                                                        gtaGTTCTTCAACTGTCGTAACGCTTTTCACTTCAACCGGAATAGCCGCCAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGa                                                               >  1:200307/3‑140 (MQ=255)                                                                                   aCTGTCGTAACGCTTTTCACTTCAACCGGAATAGCCGCCAGACTGCGCTCAAGTCGGGCATCTCCTTGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGc                                                    <  1:500870/140‑1 (MQ=255)                                                                                   aCTGTCGTAACGCTTTTCACTTCAACCGGAATAGCCGCCAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGc                                                    <  2:689682/140‑1 (MQ=255)                                                                                   aCTGTCGTAACGCTTTTCACTTCAACCGGAATAGCCGCCAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGc                                                    <  1:401152/140‑1 (MQ=255)                                                                                   aCTGTCGTAACGCTTTTCACTTCAACCGGAATAGCCGCCAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGc                                                    <  1:502198/140‑1 (MQ=255)                                                                                   aCTGTCGTAACGCTTTTCACTTCAACCGGAATAGCCGCCAGACTGCGCTCAAGTAGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGc                                                    <  2:82823/140‑1 (MQ=255)                                                                                   aCTGTCGTAACGCTTTTCACTTCAACCGGAATAGCCGCCAGACTGCGCTCAAGTAGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGc                                                    <  2:83938/140‑1 (MQ=255)                                                                                    ctctcctAACGCTATACACTACAACCGGAATAGCCGCCACACTGCGCACAAGTCCGGCATCTCCAGGGTCACCTTTCAGTTCATCCACAAGCAACACCCGCGCCCGGTCATCATCACCAAAATGGCGATAGGCTTCCGCa                                                   <  2:122248/134‑1 (MQ=255)                                                                                    cTGTCGTAACGCTTTTCACTTCAACCGGAATAGCCGCCAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCa                                                   <  1:684821/140‑1 (MQ=255)                                                                                    cTGTCGTAACGCTTTTCACTTCAACCGGAATAGCCGCCAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCa                                                   <  2:90251/140‑1 (MQ=255)                                                                                     tGTCGTAACGCTTTTCACTTCAACCGGAATAGCCGCCAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCaa                                                  <  2:1174428/140‑1 (MQ=255)                                                                                     tGTCGTAACGCTTTTCACTTCAACCGGAATAGCCGCCAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCaa                                                  <  2:1174565/140‑1 (MQ=255)                                                                                     tGTCGTAACGCTTTTCACTTCAACCGGAATAGCCGCCAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCaa                                                  <  1:1089511/140‑1 (MQ=255)                                                                                     tGTCGTAACGCTTTTCACTTCAACCGGAATAGCCGCCAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCaa                                                  <  1:445065/140‑1 (MQ=255)                                                                                      gTCGTAACGCTTTTCACTTCAACCGGAATAGCCGCCAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAg                                                 >  1:1042812/1‑140 (MQ=255)                                                                                      gTCGTAACGCTTTTCACTTCAACCGGAATAGCCGCCAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAg                                                 >  1:1042553/1‑140 (MQ=255)                                                                                      gTCGTAACGCTTTTCACTTCAACCGGAATAGCCGCCAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAg                                                 >  2:405506/1‑140 (MQ=255)                                                                                                 tttCACTTCAACCGGAATAGCCGCCAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTCa                                      >  2:714772/1‑140 (MQ=255)                                                                                                         cAACCGGAATAGCCGCCAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTCAGCGGAata                              <  2:922608/140‑1 (MQ=255)                                                                                                                aaTAGCCGCCAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTCAGCGGAATATTATCCg                       <  2:669080/140‑1 (MQ=255)                                                                                                                aaTAGCCGCCAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTCAGCGGAATATTATCCg                       <  2:881662/140‑1 (MQ=255)                                                                                                                  gcgccgccAGACCGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGAACCTCAAGCACCAGCCGCGCCCGGTCATCATGACCAACATGGCGATCGGCTTCCGCTCGGTATAAAGTCAGCGGAATATTATCCGGc                     <  2:1154403/138‑1 (MQ=255)                                                                                                                      cgccAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTCAGCGGAATATTATCCGGCACct                 <  2:577615/140‑1 (MQ=255)                                                                                                                      cgccAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTCAGCGGAATATTATCCGGCACct                 <  2:577623/140‑1 (MQ=255)                                                                                                                      cgccAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTCAGCGGAATATTATCCGGCACct                 <  2:1227820/140‑1 (MQ=255)                                                                                                                       gccAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTCAGCGGAATATTATCCGGCACctg                <  1:204557/140‑1 (MQ=255)                                                                                                                       gccAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTCAGCGGAATATTATCCGGCACctg                <  1:1232961/140‑1 (MQ=255)                                                                                                                          aGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTCAGCGGAATATTATCCGGCACctgctg             >  2:861422/1‑140 (MQ=255)                                                                                                                           gACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTCAGCGGAATATTATCCGGCACCTGCTgg            <  1:867798/140‑1 (MQ=255)                                                                                                                               gcgcTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTCAGCGGAATATTATCCGGCACCTGCTGGTGta        <  2:22023/140‑1 (MQ=255)                                                                                                                                     aaGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTCAGCGGAATATTATCCGGCACCTGCTGGTGTATATAtt  <  1:1069948/140‑1 (MQ=255)                                                                                                                                     aaGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTCAGCGGAATATTATCCGGCACCTGCTGGTGTATATAtt  <  1:1069615/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                        CAGGGCATTCTGCCCGACTTCACTGCGACAACGCAGGGTCGGCGCAGCATCGCACGCTTTTTGCTGGGCAAGCAGTTCTTCAACTGTCGTAACGCTTTTCACTTCAACCGGAATAGCCGCCAGACTGCGCTCAAGTCGGGCATCTCCTGGGTGACGTTTCAGTTGATCCTCAAGCAACAGCCGCGCCCGGTCATCATGACCAAAATGGCGATAGGCTTCCGCAAGGTATAAAGTCAGCGGAATATTATCCGGCACCTGCTGGTGTATATATT  >  NC_000913/590472‑590743

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