Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	10,696	C→T	A221T (GCC→ACC)	yaaW ←	putative enzyme‑specific chaperone YaaW

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	10,696	0	C	T	98.0%	151.1 / NA	49	A221T (GCC→ACC)	yaaW	putative enzyme‑specific chaperone YaaW
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (25/23);  minor base G (0/1);  total (25/24)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-01

CCCAGCGGTGCCGGATTAGCCAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCG  >  NC_000913/10457‑10938                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                   cccAGCGGTGCCGGATTAGCCAAATTAGTGTTGCCCATAATTCCTCAAAACTCATCATCGAATGACTGGTGAAATAATTTACCTGCATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTCCGGTTGGAATCGTcc                                                                                                                                                                                                                                         >  2:271676/1‑250 (MQ=255)             ggATTAGCCAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGgcgc                                                                                                                                                                                                                             >  1:294468/1‑251 (MQ=255)                          tAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCg                                                                                                                                                                                                                 <  1:66934/250‑1 (MQ=255)                              gtTGCCTATAATTCCTCAAAAATAATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGGTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTa                                                                                                                                                                                                            <  2:101230/251‑1 (MQ=255)                                      tAATTCCTCAAAACTCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGAGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGAAGTGTCTACGATCTTTTGAGCGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGAGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTACCGATGACCGCTTTAAacccctt                                                                                                                                                                                                    >  1:36071/1‑248 (MQ=255)                                         ttCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGGGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCAttt                                                                                                                                                                                                   >  2:297610/1‑249 (MQ=255)                                            cTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGc                                                                                                                                                                                              >  2:227183/1‑251 (MQ=255)                                                     catcatcGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTAc                                                                                                                                                                                     >  2:175143/1‑251 (MQ=255)                                                      atcatcGCATGAATGGTGAAATCATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAACTTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTAc                                                                                                                                                                                     >  1:341027/1‑250 (MQ=255)                                                      atcatcGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACa                                                                                                                                                                                    <  1:107404/251‑1 (MQ=255)                                                            gaatgaatGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCt                                                                                                                                                                              >  1:268476/1‑251 (MQ=255)                                                            gaatgaatGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCt                                                                                                                                                                              <  1:400502/251‑1 (MQ=255)                                                                   tgGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTa                                                                                                                                                                                                            <  2:33728/214‑1 (MQ=255)                                                                   tgGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTa                                                                                                                                                                                                            >  1:33728/1‑214 (MQ=255)                                                                     gcgaactcatTTCCCTGAATACCTGTAGTGTTTTCAGGGCGCGGCATCATCATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCAACAGGGCCTCCCAGcccc                                                                                                                                                                     >  2:57923/9‑251 (MQ=255)                                                                          ataatTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAACCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAAGGAAAATGTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCAGCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCccgcgcc                                                                                                                                                                <  1:271676/251‑1 (MQ=255)                                                                           taatTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCccgcgccc                                                                                                                                                               >  2:312085/1‑250 (MQ=255)                                                                            aagttCCCGGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGACTACGATCTATTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGAGAGGCGCTGCCGCTTACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCccgcgccgc                                                                                                                                                              >  2:17642/4‑251 (MQ=255)                                                                                             gTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGcc                                                                                                                                               <  1:56466/249‑1 (MQ=255)                                                                                                               ggCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCg                                                                                                                           <  2:294468/251‑1 (MQ=255)                                                                                                                       aatCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTACCCATACGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGCTAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCACGTACGCTCATTGCAGCGTGGGTGCg                                                                                                                   >  2:250076/1‑251 (MQ=255)                                                                                                                                  ggggCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCggg                                                                                                        >  1:389145/1‑251 (MQ=255)                                                                                                                                      cAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAg                                                                                                    <  1:175143/251‑1 (MQ=255)                                                                                                                                                     tttGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCa                                                                                      <  1:375329/250‑1 (MQ=255)                                                                                                                                                     tttGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCGCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCa                                                                                      <  1:370845/250‑1 (MQ=255)                                                                                                                                                       tGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCt                                                                                                                                                                                                                      <  2:424071/120‑1 (MQ=255)                                                                                                                                                       tGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCt                                                                                                                                                                                                                      >  1:424071/1‑120 (MQ=255)                                                                                                                                                        gAGGGGAAAATGAAAATTTTCCCCCGTTTCCGGTACCAGACCTGAGTGGCGCTAACCATCCGGCGCCGGCCGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCTCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAAt                                                                                   <  2:389145/250‑1 (MQ=255)                                                                                                                                                              aaaatgaaaatTTTCCCCGGTTCCCGGTATCAGACCTGAGTGGCGCTAATCATCCCGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGTCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGTGCCTCCCATCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTAAGTTGGCTGGAAAGCAAATGCGacac                                                                            <  2:293042/251‑1 (MQ=255)                                                                                                                                                               aaatgaaaatTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCAGCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTTGGTGCGTAAAATGCGGGTCATTTGGCTGGAAAGCAAATGCGacaca                                                                           <  2:55979/251‑1 (MQ=255)                                                                                                                                                                 atgaaaatGTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACAcc                                                                         <  2:123149/251‑1 (MQ=255)                                                                                                                                                                         tttCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGcc                                                                  <  1:322302/250‑1 (MQ=255)                                                                                                                                                                          ttCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCt                                                                                                                                                                                                               >  2:225890/1‑108 (MQ=255)                                                                                                                                                                          ttCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCt                                                                                                                                                                                                               <  1:225890/108‑1 (MQ=255)                                                                                                                                                                            ccccGGTTTCCGGTATCATACCTGAGTGGCGCTAACCATCCGGCGCAGGCATGCGATTTGCAGTACGGTTGGAAGCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACATGTCCTCCCAGTCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTTATTGCTGCGTGGGTTCGTAAAATGCGGGTCAGTGGTCTGGAAAGCAAATTCGACACACCTTTTGCCaat                                                               <  2:275697/250‑1 (MQ=255)                                                                                                                                                                              ccGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCaataat                                                            >  2:30102/1‑251 (MQ=255)                                                                                                                                                                              ccGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCaataat                                                            >  1:151691/1‑251 (MQ=255)                                                                                                                                                                                       ggTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCAATAATTTGTCTTtc                                                   >  1:380903/1‑251 (MQ=255)                                                                                                                                                                                             aGACCTGAGTGGCGCTAACCCTCCGGCGAAGGCAGGCGATTTGCCGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGATTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAACATGCGGGTCAGTTGGCTGGAAAGCAAATGAGACAAACCTTTTGCCAGTAATTTGTCTTTCATcagc                                             >  2:265066/1‑251 (MQ=255)                                                                                                                                                                                                ccTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCAATAATTTGTCTTTCATCagcagc                                          >  2:325265/1‑251 (MQ=255)                                                                                                                                                                                                                cATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCa                          <  2:394621/251‑1 (MQ=255)                                                                                                                                                                                                                        gcaggcagGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGAGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCTCGCAGCAAACCATGCCCAAGTCCGCTCCTTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGATGGAAAGCAAATGCGACACACCTTTTGCCAATAATTTGTCTTTCATCAGCAGAGGCAGCAGATCTTCCAGCTCATTc                  >  1:172690/1‑251 (MQ=255)                                                                                                                                                                                                                            gcaggcGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCAccc              <  2:328023/251‑1 (MQ=255)                                                                                                                                                                                                                              aggcGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTg            >  1:388640/1‑251 (MQ=255)                                                                                                                                                                                                                              aggcGATTTGCAGGACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTg            >  1:59161/1‑251 (MQ=255)                                                                                                                                                                                                                                   aTTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAACATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATc       >  2:39210/1‑251 (MQ=255)                                                                                                                                                                                                                                      tGCAGTACGGTTGGAATCGTCTCGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCTCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGa     <  1:289551/250‑1 (MQ=255)                                                                                                                                                                                                                                       gCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAAGGCGACACACCGTTTGCCAGTAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGAc    <  1:265066/250‑1 (MQ=255)                                                                                                                                                                                                                                        cAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGGGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACcg  <  1:229262/251‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                   CCCAGCGGTGCCGGATTAGCCAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCG  >  NC_000913/10457‑10938

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