Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	664,907	G→A	T95M (ACG→ATG)	rlpA ←	rare lipoprotein RlpA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	664,907	0	G	A	98.0%	149.3 / ‑5.9	51	T95M (ACG→ATG)	rlpA	rare lipoprotein RlpA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (28/21);  minor base C (0/1);  total (28/23)
Fisher's exact test for biased strand distribution p-value = 4.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

AACCATCCTGGGCAACAATAATCGGATCGATACGAACTTTGGTGTTGTTTGACGTGTTAAGACGGTCAGCTGCCGCGCGAGAAAGTGAAATAACGCGGTCGTTGCCGTAAGGGCCGCGATCATTAATGCGCACCACGATCATTCGCCCGTTAGCCAGGTTAGTG‑ATTCTGGCGTAGCTGGGGATCGGAAGCGTTGGATGGGCTGCCGTCAGCTGTGTCGGATCGAATGCTTCGCCAGAGGCCGTCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTAC  >  NC_000913/664665‑665138                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                        aaCCATCCTGGGCAACAATAATCGGATCGATACGAACTTTGGTGTTGTTTGACGTGTTAAGACGGTCAGCTGCCGCGCGAGAAAGTGAAATAACGCGGTCGTTGCCGTAAGGGCCGCGATCATTAATGCGCACCACGATCATTCGCCCGTTAGCCAGGTTAGTG‑ATTCTGGCGTAGCTGGGGATCGGAAGCGTTGGATGGGCTGCCGTCAGCTGTGTCGGATCGAATGCTTCGCCAGAGGCCATCAGGTTa                                                                                                                                                                                                                                 <  1:140720/251‑1 (MQ=255)    cATCCTGGGCAACAATAATCGGATCGATACGAACTTTGGTGTTGTTTGACGTGTTAAGACGGTCAGCTGCCGCGCGAGAAAGTGAAATAACGCGGTCGTTGCCGTAAGGGCCGCGATCATTAATGCGCACCACGATCATTCGCCCGTTAGCCAGGTTAGTG‑ATTCTGGCGTAGCTGGGGATCGGAAGCGTTGGATGGGCTGCCGTCAGCTGTGTCGGATCGAATGCTTCGCCAGAGGCCATCAGGTTACTg                                                                                                                                                                                                                              <  2:144712/251‑1 (MQ=255)    cATCCTGGGCAACAATAATCGGATCGATACGAACTTAGGAGTTGTTTGACGTGTTAAGAAGGTCAGCTGCCGCGCGAGAAAGTGAAATAACGCGGTCGTTGCCGTAAGGGACGAGATCATTAATGCGCACCACGATCATTCGCCCGTTAGCCAGGTTAGTG‑ATTCTGGAGTAGCTGGGGATCGGAAGCGTTGGATGGGCAGCCGTCAGCTGTGTCGGATCGAATGCTTCGCCAGAGGCCATCAGGTTACTg                                                                                                                                                                                                                              >  1:343723/1‑251 (MQ=255)                 aataatCGGATCGATACGAACTTTGGTGTTGTTTGACGTGTGAAGACGGTCAGCTGCCGCGCGAGAAAGTGAAATAACGCGGTCGTTGCCGTAAGGGCCGCGATCATTAATGCGCACCACGATCATTCGCCCGTTAGCCAGGTTAGTG‑ATTGTGGCGTAGCTGGGGATCGGAAGCGTTGGATGGGCTGCCGTCAGCTGGGTCGGATCGAATGCTTCGCCAGAGGCCATCAGGTTACTGCCGGGTTCGGcat                                                                                                                                                                                                                 <  2:246055/251‑1 (MQ=255)                       cGGATCGATACGAACTTTGGTGTTGTTTGACGTGTTAAGACGGTCAGCTGCCGCGCGAGAAAGTGAAATAACGCGGTCGTTGCCGTAAGGGCCGCGATCATTAATGCGCACCACGATCATTCGCCCGTTAGCCAGGTTAGTG‑ATTCTGGCGTAGCTGGGGATCGGAAGCGTTGGATGGGCTGCCGTCAGCTGTGTCGGATCGAATGCTTCGCCAGAGGCCATCAGGTTACTGCCGGGTTCGGcatcat                                                                                                                                                                                                              >  1:434529/1‑248 (MQ=255)                                          gagttgttTGACGTGTTAAGACGGTCAGCTGCCGCGCGAGAAAGTGAAATAACGCGGTCGTTGCCGTAAGGGCCGCGATCATTAATGCGCACCACGATCATTCGCCCGTTAGCCAGGTTAGTG‑ATTCTGGCGTAGCTGGGGATCGGAAGCGTTGGATGGGCTGCCGTCAGCTGTGTCGGATCGAATGCTTCGCCAGAGGCCATCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTg                                                                                                                                                                                        >  1:227552/3‑251 (MQ=255)                                                 ttGACGTGTTAAGACGGTCAGCTGCCGCGCGAGAAAGTGAAATAACGCGGTCGTTGCCGTAAGGGCCGCGATCATTAATGCGCACCACGATCATTCGCCCGTTAGCCAGGTTAGTG‑ATTCTGGCGTAGCTGGGGATCGGAAGCGTTGGATGGGCTGCCGTCAGCTGTGTCGGATCGAATGCTTCGCCAGAGGCCATCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAAt                                                                                                                                                                                 >  1:89154/1‑251 (MQ=255)                                                        gtTAAGACGGTCAGCTGCCGCGCGAGAAAGTGAAATAACGCGGTCGTTGCCGGAAGGGCCGCGATCATTAATGCGCACCACGATCATTCGCCCGTTAGCCAGGTTA‑TGTATTCTGGCGTAGCGGGGGATCGGAAGCGTTGGATGGGCTGCCGTCAGCTGTGTCGGATCGAATGCTTCGCCAGAGGCCATCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGAc                                                                                                                                                                           <  1:415742/250‑1 (MQ=255)                                                           aaGACGGTCAGCTGCCGCGCGAGAAAGTGAAATAACGCGGTCGTTGCCGTAAGGGCCGCGATCATTAATGCGCACCACGATCATTCGCCCGTTAGCCAGGTTAGTG‑ATTCTGGCGTAGCTGGGGATCGGAAGCGTTGGATGGGCTGCCGTCAGCTGTGTCGGATCGAATGCTTCGCCAGAGGCCATAAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGAt                                                                                                                                                                       >  1:362518/1‑251 (MQ=255)                                                            atgngttCAGCTCCCGCGCAAGAAATTGAAATACCGCGGGCGTGGCCGTAAGGGCCGCGATTATTAATGCGCCCCACGATCATTCGCCCGTTAGCCAGGTTAGTG‑ATTCTGGCGTAGCTGGGGATCGGAAGCTTTGGATGGGCTGCCGGCAGCTGTGTCGGATCGAATGCTTCGCCAGAGGCCATCAGGTTACTGCCGGTTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATc                                                                                                                                                                      <  2:405709/245‑1 (MQ=255)                                                             gacgggaaGCTGCCGCGCGAGAAAGTGAAAAAACGCGGTCGTTGCCGTAAGGGCCGCGATCATTAAGGCGCACCCCGATCATTCGCCCGTTAGCCAGGTTAGTG‑ATTCTGGCGTAGCTGGGGATCGGAGGCGTTGGATGGGCTGCCGTCAGCTGTGTCGGATCGAATGCTTCGCCAGAGGCCATCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATcc                                                                                                                                                                     <  2:340155/244‑1 (MQ=255)                                                             gACGGTCGGCTGCCGCGCGAGAAAGTGAAATACCGCGGCCGTTGCCGTAAGGGCCGGGATCATTAATGCGCCCCACGATCATCCGCCCGTTAGCCAGGTTAGTG‑ATTCTGGCGTAGCTGGGGATCGGAAGCTTTGGATGGGCTGCCGTCAGCTGTGTCGGATCGAATGCTTCGCCAGAGGCCATCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATcc                                                                                                                                                                     <  2:89154/251‑1 (MQ=255)                                                                ggTCAGCTGCCGCGCGAGAAAGTGAAATAACGCGGTCGTTGCCGTAAGGGCCGCGATCATTAATGCGCACCACGATCATTCGCCCGTTAGCCAGGTTAGTG‑ATTCTGGCGTAGCTGGGGATCGGAAGCGTTGGATGGGCTGCCGTCAGCTGTGTCGGATCGAATGCTTCGCCAGAGGCCATCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGc                                                                                                                                                                  >  2:397451/1‑251 (MQ=255)                                                                   cAGCTGCCGCGCGAGAAAGTGAAATAACGCGGTCGTTGCCGTAAGGGCCGCGATCATTAATGCGCACCACGATCATTCGCCCGTTAGCCAGGTTAGTG‑ATTCTGGCGTAGCTGGGGATCGGAAGCGTTGGATGGGCTGCCGTCAGCTGTGTCGGATCGAATGCTTCGCCAGAGGCCATCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACg                                                                                                                                                               >  1:179816/1‑251 (MQ=255)                                                                    aGCTGCCGCGCGAGAAAGTGAAATAACGCGGTCGTTGCCGTAAGGGCCGCGATCATTCATGCGCACCACGATCATTCGCCCGTTAGCCAGGTTAGTG‑ATTCTGGCGTAGCTGGGGATCGGTAGCGTTGGATGGGCTGCCGTCAGCTGTGTCGGATCGAATGCTTCGCCAGAGGCCATCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACg                                                                                                                                                               >  1:30762/1‑250 (MQ=255)                                                                        gCCGCGCGAGAAAGTGAAATAACGCGGTCGTTGCCGTAAGGGCCGCGATCATTAATGCGCCCCACGATCATTCGCCCGTTAGCCAGGTTAGTG‑ATTCTGGCGTAGCTGGGGATCGGAAGCGTTGGATGGGCTGCCGTCAGCTGTGTCGGATCGAATGCTTCGCCAGAGGCCCTCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGAtt                                                                                                                                                            <  1:376603/249‑1 (MQ=255)                                                                        gCCGCGCGAGAAAGTGAAATAACGCGGTCGTGGCCGTAAGGGCCGCGATCATTAATGCGCACCACGATCATTCGCCCGTTAGCCAGGTTAGTG‑ATTCTGGCGTAGCTGGGGATCGGAAGCGTTGGATGGGCTGCCGTCAGCTGTGTCGGATCGAATGCTTCGCCAGAGGCCATCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGAtttt                                                                                                                                                          <  2:398075/251‑1 (MQ=255)                                                                                 gaAAGTGAAATAAAGCGGTCGTTGCCGTAAGGGCCGCGATCATTAATGCGCACCACGATCATTCGCCCGTTAGCCAGGTTAGTG‑ATTCTGGCGTAGCTGGGGATCGGAAGCGTTGGATGGGCTGCCGTCAGCTGTGTCGGATCGAATGCTTCGCCAGAGGCCATCAGGTTACTGCCGGGTTCGGCATCATCGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGCCGGAACCTGCACGATTTTGTAGCtttt                                                                                                                                                 >  2:430695/1‑251 (MQ=255)                                                                                  aaggggAAATAACGCGGGCGTTGCCGTAAGGGCCGCGATCATTAATGCGCACCACGATCATTCGCCCGTTAGCCAGGTTAGTG‑ATTCGGGCGTAGCTGGGGATCGGAAGTGTGGGATGGGCTGCCGTCAGCTGTTGCGGATCGAATGCTTCGCCAGAGGCCATCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTa                                                                                                                                                <  2:434529/246‑1 (MQ=255)                                                                                     gTGAAATAACGCGGTCGTTGCCGTAAGGGCCGCGATCATTAATGCGCACCACGATCATTCGCCCGTTAGCCAGGTTAGTG‑ATTCTGGCGTAGCTGGGGATCGGAAGCGTTGGATGGGCTGCCGTCAGCTGTGTCGGATCGAATGCTTCGCCAGAGGCCATCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCg                                                                                                                                             >  1:276009/1‑251 (MQ=255)                                                                                            aaCGCGGTCGTTGCCGTAAGGGCCGCGATCATTAATGCGCACCACGATCATTCGCCCGTTAGCCAGGTTAGTG‑ATTCTGGCGTAGCTGGGGATCGGAAGCGTTGGATGGGCTGCCGTCAGCTGTGTCGGATCGAATGCTTCGCCAGAGGCCATCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGCTTTTGTAGCTTTTACCGTCGCGCt                                                                                                                                      >  1:354786/1‑251 (MQ=255)                                                                                                 ggTCGTTGCCGTAAGGGCCGCGATCATTAATGCGCACCACGATCATTCGCCCGTTAGCCAGGTTAGTG‑ATTCTGGCGTAGCTGGGGATCGGAAGCGTTGGATGGGCTGCCGTCAGCTGTGTCGGATCGAATGCTTCGCCAGAGGCCATCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTa                                                                                                                                  >  1:345694/1‑250 (MQ=255)                                                                                                    cGTTGCCGTAAGGGCCGCGATCAATAATGCGCACCACGATCATTCGCCCGTTAGCCAGGTTAGTG‑ATTCTGGCGTAGCTGGGGATCGGAAGCGTTGGATGGGCTGCCGTCAGCTGTGTCGGATCGAATGCTTCGCCAGAGGCCATCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATcc                                                                                                                              >  2:290292/1‑251 (MQ=255)                                                                                                         ccGTAAGGGCCGCGATCATTAATGCGCACCACGATCATTCGCCCGTTAGCCAGGTTAGTG‑ATTCTGGCGTAGCTGGGGATCGGAAGCGTTGGATGGGCTGCCGTCAGCTGTGTCGGATCGAATGCTTCGCCAGAGGCCATCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGAtt                                                                                                                         <  1:136626/251‑1 (MQ=255)                                                                                                           gTAAGGGCCGCGATCATTAATGCGCACCACGATCATTCGCCCGTTAGCCAGGTTAGTG‑ATTCTGGCGTAGCTGGGGATCGGAAGCGTTGGATGGGCTGCCGTCAGCTGTGTCGGATCGAATGCTTCGCCAGAGGCCATCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTg                                                                                                                       >  2:92804/1‑251 (MQ=255)                                                                                                                     cgATCATTAATGCGCACCACGATCATTCGCCCGTTAGCCAGGTTAGTG‑ATTCTGGCGTAGCTGGGGATCGGAAGCGTTGGATGGGCTGCCGTCAGCTGTGTCGGATCGAATGCTTCGCCAGAGGCCATCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGGAATCCTGATTTGCCGTCGCGtt                                                                                                             >  2:156617/1‑251 (MQ=255)                                                                                                                        tCATTAATGCGCACCACGATCATTCGCCCGTTAGCCAGGTTAGTG‑ATTCTGGCGTAGCTGGGGATCGGAAGCGTTGGATGGGCTGCCGTCAGCTGTGTCGGATCGAATGCTTCGCCAGAGGCCATCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTc                                                                                                            >  2:314105/1‑249 (MQ=255)                                                                                                                             aaTGCGCACCACGATCATTCGCCCGTTAGCCAGGTTAGTG‑ATTCTGGCGTAGCTGGGGATCGGAAGCGTTGGATGGGCTGCCGTCAGCTGGGTCGGATCGAATGCTTCGCCAGAGGCCATCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGt                                                                                                      <  1:212363/250‑1 (MQ=255)                                                                                                                                      cacGATCATTCGCCCGTTAGCCAGGTTAGTG‑ATTCTGGCGTAGCTGGGGATCGGAAGCGTTGGATGGGCTGCCGTCAGCTGTGTCGGATCGAATGCTTCGCCAGAGGCCATCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAAcgcg                                                                                            >  1:409847/1‑251 (MQ=255)                                                                                                                                            cATTCGCCCTTTAGCCAGGTTAGTG‑ATTCTGGCGTAGCTGGGGATCGGAAGCGTTGGATGGGCTGCCGTCAGCTGTGTCGGATCGAATGCTTCGCCAGAGGCCATCAGGTTCCTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGAGGGTcc                                                                                       <  1:308591/250‑1 (MQ=255)                                                                                                                                                 gCCCGTTAGCCAGGTTAGTG‑ATTCTGGCGTAGCTGGGGATCGGAAGCGTTGGATGGGCTGCCGTCAGCTGTGTCGGATCGAATGCTTCGCCAGAGGCCATCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCg                                                                                 >  2:57094/1‑251 (MQ=255)                                                                                                                                                    cGTTAGCCAGGTTAGTG‑ATTCTGGCGTAGCTGGGGATCGGAAGCGTTGGATGGGCTGCCGTCAGCTGTGTCGGATCGAATGCTTCGCCAGAGGCCATCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCt                                                                               >  1:53925/1‑250 (MQ=255)                                                                                                                                                         gCCAGGTTAGTG‑ATTCTGGCGTAGCTGGGGATCGGAAGCGTTGGATGGGCTGCCGTCAGCTGTGTCGGATCGAATGCTTCGCCAGAGGCCATCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTc                                                                         >  1:222141/1‑251 (MQ=255)                                                                                                                                                          ccAGGTTAGTG‑ATTCTGGCGTAGCTGGGGATCGGAAGCGTTGGATGGGCTGCCGTCAGCTGTGTCGGATCGAATGCTTCGCCAGAGGCCATCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTc                                                                         >  1:19361/1‑250 (MQ=255)                                                                                                                                                            aGGTTAGTG‑ATTCTGGCGTAGCTGGGGATCGGAAGCGTTGGATGGGCTGCCGTCAGCTGTGTCGGATCGAATGCTTCGCCAGAGGCCATCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAAc                                                                      >  2:197203/1‑251 (MQ=255)                                                                                                                                                                     attttggggttggttGGGGATAGGATGCGTTGGAGGGGCTGCTGTCATCTGTGTCGTTACAAATGCTTCGCAAGAGTCCTACAGGTTCCTGCGGGGTTCGGCATCATAGATTGCCGCTAGTCCCGCCTGTCTAAATCGAGACGGTTCCTACTCGATTTTGTAGCTGTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCTGCCCTGCTTATTTCAACTATAGGGc                                                              <  2:30762/237‑1 (MQ=255)                                                                                                                                                                       ttCTGGCGTAGCTGGGGATCGGAAGCGTTGGATGGGCTGCCGTCAGCTGTGTCGGATCGAATGCTTCGCCAGAGGCCATCAGTTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGcc                                                             <  1:156617/250‑1 (MQ=255)                                                                                                                                                                                     gggATCGGAAGCGTTGGAGCGCCGGCAGTCAGCTGTGTCGGATCGAATGATTCGCCAGAGGCCATCAGGTTACTGCCGGGTTCGGCATCCTAGATTGCAGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGCTTTTCTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGACGTCGCGTTCCGTGATTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAgg                                              >  2:176227/1‑251 (MQ=255)                                                                                                                                                                                              aGCGTTGGATGGGCTGCCGTCAGCTGTGTCGGATCGAATGCTTCGCCAGAGGCCATCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTac                                     >  2:399712/1‑251 (MQ=255)                                                                                                                                                                                                    ggATGGGCTGCCGTCAGCTGTGTCGGATCGAATGCTTCGCCAGAGGCCATCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTAc                               <  1:397451/251‑1 (MQ=255)                                                                                                                                                                                                               cGTCAGCTGTGTCGGATCGAATGCTTCGCCAGAGGCCATCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGAGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCt                    >  1:65493/1‑251 (MQ=255)                                                                                                                                                                                                                gTCAGCTGTGTCGGATCGAATGCTTCGCCAGAGGCCATCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTg                   >  1:52695/1‑251 (MQ=255)                                                                                                                                                                                                                  cAGCTGTGTCGGATCGATTGCTTCGCCAGAGGCCATCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGAc                 <  1:368513/251‑1 (MQ=255)                                                                                                                                                                                                                     gtgtgTCGGATCGATGGCTTCGCCAGAGGCCATCAGTTTACTGCCGGGTTCGGCATAATAATTTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACcat              <  2:345694/250‑1 (MQ=255)                                                                                                                                                                                                                       ctgtcggcctgAATGCTTCGCCATAGGCCATCCGGTTACTGCCGGGTTCGGCATCAGATGTTGCCTCCAGTCCCGCCTGGCTAAATCGAGTCGGATCCTGCACGATTCTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCATTGGTTCGAAACGCCGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACcat              <  2:77242/239‑1 (MQ=255)                                                                                                                                                                                                                           cGGATCAAATGCTTCCCCAGATGCCATCAGTTTTCTGCCGGGTTCGGCTTCATAGATGGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTAGAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCg         <  1:196685/250‑1 (MQ=255)                                                                                                                                                                                                                             gATCGAATGCTTCGCCAGAGGCCATCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCt                                                                               <  1:75002/178‑1 (MQ=255)                                                                                                                                                                                                                                 gAATGCTTCGCCAGAGGCCATCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTAc  <  1:311337/251‑1 (MQ=255)                                                                                                                                                                                                                                 gAATGCTTCGCCAGAGGCCATCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTAc  >  1:120430/1‑251 (MQ=255)                                                                                                                                                                                                                                               gggCCATCAGGTTACTGCCGGGTCCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCt                          <  1:204814/212‑1 (MQ=255)                                                                                                                                                                                                                                               agGCCATCAGGTTACTGCCGGGTTCGGCATCATAGATTGCAGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGCTTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATCGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCt                          >  2:204814/1‑213 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                        AACCATCCTGGGCAACAATAATCGGATCGATACGAACTTTGGTGTTGTTTGACGTGTTAAGACGGTCAGCTGCCGCGCGAGAAAGTGAAATAACGCGGTCGTTGCCGTAAGGGCCGCGATCATTAATGCGCACCACGATCATTCGCCCGTTAGCCAGGTTAGTG‑ATTCTGGCGTAGCTGGGGATCGGAAGCGTTGGATGGGCTGCCGTCAGCTGTGTCGGATCGAATGCTTCGCCAGAGGCCGTCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTAC  >  NC_000913/664665‑665138

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