Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,914,184	G→A	F211F (TTC→TTT)	prc ←	tail‑specific protease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,914,184	0	G	A	98.2%	186.5 / ‑2.8	54	coding (633/2049 nt)	prc	tail‑specific protease
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (25/28);  minor base . (1/0);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.33e-01

ACCATCGGCTTGCCTGTTTGACCAACACCGACAATTTTGTCACCAACGCTGATAGCTTTACTCTTCGCTGCCGGACCACCTGCCACCATCGAATTGATAACGGTGTAGTCATCATCCATTTGCAGCACTGCGCCAATACCTTCCAGCGACAAACTCATTTCAGTGTTGAACTGTTCGGTATTACGCGGGGAAAGATAGTTGGTATGCGGGTCGATTTCACGCGCAAACGCCGTCATTGCCAGCGAGAAAACA‑TCTTCGCTGTTGGTTTGCGCCAGACGACGAATGGCAAATTTGTAGCGGCGAGTCAGGGTTTCACGAATTTCTTTATCCGTTTTTCCTGTCAGCTTCAGGCTTAACTCGTCGAATTTGACTTTACTGTCCCACAGCGCGTTCAACTCAGCCTCGTTTTTCGGCCAGGGCGCTTTGCTGCGGTCAAGGTTATAAGTGTCGTTGCCGGTGAAATCCATCGGCTTTTCCAGTACCGA  >  NC_000913/1913939‑1914423                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 aCCATCGGCTTTCCTGTTTGACCAACACCGACAATTTTGTCACCAACGCTGATAGCTTTACTCTTCGCTGCCGGACCACCTGCCACCATCGAATTGATAACGGTGTAGTCATCATCCATTTGCAGCACTGCGCCAATACCTTCCAGCGACAAACTCATTTCAGTGTTGAACTGTTCGGTATTACGCGGGGAAAGATAGTTGGTATGCGGGTCGATTTCACGCGCAAACGCCGTCATTGCCAGCGAAAAAAc                                                                                                                                                                                                                                             <  1:605581/251‑1 (MQ=255)         cTTGCCTGTTTGACCAACACCGACAATTTTGTCACCAACGCTGATAGCTTTACTCTTCGCTGCCGGACCACCTGCCACCATCGAATTGATAACGGTGTAGTCATCATCCATTTGCAGCACTGCGCCAATACCTTCCAGCGACAAACTCATTTCAGTGTTGAACTGTTCGGTATTACGCGGGGAAAGATAGTTGGTATGCGGGTCGATTTCACGCGCAAACGCCGTCATTGCCAGCGAAAAAACA‑TCTTCGc                                                                                                                                                                                                                                    >  1:17479/1‑251 (MQ=255)            gCCTGTTTGACCAACACCGACAATTTTGTCACCAACGCTGATAGCTTTACTCTTCGCTGCCGGACCACCTGCCACCATCGAATTGATAACGGTGTAGTCATCATCCATTTGCAGCACTGCGCCAATACCTTCCAGCGACAAACTCATTTCAGTGTTGAACTGTTCGGTATTACGCGGGGAAAGATAGTTGGTATGCGGGTCGATTTCACGCGCAAACGCCGTCATTGCCAGCGAAAAAACA‑TCTTCGCTg                                                                                                                                                                                                                                  <  2:328806/250‑1 (MQ=255)                   ggACCAACTCCGACAAATTTGTCACCAACGCTGATAGCTTTACTCTTCGCTGCCGGACCTCCTGCCACCATCGAATTGAGAACGGTGTAGTCATCATCCATTTGCATCACTTCGCCAAGACCTTCCATCGACAAACTCATTTCAGTTTTGAACTTTTCGGTATTACGCAGGGAAATATAGTTGGTATGCGTGTGGATTTCACGCTCAAACGCCGTCATTGCCAGCGAAAAAACA‑TCTTCGCTGTTGGttt                                                                                                                                                                                                                           <  1:600114/249‑1 (MQ=255)                                   ttttGTCACCAACGCTGATAGCTTTACTCTTCGCTGCCGGACCACCTGCCACCATCGAATTGATAACGGTGTAGTCATCATCCATTTGCAGCACTGCGCCAATACCTTCCAGCGACAAACTCATTTCAGTGTTGAACTGTTCGGTATTACGCGGGGAAAGATAGTTGGTATGCGGGTCGATTTCACGCGCAAACGCCGTCATTGCCAGCGAAAAAACA‑TCTTCGCTGTTGGTTTGCGCCAGACGACGAATg                                                                                                                                                                                                          >  2:53764/1‑251 (MQ=255)                                        tCACCAACGCTGATAGCTTTACTCTTCGCTGCCGGACCACCTGCCACCATCGAATTGATAACCGTGTAGTCATCATCCATTTTCAGCACTGCGCCAATACCTTCCAGCGACAAACTCATTGCAGTGTTGAACTGTTCGGTATTACGCGGGGAAAGATAGTTGGTATGCGGGTCGATTTCACGCGCAAACGCCGTCATTGCCAGCGAAAAAACA‑TCTTCGCTGTTGGTTTGCGCCAGACGACGAATGGCaa                                                                                                                                                                                                      <  1:453960/250‑1 (MQ=255)                                          aCCAACGCTGATAGCTTTACTCTTCGCTGCCGGACCACCTGCCACCATCGAATTGATAACGGTGTAGTCATCATCCAGTTGCAGCACTGCGCCAATACCTTCCAGCGACAAACTCATTTCAGTGTTGAACTGTTCGGTATTACGCGGGGAAAGATAGTTGGTATGCGGGTCGATTTCACGCGCAAACGCCGTCATTGCCAGCGAAAAAACA‑TCTTCGCTGTTGGTTTGCGCCAGACGACGAATGGCAAAtt                                                                                                                                                                                                   <  2:441854/251‑1 (MQ=255)                                             aaCGCTGATAGCTTTACTCTTCGCTGCCGGACCACCTGCCACCATCGAATTGATAACGGTGTAGTCATCATCCATTTGCAGCACTGCGCCAATACCTTCCAGCGACAAACTCATTTCAGTGTTGAACTGTTCGGTATTACGCGGGGAAAGATAGTTGGTATGCGGGTCGATTTCACGCGCAAACGCCGTCATTGCCAGCGAAAAAACA‑TCTTCGCTGTTGGTTTGCGCCAGACGACGAATGGCAAATTTGt                                                                                                                                                                                                <  1:331701/251‑1 (MQ=255)                                              aCGCTGATAGCTTTACTCTTCGCTGCCGGACCACCTGCCACCATCGAATTGATAACGGTGTAGTCATCATCCATTTGCAGCACTGCGCCAATACCTTCCAGCGACAAACTCATTTCAGTGTTGAACTGTTCGGTATTACGCGGGGAAAGATAGTTGGTATGCGGGTCGATTTCACGCGCAAACGCCGTCATTGCCAGCGAAAAAACA‑TCTTCGCTGTTGGTTTGCGCCAGACGACGAATGGCAAATTTGTa                                                                                                                                                                                               <  1:73776/251‑1 (MQ=255)                                                gCTGATAGCTTTACTCTTCGCTGCCGGACCACCTGCCACCATCGAATTGATAACGGTGTAGTCATCATCCATTTGCAGCACTGCGCCAATACCTTCCAGCGACAAACTCATTTCAGTGTTGAACTGTTCGGTATTACGCGGGGAAAGATAGTTGGTATGCGGGTCGATTTCACGCGCAAACGCCGTCATTGCCAGCGAAAAAACA‑TCTTCGCTGTTGGTTTGCGCCAGACGACGAATGGCAAATTTGTAgc                                                                                                                                                                                             >  2:165106/1‑251 (MQ=255)                                                      aGCTTTACTCTTCGCTGCCGGACCACCTGCCACCATCGAATTGATAACGGTGTAGTCATCATCCATTTGCAGCACTGCGCCAATACCTTCCAGCGACAAACTCATTTCAGTGTTGAACTGTTCGGTATTACGCGGGGAAAGATAGTTGGTATGCGGGTCGATTTCACGCGCAAACGCCGTCATTGCCAGCGAAAAAACA‑TCTTCGCTGTTGGTTTGCGCCAGACGACGAATGGCAAATTTGTAGCGGCGa                                                                                                                                                                                        <  2:246049/250‑1 (MQ=255)                                                      aGCTTTACTCTTCGCTGCCGGACCACCTGCCACCATCGAATTGATAACGGTGTAGTCATCATCCATTTGCAGCACTGCGCCAATACCTTCCAGCGACAAACTCATTTCAGTGTTGAACTGTTCGGTATTACGCGGGGAAAGATAGTTGGTATGCGGGTCGATTTCACGCGCAAACGCCGTCATTGCCAGCGAAAAAACA‑TCTTCGCTGTTGGTTTGCGCCAGACGACGAAGGGCAAATTTGTAGCGGCGAg                                                                                                                                                                                       <  2:51631/251‑1 (MQ=255)                                                            aCTCTTCGCTGCCGGACCACCTGCCACCATCGAATTGATAACGGTGTAGTCATCATCCATTTGCAGCACTGCGCCAATACCTTCCAGCGACAAACTCATTTCAGTGTTGAACTGTTCGGTATTACGCGGGGAAAGATAGTTGGTATGCGGGTCGATTTCACGCGCAAACGCCGTCATTGCCAGCGAAAAAACA‑TCTTCGCTGTTGGTTTGCGCCAGACGACGAATGGCAAATTTGTAGCGGCGAGTCAggg                                                                                                                                                                                 >  2:151264/1‑251 (MQ=255)                                                                 tcggtGCCGGACCACCTGCCTCCATCGAATTGATAACGGTGTAGTCATCATCCATTTTCAGCACTGCGCCAATACCTTCCAGCGACAAACTCATTTCAGTGTTGAACTGTTCGGTATTACGCGGGGAAATATAGTTGGTATGCGGGTCGATTTCACGCGCAAACGCCGTCATTGCCAGCGAAAAAACA‑TCTTCGCTGTTGTTTTGCGCCAGACGACTAATGACAAATTTGTAGCTGCGAGTCAGGGttt                                                                                                                                                                              <  2:129975/245‑1 (MQ=255)                                                                  cGCTGCCGGACCACCTGCCACCATCGAATTGATAACGGTGTAGTCATCATCCATTTGCAGCACTGCGCCAATACCTTCCAGCGACAAACTCATTTCAGTGTTGAACTGTTCGGTATTACGCGGGGAAAGATAGTTGGTATGCGGGTCGATTTCACGCGCAAACGCCGTCATTGCCAGCGAAAAAACA‑TCTTCGCTGTTGGTTTGCGCCAGACGACGAATGGCAAATTTGTAGCGGCGAGTCAGGGTTTCAc                                                                                                                                                                           <  1:54247/251‑1 (MQ=255)                                                                            ccaacTGCCACCAGCGAATTGATAACGGTGTAGTCATCATCCATTTGCAGCACTGCGCCAATACCTTCCATCGACAAACTCATTTCAGTGTTGAACTGTTCGGTATTACGCGGGGAAAGATAGTTGGTATGCGGGTCGATTTCACGCGCAAACGCCGTCATTGCCAGCGAAAAAACA‑TCTTCGCTGTTGGTTTGCGCCAGACGACGAATGGCAAATTTGTAGCGGCGAGTCAGGGTTTCACGAATTTCttt                                                                                                                                                                 <  2:233071/247‑1 (MQ=255)                                                                                   ccaccaTCGAATTGATAACGGTGTAGTCATCATCCATTTGCAGCACTGCGCCAATACCTTCCAGCGACAAACTCATTTCAGTGTTGAACTGTTCGGTATTACGCGGGGAAAGATAGTTGGTATGCGGGTCGATTTCACGCGCAAACGCCGTCATTGCCAGCGAAAAAACA‑TCTTCGCTGTTGGTTTGCGCCAGACGACGAATGGCAAATTTGTAGCGGCGAGTCAGGGTTTCACGAATTTCTTTATCCGtt                                                                                                                                                          >  1:148194/1‑251 (MQ=255)                                                                                    caccaTCGAATTGATAACGGTGTAGTCATCATCCATTTGCAGCACTGCGCCAATACCTTCCAGCGACAAACTCATTTCAGTGTTGAACTGTTCGGTATTACGCGGGGAAAGATAGTTGGTATGCGGGTCGATTTCACGCGCAAACGCCGTCATTGCCAGCGAAAAAACA‑TCTTCGCTGTTGGTTTGCGCCAGACGACGAATGGCAAATTTGTAGCGGCGAGTCAGGGTTTCACGAATTTCTTTATCCGttt                                                                                                                                                         <  1:602982/251‑1 (MQ=255)                                                                                      ccaTCGAATTGATAACGGTGTAGTCATCATCCATTTGCAGCACTGCGCCAATACCTTCCAGCGACAAACTCATTTCAGTGTTGAACTGTTCGGTATTACGCGGGGAAAGATAGTTGGTATGCGGGTCGATTTCACGCGCAAACGCCGTCATTGCCAGCGAAAAAACA‑TCTTCGCTGTTGGTTTGCGCCAGACGACGAATGGCAAATTTGTAGCGGCGAGTCAGGGTTTCACGAATTTCTTTATCCGttttt                                                                                                                                                       >  2:234549/1‑251 (MQ=255)                                                                                                       gtgtAGTCATCATCCATTTGCAGCACTGCGCCAATACCTTCCAGCGACAAACTCATTTCAGTGTTGAACTGTTCGGTATTACGCGGGGAAAGATAGTTGGTATGCGGGTCGATTTCACGCGCAAACGCCGTCATTGCCAGCGAAAAAACA‑TCTTCGCTGTTGGTTTGCGCCAGACGACGAATGGCAAATTTGTAGCGGCGAGTCAGGGTTTCACGAATTTCTTTATCCGTTTTTCCTGTCAGCTTCAGGCt                                                                                                                                      <  2:347660/251‑1 (MQ=255)                                                                                                        tgtAGTCATCATCCATTTGCAGCACTGCGCCAATGCCTTCCAGCGACAAACTCATTTCAGGGTTGAACTGGTCGGTAGTACGCGGGGAAAGATAGTTGGTATGCGGGTCGATTTCACGCGCAAACGCCGTCATTGCCAGCGAAAAAACA‑TCTTCGCTGTTGGTTTGCGCCAGACGACGAATGGCAAATTTGTAGCGGCGAGTCAGGGTTTCACGAATTTCGTTATCCGTTTTTCCTGTCAGCTTCAGGCt                                                                                                                                      <  2:573080/250‑1 (MQ=255)                                                                                                                    ccATTTGCAGCACTGCGCCAATACCTTCCAGCGACAAACTCATTTCAGTGTTGAACTGTTCGGTATTACGCGGGGAAAGATAGTTGGTATGCGGGTCGATTTCACGCGCAAACGCCGTCATTGCCAGCGAAAAAACA‑TCTTCGCTGTTGGTTTGCGCCAGACGACGAATGGCAAATTTGTAGCGGCGAGTCAGGGTTTCACGAATTTCTTTATCCGTTTTTCCTGTCAGCTTCAGGCTTAACTCGTCGAAt                                                                                                                         <  1:17675/251‑1 (MQ=255)                                                                                                                          gcagcaCTGCGCCAATACCTTCCAGCGACAAACTCATTTCAGCGTTGAACTGTTCGGAATTACGCGGGGAAAGATAGTTGGTATGCGGGTCGATTTCACGCGCAAACGCCGTCATTGCCAGCGAAAAAACA‑TATTCGCTGTTGGTTTGCGCCAGACGACGAATGGCAAATTTGTAGCGGCGAGTCAGGGTTTTACGAATTTCTTTATTCGTTTTTCCTTTAAGCTTCAGGCTTAACTCGGCGAATTTGACt                                                                                                                   >  2:166029/1‑251 (MQ=255)                                                                                                                              caCTGCGCCAATACCTTCCAGCGACAAACTCATTTCAGTGTTGAACTGTTCGGTATTACGCGGGGAAAGATAGTTGGTATGCGGGTCGATTTCACGCGCAAACGCCGTCATTGCCAGCGAAAAAACA‑TCTTCGCTGTTGGTTTGCGCCAGACGACGAATGGCAAATTTGTAGCGGCGAGTCAGGGTTTCACGAATTTCTTTATCCGTTTTTCCTGTCAGCTTCAGGCTTAACTCGTCGAATTTGACTTTAc                                                                                                               >  2:113211/1‑251 (MQ=255)                                                                                                                                  gcgcCAATACCTTCCAGCGACAAACTCATTTCAGTGTTGAACTGTTCGGTATTACGCGGGGAAAGATAGTTGGTATGCGGGTCGATTTCACGCGCAAACGCCGTCATTGCCAGCGAAAAAACA‑TCTTCGCTGTTGGTTTGCGCCAGACGACGAATGGCAAATTTGTAGCGGCGAGTCAGGGTTTCACGAATTTCTTTATCCGTTTTTCCTGTCAGCTTCAGGCTTAACTCGTCGAATTTGACTTTACTGt                                                                                                            <  1:308418/250‑1 (MQ=255)                                                                                                                                       aaTACCTTCCAGCGACAAACTCATTTCAGTGTTGAACTGTTCGGTATTACGCGGGGAAAGATAGTTGGTATGCGGGTCGATTTCACGCGCAAACGCCGTCATTGCCAGCGAAAAAACA‑TCTTCGCTGTTGGTTTGCGCCAGACGACGAATGGCAAATTTGTAGCGGCGAGTCAGGGTTTCACGAATTTCTTTATCCGTTTTTCCTGTCAGCTTCAGGCTTAACTCGTCGAATTTGACTTTACTGTCCcaca                                                                                                      <  2:502347/251‑1 (MQ=255)                                                                                                                                            cTTCCAGCGACAAACTCATTTCAGTGTTGAACTGTTCGGTATTACGCGGGGAAAGATAGTTGGTATGCGGGTCGATTTCACGCGCAAACGCCGTCATTGCCAGCGAAAAAACA‑TCTTCGCTGTTGGTTTGAGCCAGACGACGAATGGCAAATTTGTAGCGGCGAGTCAGGGTTTCACGAATTTCTTTATCCGTTTTTCCTTTCAGCTTCAGGCTTAACTCGTCGAATTTGACTTTACTGTCCCACAGcgcg                                                                                                 >  1:631548/1‑251 (MQ=255)                                                                                                                                               ccAGCGACAAACTCATTTCAGTGTTGAACTGTTCGGTCTTACGCGGGGAAAGATAGTTGGTATGCGGGTCGATTTCACGCGCAAACGCCGTCATTGCCAGCGAAAAAACA‑TCTTCGCTGTTGGTTTGCGCCAGACGCCGAATGGCAAATTTGTAGCGGCGAGTCAGGGCTTCACGAATTTCTTTATCCGTTTTTCCTGTTAGCTTAAGGCTTAACTCGCCGAATTTGAATTTACTTTCCCACAGCGCGTTc                                                                                              >  1:309456/1‑251 (MQ=255)                                                                                                                                                  gCGACAAACTCATTTAAGTGTTGAACTGTTCGGTATTACGCGGGGAAAGATAGTTGGTATGCGGGTCGATTTCACGCGCAAACGCCGTCATTGCCAGCGAACAACCA‑TCTTCGCTGTTGGTTTGCGCCAGACGACGAATGGCAAATTTGTAGCGGCGAGTCAGGGTTTCACGAATTTCTTTATCCGTTTTTCCTGTCAGCTTCAGGCTTAACTCGTCGAATTTGACTTTACTGTCCCCCAGCGCGTTCAAc                                                                                           >  1:537634/1‑251 (MQ=255)                                                                                                                                                          cTCATTTCAGTGTTGAACTGTTCGGTATTACGCGGGGAAAGATAGTTGGTATGCGGGTCGATTTCACGCGCAAACGCCGTCATTGCCAGCGAAAAAACA‑TCTTCGCTGTTGGTTTGCGCCAGACGACGAATGGCAAATTTGTAGCGGCGAGTTAGGGTTTCACGAACTTCTTTATCCGTTTTTCCTGTCAGCTTCAGGCTTAACTCGTCGAATTTTACTTTACTGTCCCACCGCGCGTTCAACTCCGCCTc                                                                                   >  2:69404/1‑251 (MQ=255)                                                                                                                                                               ttCAGTGTTGAACTGTTCGGTATTACGCGGGGAAAGATAGTTGGTATGCGGGTCGATTTCACGCGCAAACGCCGTCATTGCCAGCGAAAAAACA‑TCTTCGCTGTTGGTTTGCGCCAGACGACGAATGGCAAATTTGTAGCGGCGAGTCAGGGTTTCACGAATTTCTTTATCCGTTTTTCCTGTCAGCTTCAGGCTTAACTCGTCGAATTTGACTTTACTGTCCCACAGCGCGTTCAACTCAGCCTCGttt                                                                               <  2:117468/250‑1 (MQ=255)                                                                                                                                                                 cAGTGTTGAACTGTTCGGTATTACGCGGGGAAAGATAGTTGGTATGCGGGTCGATTTCACGCGCAAACGCCGTCATTGCCAGCGAAAAAACA‑TCTTCGCTGTTGGTTTGCGCCAGACGACGAATGGCAAATTTGTAGCGGCGAGTCAGGGTTTCACGAATTTCTTTATCCGTTTTTCCTGTCAGCTTCAGGCTTAACTCGTCGAATTTGACTTTACTGTCCCACAGCGCGTTCAACTCAGCCTCGTTTTTc                                                                            <  2:234777/251‑1 (MQ=255)                                                                                                                                                                       tGAACTGTTCGGTATTACGCGGGGCAAGATAGTTGGTATGCGGGTCGATTTCACGCGCAAACGCCGTCATTGCCAGCGA‑AAAACAATCTTCGCAGTTGGTTTGCGCCAGACGACGAATGGCAAATTTGTAGCGGCGAGTCAGGGTTTCACGAATTTCTTTCTCCGTTTTTCCTGTCAGCTTCAGGCTTAACTCGTCGAATTTGACTTTACTGTCCCACAACTCGTTCAACTCCGCCTAGTTTTTCGGca                                                                        >  2:40416/1‑248 (MQ=255)                                                                                                                                                                       tGAACTGTTCGGTATTACGCGGGGAAAGATAGTTGGTATGCGGGTCGATTTCACGCGCAAACGCCGTCATTGCCAGCGAAAAAACA‑TCTTCGCTGTTGGTTTGCGCCAGACGACGAATGGCAAATTTGTAGCGGCGAGTCAGGGTTTCACGAATTTCTTTATCCGTTTTTCCTGTCAGCTTCAGGCTTAACTCGTCGAATTTGACTTTACTGTCCCACAGCGCGTTAAACTCAGCCTCGTTTTTCGGcc                                                                        >  1:423276/1‑249 (MQ=255)                                                                                                                                                                          aCTGTTCGGTATTACGCGGGGAAATATAGTTTGTATGCGGGTCGATTTCACGCGCAAACGCCGTCATTGCCAGCGAAAAAACA‑TCTTCGCTGTTGGTTTGCGCCAGACGACGAATGGCAAATTTGTAGCGGCGAGTCAGGTTTTCACGAATTTCTTTATCCGTTTTTCCTGTCAGCTTCAGGCTTAACTCTTCGAATTTGTCTTTACTGTCCCACAGCGCGTTCAACTCAGCCTCGTTTTTCGGCCAGGgc                                                                   <  1:617627/251‑1 (MQ=255)                                                                                                                                                                                 ggTATTACGCGGGGAAAGATAGTTGGTATGCGGGTCGATTTCACGCGCAAACGCCGTCATTGCCAGCGAAAAAACA‑TCTTCGCTGTTGGTTTGCGCCAGACGACGAATGGCAAATTTGTAGCGGCGAGTCAGGGTTTCACGAATTTCTTTATCCGTTTTTCCTGTCAGCTTCAGGCTTAACTCGTCGAATTTGACTTTACTGTCCCACAGCGCGTTCAACTCAGCCTCGTTTTTCGGCCAGGGCGCTTtgc                                                            >  1:5861/1‑251 (MQ=255)                                                                                                                                                                                 ggTATTACGCGGGGAAAGATAGTTGGTATGCGGGTCGATTTCACGCGCAAACGCCGTCATTGCCAGCGAAAAAACA‑TCTTCGCTGTTGGTTTGCGCCAGACGACGAATGGCAAATTTGTAGCGGCGAGTCAGGGGTTCACGAATTTCTTTATCCGTTTTTCCTGTCAGCTTCAGGCTTAACTCGTCGAATTTGACTTTACTGTCCCACAGCGCGTTCAACTCAGCCTCGTTTTTCGGCCAGGGCGCTTtgc                                                            >  1:403996/1‑251 (MQ=255)                                                                                                                                                                                  gTATTACGCGGGGAAAGATAGTTGGTATGCGGGTCGATTTCACGCGCAAACGCCGTCATTGCCAGCGAAAAAACA‑TCTTCGCTGTTGGTTTGCGCCAGACGACGAATGGCAAATTTGTAGCGGCGAGTCAGGGTTTCACGAATTTCTTTATCCGTTTTTCCTGTCAGCTTCAGGCTTAACTCGTCGAATTTGACTTTACTGTCCCACAGCGCGTTCAACTCAGCCTCGTTTTTCGGCCAGGGCGCTTtgct                                                           >  1:348122/1‑251 (MQ=255)                                                                                                                                                                                   tATTACGCGGGGAAAGATAGTTGGTATGCGGGTCGATTTCACGCGCAAACGCCGTCATTGCCAGCGAAAAAACA‑TCTTCGCTGTTGGTTTGCGCCAGACGACGAATGGCAAATTTGTAGCGGCGAGTCAGGGTTTCACGAATTTCTTTTTCCGTTTTTCCTGTCCGCTTCAGGCTTAACTCGTCGCATTTGACTTTAc                                                                                                               >  2:561734/1‑198 (MQ=255)                                                                                                                                                                                   tATTACGCGGGGAAAGATAGTTGGTATGCGGGTCGATTTCACGCGCAAACGCCGTCATTGCCAGCGAAAAAACA‑TCTTCGCTGTTGGTTTGCGCCAGACGACGAATGGCAAATTTGTAGCGGCGAGTCAGGGTTTCACGAATTTCTTTATCCGTTTTTCCTGTCAGCTTCAGGCTTAACTCGTCGAATTTGACTTTAc                                                                                                               <  1:561734/198‑1 (MQ=255)                                                                                                                                                                                   tATTACGCGGGGAAAGATAGTTGGTATGCGGGTCGATTTCACGCGCAAACGCCGTCATTGCCAGCGAAAAAACA‑TCTTCGCTGTTGGTTTGCGCCAGACGACGAATGGCAAATTTGTAGCGGCGAGTCAGGGTTTCACGAATTTCTTTATCCGTTTTTCCTGTCAGCTTCAGGCTTAACTAGTCGAATTTGACTTTACTTTCCCACAGCGCGTTCAACCAACCCTCGTTTTTCGGCCAGGGCGCttttctt                                                          >  2:546892/1‑247 (MQ=255)                                                                                                                                                                                           ggggAAAGATAGTTGGTCTGCGGGTCGATTTCACGCGCAAACGCCGCCATTGCCAGCGAAAAAACA‑TCTTCGCTGTTGGTTTGCGCCAGACGACGAATGGCAAATTTGTAGCGGCGAGTCAGGGGTTCACGAATTTCTTTATCCGTTTTTCCTGTCAGCTTCAGGCTTAACTCGGCCAATTTGACTTTACTGTCCCACAGCGCGTTCAAATCAGCCTCGTTTTTCGGCCAGGGTCCTTTGCTGCGGTCa                                                    >  2:89588/1‑249 (MQ=255)                                                                                                                                                                                                     aGTTGGTATGCGGGTCGATTTCACGCGCAAACGCCGTCATTGCCAGCGAAAAAACA‑TCTTCGCTGTTGGTTTGCGCCAGACGACGAATGGCAAATTTGTAGCGGCGAGTCAGGGTTTCACGAATTTCTTTATCCGTTTTTCCTGTCAGCTTCAGGCTTAACTCGTCGAATTTGACTTTACTGTCCCACAGCGCGTTCAACTCAGCCTCGTTTTTCGGCCAGGGCGCTTTGCTGCGGTCAAGGTTATAAgtg                                        <  1:598097/251‑1 (MQ=255)                                                                                                                                                                                                     aGTTGGTATGCGGGTCGATTTCACGCGCAAACGCCGTCATTGCCAGCGAAAAAACA‑TCTTCGCTGTTGGTTTGCGCCAGACGACGAATGGCAAATTTGTAGCGGCGAGTCAGGGTTTCACGAATTTCTTTATCCGTTTTTCCTGTCAGCTTCAGGCTTAACTCGTCGAATTTGACTTTACTGTCCCACAGCGCGTTCAACTCAGCCTCGTTTTTCGGCCAGGGCGCTTTGCTGCGGTCAAGGTTATAAgtg                                        <  1:165106/251‑1 (MQ=255)                                                                                                                                                                                                     aGTTGGTATGAGGGCCGATTTCACGCGCAAACGCCGTCCTTGCCAGCGAAAAAACA‑TCTTCGCTGTTGGTTTGCGCCAGACGACGAATGGCAAATTTGAAGCGGCGAGTCCGGGGTTCACGCATTTCTTTATCCGTTTTTCCTGTCGGCTTCAGGCTTAACTCGTCGAATTTGACTTTACTGTCCCACAGCGCGTTCAACTCAGCCTCGTTTTTAGGCCAGGGCGCTTTGCTGCGGGCAAGGGTATANgtg                                        >  2:589616/1‑251 (MQ=255)                                                                                                                                                                                                      gTTGGTATGCGGGTCGATTTCACGCGCAAACGCCGTCATTGCCAGCGAAAAAACA‑TCTTCGCTGGTGGTTTGCGCCAGACGACGAATGGCAAATTTGTAGCGGCGAGTCAGGGTTTCACGAATTTATTTATCCGGTTTTCCTGTCAGCTTCAGGCTTAACTCCACGAATTTGACTTTACTGTCCCACAGCGAGTTCAACTCAGCCTCGTTTTTCGGCCAGGGCGCTTTGCTGCGGGCAAGGTTATAAgtgt                                       >  2:632780/1‑251 (MQ=255)                                                                                                                                                                                                                 ggTCGATTTCACGCGCAAACGCCGTCATTGCCAGCGAAAAAACA‑TCTTCGCTGTTGGTTTGCGCCAGACGACGAATGGCAAATTTGTAGCGGCGAGTCAGGGTTTCACGAATTTCTTTATCCGTTTTTCCTGTCAGCTTCAGGCTTAACTCGTCGAATTTGACTTTACTGTCCCACAGCGCGTGCAACTCAGCCTCGTTTTTCGGCCAGGGCGCTTTGCTGCTGTCAAGGTTATAAGTGTCGTTGCCGGTg                            <  1:478489/251‑1 (MQ=255)                                                                                                                                                                                                                   tCGATTTCACGCGCAAACGCCGTCATTGCCAGCGAAAAAACA‑TCTTCGCTGTTGGTTTGCGCCAGACGACGAATGGCAAATTTGTAGCGGCGAGTCAGGGGTTCACGAATTTCTTTATCCGTTTTTCCTGTCAGCTTCAGGCTTAACTCGTCGAATTTGACTTTACTGTCCCACAGCGCGTTCAACTCAGCCTCGTTTTTCGGCCAGGGCGCTTTGCTGCGGTCAAGGTTATAAGTGTCGTTGCCGGTGa                           >  2:116550/1‑250 (MQ=255)                                                                                                                                                                                                                     gATTTCACGCGCAAACGCCGTCATTGCCAGCGAAAAAACA‑TCTTCGCTGTTGGTTTGCGCCAGACGACGAATGGCAAATTTGTAGCGGCGAGTCAGGGTTTCACGAATTTCTTTATCCGTTTTTCCTGTCAGCTTCAGGCTTAACTCGTCGAATTTGACTTTACTGTCCCACAGCGCGTTCAACTCAGCCTCGTTTTTCGGCCAGGGCGCTTTGCTGCGGTAAAGGTTATAAGTGTAGTTGCCGGGGAAAt                        >  2:637355/1‑251 (MQ=255)                                                                                                                                                                                                                        ttCACGCGCAAACGCCGTCATTGCCAGCGAAAAAACA‑TCTTCGCTGTTGGTTTGCGCCAGACGACGAATGGCAAATTTGTAGCGGCGAGTCAGGGTTTCACGAATTTCTTTAGCCGTTTTTCCTGTCAGCTTCAGGCTTGACTCGTCGAATTTGACTTTACTGTCCCACAGCGCGTTCAACTCAGCCTCGTTTTTCGGCCAGGGCGCTTTGCTGCGGTCAAGGTTATAAGTGTCGTTGCCGGTGAAATcc                      <  2:174423/250‑1 (MQ=255)                                                                                                                                                                                                                                 aaaCGCCGTCATTGCCAGCGAAAAAACA‑TCTTCGCTGTTGGTTTGCGCCAGACGACGAATGGCAAATTTGTAGCGGCGAGTCAGGGTTTCACGAATTTCTTTATCCGTTTTTCCTGTCAGCTTCAGGCTTAACTCGTCGAATTTGACTTTAATGTCCCACAGCGAGTTCAACTCAGCCTCGTTTTTCGGCCAGGGCGCTTTGCTGCGGTCAAGGTTATAAGGGTCGTTGCCGGGGAAATCCATCGGCtttt            >  1:452069/1‑251 (MQ=255)                                                                                                                                                                                                                                     gCCGTCATTGCCAGCGAAAAAACA‑TCTTCGCTGTTGGTTTGCGCCAGACGACGAATGGCAAATTTGTAGCGGCGAGTCAGGGTTTCACGAATTTCTTTATCCGTTTTTCCTGTCAGCTTCAGGCTTAACTCGTCGAATTTGACTTTACTGTCCCACAgc                                                                                                    <  2:314125/159‑1 (MQ=255)                                                                                                                                                                                                                                     gCCGTCATTGCCAGCGAAAAAACA‑TCTTCGCTGTTGGTTTGCGCCAGACGACGAATGGCAAATTTGTAGCGGCGAGTCAGGGTTTCACGAATTTCTTTATCCGTTTTTCCTGTCAGCTTCAGGCTTAACTCGTCGAATTTGACTTTACTGTCCCACAgc                                                                                                    >  1:314125/1‑159 (MQ=255)                                                                                                                                                                                                                                           aTTGCCAGCGAAAAAACA‑TCTTCGCTGTTGGTTTGCGCCAGACGACGAATGGCAAATTTGTAGCGGCGAGTCAGGGTTTCACGAATTTCTTTATCCGTTTTTCCTGTCAGCTTCAGGCTTAACTCGTCGAATTTGACTTTACTGTCCCACAGCGCGTTCAACTCAGCCTCGTTTTTCGGCCAGGGCGCTTTGCTGCGGTCAAGGTTATAAGTGTCGTTGCCGGTGAAATCCATCGGCTTTTCCAGTACCGa  <  2:348122/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 ACCATCGGCTTGCCTGTTTGACCAACACCGACAATTTTGTCACCAACGCTGATAGCTTTACTCTTCGCTGCCGGACCACCTGCCACCATCGAATTGATAACGGTGTAGTCATCATCCATTTGCAGCACTGCGCCAATACCTTCCAGCGACAAACTCATTTCAGTGTTGAACTGTTCGGTATTACGCGGGGAAAGATAGTTGGTATGCGGGTCGATTTCACGCGCAAACGCCGTCATTGCCAGCGAGAAAACA‑TCTTCGCTGTTGGTTTGCGCCAGACGACGAATGGCAAATTTGTAGCGGCGAGTCAGGGTTTCACGAATTTCTTTATCCGTTTTTCCTGTCAGCTTCAGGCTTAACTCGTCGAATTTGACTTTACTGTCCCACAGCGCGTTCAACTCAGCCTCGTTTTTCGGCCAGGGCGCTTTGCTGCGGTCAAGGTTATAAGTGTCGTTGCCGGTGAAATCCATCGGCTTTTCCAGTACCGA  >  NC_000913/1913939‑1914423

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