Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	10,703	C→T	T218T (ACG→ACA)	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,703	0	C	T	97.8%	137.7 / ‑4.4	46	T218T (ACG→ACA)	yaaW	putative enzyme‑specific chaperone YaaW
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (20/25);  total (20/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.94e-01

ATTAGCCAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGC‑ACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAA  >  NC_000913/10471‑10947                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                      aTTAGCCAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACAGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTg                                                                                                                                                                                                                                     <  2:280365/251‑1 (MQ=255)          attggtgtAGCCCATAATGCCTCAAAAATCATCATTGAATGAATGTTGAAATAATTTCCCTTAATAACAGTATTGTTTTCAGGGTGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTGGAGGAGAAAATTAAAATTTTCCCCGTTTTACGGTGTCAGCCCTTAGTGGCGCTATCCATCCGGCGCAGGCAGGCGATTGGCAGTACGGCTGGACTTTTCACGCGATAGGCTCTGCCGCTGAcc                                                                                                                                                                                                                            <  2:358446/247‑1 (MQ=255)             aGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACAGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGc                                                                                                                                                                                                                          <  2:150242/250‑1 (MQ=255)                 ttGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACAGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGACCCGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGAAAGGCACTGCCGCTGACCACtttt                                                                                                                                                                                                                      >  2:291269/1‑249 (MQ=255)                          aTTCCTCAAAAATCATCGTCGAATGAATGGTGAAATAATTTCCCTGAATAACAGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCGGGAATTGTCCCGGGATAGGCACTGCGGCTGAGCGCTTTAACCCCAtttg                                                                                                                                                                                                            >  2:153561/1‑250 (MQ=255)                          aTTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACAGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCAAttt                                                                                                                                                                                                             >  1:327721/1‑250 (MQ=255)                                      tcatcatcGAATGAATGGTGAAATAATTTCCCTGAATAACAGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGctgcgccc                                                                                                                                                                                                 >  1:379918/1‑243 (MQ=255)                                          catcGAATGAATGGTGAAATAATTTCCCTGAATAACAGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGC‑GCCCACAgg                                                                                                                                                                                           <  2:379918/251‑1 (MQ=255)                                              gaatgaatGGTGAAATAATTTCCCTGACGAACAGGAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGC‑GCCCACAGGGCCt                                                                                                                                                                                       <  2:319641/251‑1 (MQ=255)                                              gaatgaatGGTGAAATAATTTCCCTGAATAACAGTAGTGTTTGCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTGGAGGGGAAAATGAAAATTTGCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGC‑GCCCACAGGGCCt                                                                                                                                                                                       <  1:450398/251‑1 (MQ=255)                                                         gAAAAAATTTCCCGGAATAACAGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGC‑GCCCACAGGGCCTCCCAGACccg                                                                                                                                                                             <  2:327721/250‑1 (MQ=255)                                                          aaaTAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGC‑ACCTACAGGGCCTCCCAGCCccgcg                                                                                                                                                                           <  1:362727/251‑1 (MQ=255)                                                            ataatTTCCCTGAATAACAGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCCTTTAATGCTGCGACC‑ACAGGGCCTCCCAGACccgcgcc                                                                                                                                                                         >  1:126953/1‑251 (MQ=255)                                                                  tCCCTGAATAACAGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAACTTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGAAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCCTTTAATGCTGC‑GCCCCCAGGTCCTCCCAGACccgcgccgcccac                                                                                                                                                                   >  2:399500/1‑247 (MQ=255)                                                                      tGAATAACAGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGC‑GCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCaa                                                                                                                                                                >  1:243747/1‑250 (MQ=255)                                                                                    gtTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGC‑GCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAgcac                                                                                                                                                 >  1:154412/1‑248 (MQ=255)                                                                                        tCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATTTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCTCAGGCTGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGC‑GCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCt                                                                                                                                              <  2:394315/250‑1 (MQ=255)                                                                                         cAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGC‑GCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCa                                                                                                                                            <  1:447267/251‑1 (MQ=255)                                                                                          aGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGC‑GCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCAt                                                                                                                                           >  1:86522/1‑251 (MQ=255)                                                                                              cgcgGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGC‑GCCCACCGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCCCATtgct                                                                                                                                       >  2:1423/1‑251 (MQ=255)                                                                                               gcgGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGC‑GCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATtgctg                                                                                                                                      <  1:209636/251‑1 (MQ=255)                                                                                                    ataataATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGC‑GCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTgg                                                                                                                                 >  1:314361/1‑251 (MQ=255)                                                                                                              gCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGC‑GCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTaaaa                                                                                                                       <  2:154412/251‑1 (MQ=255)                                                                                                                          gtgtCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGC‑GCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCaatt                                                                                                           <  2:86816/251‑4 (MQ=255)                                                                                                                                   aTCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGC‑GCCCCCAGGGCCTCCCCAACCCGCGCCGCGCAGCAAACCATGCCCAAGCCAGCTCATTGCCGCGTGGGTGTGACAAATGCGGATCaatatgctggaa                                                                                                   >  2:198579/1‑239 (MQ=255)                                                                                                                                     cTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGC‑GCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGc                                                                                                >  1:102336/1‑251 (MQ=255)                                                                                                                                            ggggAAAATGAAAATTTTCCCCGGATTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCGGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGAAATGCCGCTGAACGCTTTTACCCCATTTAATGCTGC‑GCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCACAAGCACGCTCATTGCTGCGTGGGTGGGTAAAATGAGGGTAAATTGGCTGGAAAGGAAATg                                                                                           >  2:64544/1‑249 (MQ=255)                                                                                                                                                 aaatgaaaatTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAAcc                                                                                                                                                                                                                   <  1:436924/125‑1 (MQ=255)                                                                                                                                                 aaatgaaaatTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAAcc                                                                                                                                                                                                                   >  2:436924/1‑125 (MQ=255)                                                                                                                                                      aaaatTTTCCCCGGTTTACGGTGTCATACCGTAGTGGCTCGACCCACCCGGCTCAGGCTGGCGATTTGCAGTACGGCTGGAATTGGCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGC‑GCCCACAGGGCCTCCCAGACCGGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCttt                                                                               <  2:479824/251‑1 (MQ=255)                                                                                                                                                       aaatTTTCCCCGGTTTACGGTGTCCGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGC‑GCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCtttt                                                                              <  2:102336/251‑1 (MQ=255)                                                                                                                                                         atTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGC‑GCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCAt                                                                                                                                                           <  1:432335/172‑1 (MQ=255)                                                                                                                                                         atTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGC‑GCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCAt                                                                                                                                                           >  2:432335/1‑172 (MQ=255)                                                                                                                                                             tCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGC‑GCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCaa                                                                         <  1:33572/250‑1 (MQ=255)                                                                                                                                                             tCCCCGATTTACGGTGACAGCCCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGGTGACCGCTTTAACCCCATTTAATGCTGC‑GCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAAACATGCCCAAGCACGCTCATTGCTGCGTGGGTGAGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGcccat                                                                        >  2:15305/1‑248 (MQ=255)                                                                                                                                                              ccccGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGC‑GCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCaat                                                                        <  1:165922/250‑1 (MQ=255)                                                                                                                                                                             tCAGACCTGAGTGGCGCTAACCATCCGGCGCATGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCAGTTAATGCTGC‑GCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTtcatc                                                         <  2:304300/250‑1 (MQ=255)                                                                                                                                                                                  ccTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGC‑GCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCagcagc                                                   >  1:78823/1‑251 (MQ=255)                                                                                                                                                                                          gcgcTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGC‑GCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGcagcag                                           <  1:484727/251‑1 (MQ=255)                                                                                                                                                                                              tAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGC‑GCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCCAGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAACTGCGACAAGCCTTTTGCCAATAATTTGTCTTTCATAAGCAGCGGCAGCAGctct                                       >  1:471439/1‑251 (MQ=255)                                                                                                                                                                                                            aggcaggcGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGC‑GCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCAc                         <  1:27768/251‑1 (MQ=255)                                                                                                                                                                                                                gggcGATTTGCAGTACTGCTGGAATTGTCACGCTATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGC‑GCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTg                     <  2:314361/250‑1 (MQ=255)                                                                                                                                                                                                                  gcGATTTGCAGTACGGCTTGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGC‑GCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGc                   <  1:153561/251‑1 (MQ=255)                                                                                                                                                                                                                   cgagttGCAGTACGGATGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCAGTTAATGCTGC‑GCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCa                  <  1:399500/247‑1 (MQ=255)                                                                                                                                                                                                                      tttGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGC‑GCCCACAAGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCg               >  1:416483/1‑251 (MQ=255)                                                                                                                                                                                                                                    ggAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCGGC‑GCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCTCCCTGGCATCGACCGCGTGCAGaa  <  2:243747/250‑1 (MQ=255)                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                      ATTAGCCAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGC‑ACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAA  >  NC_000913/10471‑10947

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