Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,213,775	G→A	18.6%	F116F (TTC→TTT)	yjaB ←	GNAT‑family putative N‑acetyltransferase, acetyl coenzyme A‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,213,775	0	G	A	18.6%	99.2 / 2.5	53	F116F (TTC→TTT)	yjaB	GNAT‑family putative N‑acetyltransferase, acetyl coenzyme A‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (36/6);  new base A (2/8);  total (39/14)
Fisher's exact test for biased strand distribution p-value = 1.38e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCAC‑CGGGATGAGAAAAAATCGCCTACGCCCCCACATA‑CGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCT‑TATAGAA‑CCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCG‑CATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAA  >  NC_000913/4213558‑4213991                                                                                                                                                                                                                             |                                                                                                                                                                                                                           aTGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCAC‑CGGGATGAGAAAAAATCGCCTACGCCCCCACATA‑CGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCT‑TATAGAA‑CCCAACCg                                                                                                                                                                                                                 >  2:359851/1‑228 (MQ=255)      catccggtgtttctcgtagagtgTGCGGTCATGGCAATGCTTAAGCAGAAATAACCGTGTCCACTTGGGGGTGTAACAAACTCATAGTGCAGCCCCC‑CGGGATGAGAAAAAACCGCCTCCGCCCCCACATC‑CGACAGTTTCAGCACCGGATACGTTTTCCCCAAATCGTCCCCCCAAGACCGCCCCGAAACTTAAAACCCCACCTCCT‑TATAAAC‑CCCAACCGCCTGc                                                                                                                                                                                                            <  1:618342/206‑1 (MQ=255)      cAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCAC‑CGGGATGAGAAAAAATCGCCTACGCCCCCACATA‑CGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCT‑TATAAAA‑CCCAACCGCCt                                                                                                                                                                                                              >  2:142893/1‑226 (MQ=255)         tCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCAC‑CGGGATGAGAAAAAATCGCCTACGCCCCCACATA‑CGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCT‑TATAGAA‑CCCAACCGCCTGCTCa                                                                                                                                                                                                         >  2:241414/1‑228 (MQ=255)         tCAGCGGTACTCAGTGATAGTGCGGTCATGGAAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCAC‑CGGGATGAGAAAAAATCGCCTACGCCCCCACATA‑CGCCAGATTAAGCAACGGAGACGGCTTCCCCAAATCGTCAACCCCAAACCGTCCCGTAACCTTAAACCCCCCCTTAC‑TTTAAAA‑CCCAAACCCCTggccc                                                                                                                                                                                                         >  1:165399/1‑224 (MQ=255)                     aGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCAC‑CGGGATGAGAAAAAATCGCCTACGCCCCCACATA‑CGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCAACCTTCT‑TATAGAA‑CCCAACCGCCTGCTCATTTTCCTAATTa                                                                                                                                                                                             >  1:293730/1‑228 (MQ=255)                       tGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCAC‑CGGGATGAGAAAAAATCGCCTACGCCCCCACATA‑CGCCAGATTCAGCAACGGATACGGGTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCT‑TATAGAA‑CCCAACCGCCTCCTCCTTTTGCTCCTTAAc                                                                                                                                                                                           >  2:620765/1‑228 (MQ=255)                             tGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCAC‑CGGGATGAGAAAAAATCGCCTACGCCCCCACATA‑CGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCT‑TATAGAA‑CCCAACCGCCTGCTCATTTTGCTCATTAACgttggt                                                                                                                                                                                     >  2:316766/1‑228 (MQ=255)                               cGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCAC‑CGGGATGAGAAAAAATCGCCTACGCCCCCACATA‑CGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCT‑TATAGAA‑CCCAACCGCCTGCCAATTTTCATCATTAACgttggtt                                                                                                                                                                                    >  2:222179/1‑227 (MQ=255)                                 gTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCAC‑CGGGATGAGAAAAAATCGCCTACGCCCCCACATA‑CGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCT‑TATAGAA‑CCCAACCGCCTGCTCATTTTTCTCATTAACGTTGGTTGt                                                                                                                                                                                  >  2:52353/1‑227 (MQ=255)                                      ggaattcttaaaagaaaaaTAATCGTGCCACTAGTGGGGGGCACAAACCCAGAAG‑TCCGCCCACACGGTAGGAAAAAAAACTGCCCACGCCCCCACAGC‑ACCCAGTTCCACCAACGGATACGGTTCCCCCAAACCGTCCCCCCAAGAGCTTCCCGTAACTTAAAAACCCACTTTCT‑TATAGAA‑CCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTc                                                                                                                                                                            <  1:138912/212‑1 (MQ=255)                                       gCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCAC‑CGGGATGAGAAAAAATCGCCTACGCCCCCACATA‑CGCCAGATTCAGCAACGGATACGGTTTCCCCAAACCGTCCACCTCCGAGCGTCCCGCAACCTTAAAACCCACCTTCT‑TATAGAA‑CCCAACCCCCTGCCCCTTTTTCTCATTAACGTTAGTTGTCAGTTcc                                                                                                                                                                           >  1:171693/1‑228 (MQ=255)                                            gCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCAC‑CGGGATGAGAAAAAATCGCCTACGCCCCCACATA‑CGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCT‑TATAGAA‑CCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTg                                                                                                                                                                       >  2:390599/1‑227 (MQ=255)                                                       aataatCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCAC‑CGGGATGAGAAAAAATCGCCTACGCCCCCACATA‑CGCCAGCTTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAACCCCACCTTCT‑TATAGAA‑CCCAACCGCCTGCTCATTTTGCTCATTAACGTTGCTTGTCAGCTCCTGGGCCATCGAGAgc                                                                                                                                                            >  2:430126/1‑227 (MQ=255)                                                            tCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCAC‑CGGGATGAGAAAAAATCGCCTACGCCCCCACATA‑CGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCT‑TATAGAA‑CCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGGTGTCAGTTCCGGGGCCACGCAAAGCG‑CATg                                                                                                                                                      >  1:64686/1‑227 (MQ=255)                                                                  cACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCAC‑CGGGATGAGAAAAAATCGCCTACGCCCCCACATA‑CGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCACAACGTCCCGTAACCTTAAAACCCACCTTCT‑TATAGAA‑CCCAACCCCCTCATCCTTTTGCTCATTAACGTTGGTTTTCAGTTCCGGTCCCacccaaaccgcctccccccc                                                                                                                                                 >  1:307807/1‑208 (MQ=255)                                                                       ttGGTGGGTACTAAACCTGAAGTTCAGCCCAC‑CGGGATGAGAAAAAATCGCCTACGCCCCCACATA‑CGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCT‑TATAGAA‑CCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCAGAAGCC‑AATGCTccaccaaga                                                                                                                                           >  2:514986/1‑224 (MQ=255)                                                                               tACTAAACCTGAAGTTCAGCCCAC‑CGGGATGAGAAAAAATCGCCTACGCCCCCACATA‑CGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCT‑TATAGAA‑CCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTTTCAGTTCCGTTGCCATCGAGAGCG‑CCTGCTCCACCAGCACCCGACCt                                                                                                                                   >  2:241306/1‑227 (MQ=255)                                                                                  tAAACCTGAAGTTCAGCCCAC‑CGGGATGAGAAAAAATCGCCTACGCCCCCACATA‑CGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCT‑TATAGAA‑CCCAACCGCCTGCTCATTTTTCTCATTAACGTTGGTTGTCAGTTCCCAGGCCCTCCAGAGCG‑CATGCTCCCCCAGCACCCGACCTAccc                                                                                                                               >  1:46884/1‑226 (MQ=255)                                                                                  tAAACCTGAAGTTCAGCCCAC‑CGGGATGAGAAAAAATCGCCTACGCCCCCACATA‑CGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCT‑TATAGAA‑CCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCG‑CATGCTCCACCAGCACCCGACCTAcgc                                                                                                                               >  2:360265/1‑228 (MQ=255)                                                                                    aaCCTGAAGTTCAGCCCAC‑CGGGATGAGAAAAAATCGCCTACGCCCCCACATA‑CGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCT‑TATAGAA‑CCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAACG‑CATGCTCCACCAGCACCCCCCCTAcccca                                                                                                                             >  1:157821/1‑227 (MQ=255)                                                                                         gAAGTTCAGCCCAC‑CGGGATGAGAAAAAATCGCCTACGCCCCCACATA‑CGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCT‑TATAGAA‑CCCAACCGCCTGCTCCTTTTGCTCCTTAACCTTGGTTCGCAATTCCCGGGCCGACCA‑ATCGCCACGCTCCCCCAGCACCCGAccaccccccacacc                                                                                                                        >  1:114064/1‑216 (MQ=255)                                                                                         gAAGTTCAGCCCAC‑CGGGATGAGAAAAAATCGCCTACGCCCCCACATA‑CGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCT‑TATAGAA‑CCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCG‑CCTGCTCCACCAGCACCCGACCTAggccccgccc                                                                                                                        >  2:201667/1‑222 (MQ=255)                                                                                         gAAGTTCAGCCCAC‑CGGGATGAGAAAAAATCGCCTACGCCCCCACATA‑CGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCT‑TATAGAA‑CCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCG‑CATGCTCCACCAGCCCCCGACCTACGCCACAGcc                                                                                                                        >  2:130881/1‑228 (MQ=255)                                                                                         gAAGTTCAGCCCAC‑CGGGATGAGAAAAAATCGCCTACGCCCCCACATA‑CGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCT‑TATAGAA‑CCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCCAAAGCG‑CATGCTCCACCAGCACCCGACCTCCCCCGCAGcc                                                                                                                        >  2:524471/1‑228 (MQ=255)                                                                                         gAAGTTCAGCCCAC‑CGGGATGAGAAAAAATCGCCTACGCCCCCACATA‑CGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGAAACCTTAAAACCCACCTTCT‑TATAGAA‑CCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGCCAGTTCCGGTGCCCTCGAGATCG‑CATGCTCCACCAGCCACGGACCTACGCCCCAGc                                                                                                                         >  2:355158/1‑227 (MQ=255)                                                                                         gAAGTTCAGCCCAC‑CGGGATGAGAAAAAATCGCCTACGCCCCCACATA‑CGACAGATTCAGCAACGGATAAGGTTTCCCCAAATCGTCCACCTGAGAGCGTCCCGTAACCTTAAAAACCACCTACT‑TATACAA‑GCCAACCGCCTGCTCATTTTGATCATTACCGCCGGATGTCAGTTCCAGTGCCATCGAGCGCG‑CATACTCCACCCGCAACGGCCCT‑CCCCGCAGCcg                                                                                                                       >  2:403159/1‑228 (MQ=255)                                                                                            gtacaccccAC‑GGGGATGAGAAAAAATCGCCTACGCCCCCACTTA‑CGCCAGATTCAGCAACGGATACGGTTCCCCAAAATCGTCCACCTCAGAGCGTCCCGTAACCTAAAAACCCACCTTCT‑TATAGAA‑CCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCG‑CATGCTCCACCAGCACCCGACCTACGCCGCAGCcgcg                                                                                                                     <  2:114064/222‑1 (MQ=255)                                                                                            gTTCAGCCCAC‑CGGGATGAGAAAAAATCGCCTACGCCCCCACATA‑CGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCT‑TATAGAA‑CCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGACAGCG‑CATGCTCCACCAGCACCCGACCTACGCCGCAGCcgcg                                                                                                                     >  2:655921/1‑228 (MQ=255)                                                                                                                     aaatccccttccccccACATACCGCCGATTTCAGCAACGATTACGGTTCCCCCAAACCGTCCACCTCAGGGCGCCCCGAACCCTAAAAACCCACCTTCT‑AATAAAA‑CCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCG‑CATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAgcg                                                                                              <  1:316766/217‑1 (MQ=255)                                                                                                                     aaacccccacccccccccaccccGCCAGATACCACAACGGATCCGTTTCCCCAAAACCCCCCACCCCAGGGCGTCCCGTAACCTTAAACCCCCCCTTCT‑TAAAAAA‑CCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCG‑CATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAgcg                                                                                              <  1:620765/206‑1 (MQ=255)                                                                                                                                                 aTTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCT‑TATAGAA‑CCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCTAGAGCG‑CATGCTCCACCAGCACCCGACCTACGCCGCCGCCGCCACCATCAGGATCGCTAAACAGCGCACCCatatcgcgccccactaccacaat                                                                  >  1:519461/1‑209 (MQ=255)                                                                                                                                                 aTTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCT‑TATAGAA‑CCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGCGAGCG‑CATGCTCCACCAGCACCCGACCTACCCCGCAGCCCCCCACATCAGGACCGATCAACAGCGCACCCATCTCCTGCCCACTAACCAACaa                                                                  >  1:404918/1‑227 (MQ=255)                                                                                                                                                 aTTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCT‑TATAGAA‑CCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCG‑CATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCCCCCCACTTCGAAACaa                                                                  >  1:261302/1‑227 (MQ=255)                                                                                                                                                               aCGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCT‑TATAGAA‑CCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCG‑CATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACCGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCt                                                    >  2:383857/1‑228 (MQ=255)                                                                                                                                                                      ccccAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCT‑TATAGAA‑CCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTTCCATCGAGAGCG‑CATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAACCAGCGCATCCATATCCTGCCCACTTAg                                                                        >  2:351674/1‑201 (MQ=255)                                                                                                                                                                      cccaAACTCGCCCACCTCGGACTCCCCCGTAACATTAAACCCCACCTTCT‑TA‑AGAACCCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCG‑CATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCg                                             <  2:307807/228‑1 (MQ=255)                                                                                                                                                                        caaagcgcccccccTCAGTCCGTCCGTTAACTTTAAACCCCGCCTTCT‑TATGAAC‑CCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCG‑CATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCg                                             <  2:46884/214‑1 (MQ=255)                                                                                                                                                                             tCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCT‑TATAGAA‑CCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCG‑CATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTa                                      >  2:248638/1‑228 (MQ=255)                                                                                                                                                                                 ccccccaagcccctCCCGTAACCTAAAACCCCCCCTTCT‑AAACAAC‑CCCCACCGCCTGCTCATTGTGCTCATTAACGTTGGTTGCCCGTTCCGGGGCCATCCAGAACG‑CATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACcg                                  <  1:52353/215‑1 (MQ=255)                                                                                                                                                                                 cccccTCAGAGCTCCCCGTAACC‑TCAAACCCTCCTTTTGAATAAAA‑CCCACCCGCCTGCTCATTGTGCGCCTTAACGTTGGTTGTCAGTTCCTGTGCCATCGAGAGCG‑CATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACcg                                  <  1:130881/225‑1 (MQ=255)                                                                                                                                                                                 ccACCCCAGAGCCTCCCTTACCTTTAAACCCCCCTTTTT‑TATGAAA‑CCAAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCG‑CATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACcg                                  <  2:157821/228‑1 (MQ=255)                                                                                                                                                                                          agCGTCCCGTAACCTTAAAACCCACCTTCT‑TATAGAA‑CCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCG‑CATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGAccacca                         >  2:610335/1‑224 (MQ=255)                                                                                                                                                                                          agCGTCCCGTAACCTTAAAACCCACCTTCT‑TATAGAA‑CCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCG‑CATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTCAGCAACATAAACCCAACCGGCTGATCCCGCCCCTTAACCACCAccccc                          >  1:235097/1‑225 (MQ=255)                                                                                                                                                                                          agCGTCCCGTAACCTTAAAACCCACCTTCT‑TATAGAA‑CCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCG‑CATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCCCATCCACATCCTGCCCCCTTAGCAACATAAATCCAACCGGCCGATCCCGCCCCTTAACCCCCAccccaa                         >  1:392855/1‑225 (MQ=255)                                                                                                                                                                                                          aaaaCCCACCTTCT‑TATAGAA‑CCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCG‑CATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCa          >  2:561581/1‑227 (MQ=255)                                                                                                                                                                                                                ccctttct‑tATAGCC‑CCCAGCCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCG‑CATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCaggaagga   <  2:64686/224‑1 (MQ=255)                                                                                                                                                                                                                ccccTTTT‑TACAGAA‑CCCAACCGCCTGGTCATTTTGCTCATTAACGGTGGTTTTCAGTTCCGGTGCCATCGAGAGCG‑CATGCTCCCCCAGCCCCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGACCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCaggaagga   <  1:514986/226‑1 (MQ=255)                                                                                                                                                                                                                ccccTTTT‑TAAAAAC‑CCCAACCGCCTGCGCATTTTTCCCATTACCGGTGGTTGTCATTTCCGGTGCCATCGAGAGCG‑CATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCaggaagga   <  1:222179/226‑1 (MQ=255)                                                                                                                                                                                                                 cccTTTT‑TAAAGAC‑CCCAACCGCCTGCCCATTTTGCTCATTAACGTTGGGTGTCAGTTCCGGTGCCATCGAGAGCG‑CATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGAACGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAggaaggaa  <  2:404918/227‑1 (MQ=255)                                                                                                                                                                                                                 aCCTTCT‑TATAGAA‑CCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCG‑CATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCCTATGCTGCCCACTTAGCAACATAAATCCAACCGCCGGTTCCCGCTCATTAAACCCGCACCCCAACGGACCTTcagcacggaggga   >  2:255859/1‑215 (MQ=255)                                                                                                                                                                                                                 aCCTTCT‑TATAGAA‑CCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCG‑CATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCaggaagga   >  2:5432/1‑227 (MQ=255)                                                                                                                                                                                                                 aCCTTCT‑TATAGAA‑CCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCG‑CATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACACCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGCCTGATCCCCCTCCTTACCCGCGACCCACAACGACGCTTCCGTCgggaagga   >  2:539145/1‑227 (MQ=255)                                                                                                                                                                                                                             |                                                                                                                                                                                                                           ATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCAC‑CGGGATGAGAAAAAATCGCCTACGCCCCCACATA‑CGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCT‑TATAGAA‑CCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCG‑CATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAA  >  NC_000913/4213558‑4213991

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