Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,514,897	0	G	A	72.5%	63.9 / 28.7	52	G46S (GGC→AGC)	ydcV	putative ABC transporter membrane subunit YdcV
Reads supporting (aligned to +/- strand):  ref base G (8/6);  new base A (14/23);  total (22/30)
Fisher's exact test for biased strand distribution p-value = 3.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.80e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GTTTATTCCCAGCAGGGGGCGATTGGCAATATGCCGATGGCGGCGGCATTCACCCTGGTGCCGATTATTCTCATCGCACTGTACCTGGCGTTCGTGAAACGTCTGGGAGCGTTCGATGCACTCTGAACGCGCACCGTTTTTCCTCAAACTGGCGGCCTGGGGCGGCGTTGTTTTCCTACATTTTCCCATCCTGATAATCGCCGCCTATGCGTTTAACACTGAAGATGCGGCGTTTAGTTTTCCACCGCAGGGCCTGACGCTGCGCTGGTTTAGCGTGGCAGCACAGCGTAGTGATATTCTTGATGCCGTGACACTGTCACTTAAAGTGGCGGCGCTGGCGACATTAATTGCGCTGGTGTTAGGGACGCTGGCAGCTGCCGCGCTGTGGCGACGAGACTTTTTCGGCAAAAACGCCATTTCGCTGTTACTGCTGCTGCCCATTGCGCTGCCGGGCATTGTCACTGGTCTGGCGTTATTAACCGCCTTTAAAACCAT  >  NC_000913/1514647‑1515141                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                     gTTTATTCCCAGCAGGGGGCGATTGGCAATATGCCGATGGCGGCGGCATTCACCCTGGTGCCGATTATTCTCATCGCACTGTACCTGGCGTTCGTGAAACGTCTGGGAGCGTTCGATGCACTCTGAACGCGCACCGTTTTTCCTCAAACTGGCGGCCTGGGGCGGCGTTGTTTTCCTACATTTTCCCATCCTGATAATCGCCGCCTATGCGTTTAACACTGAAGATGCGGCGTTTAGTTTTCCCCCGCAgg                                                                                                                                                                                                                                                      >  2:483757/1‑251 (MQ=255)     aTTCCCAGCAGGGGGCGATTGGCAATATGCCGATGGCGGCGGCATTCACCCTGGTGCCGATTATTCTCATCGCACTGTACCTGGCGTTCGTGAAACGTCTGGGAGCGTTCGATGCACTCTGAACGCGCACCGTTTTTCCTCAAACTGGCGGCCTGGGGCGGCGTTGTTTTCCTACATTTTCCCATCCTGATAATCGCCGCCTATGCGTTTAACACTGAAGATGCGGCGTTTAGTTTTCCACCGCAGAGCCt                                                                                                                                                                                                                                                  >  1:109158/1‑251 (MQ=255)       tCCCAGCAGGGGGCGATTGGCAATATGCCGATGGCGGCGGCATTCACCCTGGTGCCGATTATTCTCATCGCACTGTACCTGGCGTTCGTGAAACGTCTGGGAGCGTTCGATGCACTCTGAACGCGCACCGTTTTTCCTCAAACTGGCGGCCTGGGGCGGCGTTGTTTTCCTACATTTTCCCATCCTGATAATCGCCGCCTATGCGTTTAACACTGAAGATGCGGCGTTTAGTTTTCCCCCGCCGGGCCTg                                                                                                                                                                                                                                                 >  2:503992/1‑250 (MQ=255)        cccAGCAGGGGGCGATTGGCAATATGCCGATGGCGGCGGCATTCACCCTGGTGCCGATTATTCTCATCGCACTGTACCTGGCGTTCGTGAAACGTCTGGGAGCGTTCGATGCACTCTGAACGCGCACCGTTTTTCCTCAAACTGGCGGCCTGGGGCGGCGTTGTTTTCCTACATTTTCCCATCCTGATAATCGCCGCCTATGCGTTTAACACTGAAGATGCGGCGTTTAGTTTTCCACCGCAGAGCCTGa                                                                                                                                                                                                                                                <  1:467728/250‑1 (MQ=255)          cagcagGGGGCGATTGGCAATATGCCGATGGCGGCGGCATTCACCCTGGTGCCGATTATTCTCATCGCACTGTACCTGGCGTTCGTGAAACGTCTGGGAGCGTTCGATGCACTCTGAACGCGCACCGTTTTTCCTCAAACTGGCGGCCTGGGGCGGCGTTGTTTTCCTACATTTTCCCATCCTGATAATCGCCGCCTATGCGTTTAACACTGAAGATGCGGCGTTTAGTTTTCCACCGCAGGGCCTGAcgc                                                                                                                                                                                                                                             >  1:109517/1‑251 (MQ=255)           agcagGGGGCGATTGGCAATAGGCCGATGGCGGCGGCATTCACCCTGGTGCCGATTATTCTCATCGCACTGTACCTGGCGTTCGTGAAACGTCTGGGAGCGTTCGATGCACTCGGAACGCGCACCGTTTTTCCTCAAACTGGCGGCCTGGGGCGGCGTTGTTTTCCTACATTTTCCCATCCTGATAATCGCCGCCTATGCGTTTAACACTGAAGATGCGGCGTTGAGTTTTCCACCGCAGGTCCTGAcgct                                                                                                                                                                                                                                            <  2:421903/251‑1 (MQ=255)                   gCGATTGGCAATATGCCGATGGCGGCGGCATTCACCCTGGTGCCGATTATTCTCATCGCACTGTACCTGGCGTTCGTGAAACGTCTGGGAGCGTTCGATGCACTCTGAACGCGCACCGTTTTTCCTCAAACTGGCGGCCTGGGGCGGCGTTGTTTTCCTACATTTTCCCATCCTGATAATCGCCGCCTATGCGTTTAACACTGAAGATGCGGCGTTTAGTTTTCCACCGCAGAGCCTGACGCTGCGCTGGt                                                                                                                                                                                                                                    <  1:589647/251‑1 (MQ=255)                    cGATTGGCAATATGCCGATGGCGGCGGCATTCACCCTGGTGCCGATTATTCTCATCGCACTGTACCTGGCGTTCGTGAAACGTCTGGGAGCGTTCGATGCACTCTGAACGCGCACCGTTTTTCCTCAAACTGGCGGCCTGGGGCGGCGTTGTTTTCCTACATTTTCCCATCCTGATAATCGCCGCCTATGCGTTTAACACTGAAGATGCGGCGTTTAGTTTTCCACCGCAGAGCCTGACGCTGCGCTGGtt                                                                                                                                                                                                                                   >  1:549981/1‑251 (MQ=255)                      aTTGGCAATATGCCGATGGCGGCGGCATTCACCCTGGTGCCGATTATTCTCATCGCACTGTACCTGGCGTTCGTGCAACGTCTGGGAGCGTTCGATGCACTCTGCACGCGCACCGTTTTTCCTCAAACTGGAGGCCTGGGGCGGCGTTGTTTTCCTACATTTTCCCATCCTGATAATAGCCGCCTATGCGTTTAAAACTGAAGATGCGGCGTTTAGTTTTCCACCGCAGAGCCTCACGCTGCGCTGGttt                                                                                                                                                                                                                                  >  2:311897/1‑250 (MQ=255)                        tGGCAATAGGCCGATGGCGGCGGCATTCACCCTGGTGCCGATTATTCTCATCGCACTGTACCTGGCGTTCGTGAAACGTCTGGGAGCGTTCGATGCACTCTGAACGCGCACCGTTTTTCCTCAAACTGGCGGCCTGGGGCGGCGTTGTTTTCCTACATTTTCCCATCCTGATAATCGCCGCCTATGCGTTTAACACGGAAGATGCGGCGTTTAGTTTTCCACCGCAGAGCCTGACGCTGCGCTGGTTTAg                                                                                                                                                                                                                                <  1:161995/250‑1 (MQ=255)                            aaTATGCCGATGGCGGCGGCATTCACCCTGGTGCCGATTATTCTCATCGCACTGTACCTGGCGTTCGTGAAACGTCTGGGAGCGTTCGATGCACTCTGAACGCGCACCGTTTTTCCTCAAACTGGCGGCCTGGGGCGGCGTTGTTTTCCTACATTTTCCCATCCTGATAATCGCCGCCTATGCGTTTAACACTGAAGATGCGGCGTTTAGTTTTCCACCGCAGAGCCTGACGCTGCGCTGGTTTAGCGTgg                                                                                                                                                                                                                           <  1:234600/251‑1 (MQ=255)                               atGCTGATGGCGGCGGCACTCACCCTGGTGCCGATTATTCTCATCGCACTGTCCCTGGCGTTCGTGAAACGTCTGGTAGCGTGCGATGCACTCTGAACGCGCACCTTTTTTCCTCAAACTGGCGGCCTGGGGCGGCGTTGTTTTCCTACAGTTTCCCATCCTGATAATCGCCGCCTATGCGTTTAACACTGAAGATGCGGCGTTTAGTTTTCCACCGCAGAGCCTGACGCTGCGCTGGTTTAGCGTGgcag                                                                                                                                                                                                                        <  2:382305/251‑1 (MQ=255)                               atGCCGATGGCGGCGGCATTCACCCTGGTGCCGATTATTCTCATCGCACTGTACCTGGCGTTCGTGAAACGTCTGGGAGCGTTCGATGCACTCTGAACGCGCACCGTTTTTCCTCAAACTGGCGGCCTGGGGCGGCGTTGTTTTCCTACATTTTCCCATCCTGATAATCGCCGCCTATGCGTTTAACACTGAAGATGCGGCGTTTAGTTTTCCACCGCAGAGCCTGACGCTGCGCTGGTTTAGCGTGgca                                                                                                                                                                                                                         <  2:543672/250‑1 (MQ=255)                                             ggcATTCACCATGGTGCCGATTATTCTCATCGCACTGTACCTGGCGTTCGTGAAACGTCTGGGAGCGTTCGATGCACTCTGAACGCGCACCGTTTTTCCTCAAACTGGCGGCCTGGGGCGGCGTTGTTTTCCTACATTTTCCCATCCTGATAATCGCCGCCTATGCGTTTAACACTGAAGATGCGGCGTTTAGTTTTCCCCCGCAGAGCCTGACGCTGCGCTGGTTTAGCGTGGCAGCCCAGCGCAGTggt                                                                                                                                                                                                          >  2:514169/1‑249 (MQ=255)                                                   cACCCTGGGGCCGATTATTCTCAGCGCACTGTACCTGGCGTTCGTGAAACGTCTGGGAGCGTTCGATGCATTCTGAACGCGCACCGTTTTGCCTCAAACTGGCGGCCTGGGGCGTCGTTGTTTTCCTACATGTTCCCATCCTGATAATCGCCGCCTATGTGGTTAACACGGAAGATGCGGCGGTTAGGTTTCCACCGCAGAGCCTGACGCTGCGCTGGTCTAGCGTGTCAGCACAGCGTAGTGATATTCt                                                                                                                                                                                                     <  2:259279/250‑1 (MQ=255)                                                     cccTGGTGCCGATTATTCTCATCGCACTGTACCTGGCGTTCGTGAAACGTCTGGGAGCGTTCGATGCACTCTGAACGCGCACCGTTTTTCCTCAAACTGGCGGCCTGGGGCGGCGTTGTTTTCCTACATTTTCCCATCCTGATAATCGCCGCCTATGCGTTTAACACTGAAGATGCGGCGTTTAGTTTTCCACCGCAGAGCCTGACGCTGCGCTGGTTTAGCGTGGCAGCACAGCGTAGTGATATTCTTGa                                                                                                                                                                                                  <  1:137880/251‑1 (MQ=255)                                                          gTGCCGATTATTCTCATCGCACTGTACCTGGCGTTCGTGAAACGTCTGGGAGCGTTCGATGCACTCTGAACGCGCACCGTTTTTCCTCAAACTGGCGGCCTGGGGCGGCGTTGTTTTCCTACATTTTCCCATCCTGATAATCGCCGCCTATGCGTTTAACACTGAAGATGCGGCGTTTAGTTTTCCACCGCAGAGCCTGACGCTGCTCTGGTTTAGCGTGGCAGCACAGCGTAGTGATATTCTTGATGCCg                                                                                                                                                                                             <  1:512525/251‑1 (MQ=255)                                                            gCCGATTATTCTCATCGCACTGTACCTGGCGTTCGTGAAACGTCTGGGAGCGTTCGATGCACTCTGAACGCGCACCGTTTTTCCTCAAACTGGCGGCCTGGGGCGGCGTTGTTTTCCTACATTTTCCCATCCTGATAATCGCCGCCTATGCGTTTAACACTGAAGATGCGGCGTTTAGTTTTCCACCGCAGAGCCTGACGCTGCGCTGGTTTAGCGTGGCAGCACAGCGTAGTGATATTCTTGATGCCGt                                                                                                                                                                                            <  2:438237/250‑1 (MQ=255)                                                                                       ggCGTTCGTGAAACGTCTGGGAGCGTTCGATGCACTCTGAACGCGCACCGTTTTTCCTCAAACTGGCGGCCTGGGGCGGCGTTGTTTTCCTACATTTTCCCATCCTGATAATCGCCGCCTATGCGTTTAACACTGAAGATGCGGCGTTTAGTTTTCCACCGCAGGGCCTGACGCTGCGCTGGTTTAGCGTGGCAGCACAGCGTAGTGATATTCTTGATGCCGTGACACTGTCACTTAAAGTGGCGGCGCTg                                                                                                                                                                >  1:625053/1‑251 (MQ=255)                                                                                       ggCGTTCGTGAAACGTCTGGGAGCGTTCGATGCACTCTGAACGCGCACCGTTTTTCCTCAAACTGGCGGCCTGGGGCGGCGTTGTTTTCCTACATTTTCCCATCCTGATAATCGCCGCCTATGCGTTTAACACTGAAGATGCGGCGTTTAGTTTTCCACCGCAGAGCCTGACGCTGCGCTGGTTTAGCGTGGCAGCACAGCGTAGTGATATTCTTGATGCCGTGACACTGTCACTTAAAGTGGCGGCGCt                                                                                                                                                                 <  1:163522/250‑1 (MQ=255)                                                                                          gTTCGTGAAACGTCTGGGAGCGTTCGATGCACTCTGAACGCGCACCGTTTTTCCTCAAACTGGCGGCCTGGGGCGGCGTTGTTTTCCTACATTTTCCCATCCTGATAATCGCCGCCTATGCGTTTAACACTGAAGATGCGGCGTTTAGTTTTCCACCGCAGAGCCTGACGCTGCGCTGGTTTAGCGTGGCAGCACAGCGTAGTGATATTCTTGATGCCGTGACACTGTCACTTAAAGGGGCGGGGCTGGCg                                                                                                                                                             >  2:337844/1‑251 (MQ=255)                                                                                             cGTGAAACGTCGGGGAGCGTTCGATGCGCTCTGAACGCGCACCGTTTTTCCTCAAACGGGCGGCCTGGGGCGGCGTTGTTTGCCTACATTTTCCCAGCCTGATAATCGCCGCCTATGCGTTTAACACTGAAGATGCGGCGTTTAGTTGTCCACCGCAGGGCCTGACGCTGCGCTGGTTTAGCGTGGCAGCACAGCGTAGTGATATTCTTGATGCCGTGACACTGTCACTTAAAGTGGCGGCGCTGGCGACa                                                                                                                                                          <  1:84917/251‑1 (MQ=255)                                                                                                  aaCGTCTGGGAGCGTTCGATGCACTCTGAACGCGCCCCGTTTTTCCTCAAACTGGCGGCCTGGTGCGGCGTTGTTTTCCTACATTTTCCCATCCTGATAATCGCCGCCTATGCGTTTAACACTGAAGATGCGGCGTTTAGTTTTCCACCGCAGAGCCTTACGCTGCGCTGGTTTAGCGTGGCAGCACAGCGTAGTGATATTCTTGATGCCGTGACACTGTCACTTAAAGTGGCGGCGCTGGCGACATTAAt                                                                                                                                                     <  2:611578/251‑1 (MQ=255)                                                                                                  aaCGTCGGGGAGCGTTCGATGCACTCTGAACGCGCACCGTTTTTCCTCAAACTGGCGGCCTGGGGCGGCGTTGTTTTCCTACATTTTCCCATCCTGATAATCGCCGCCTATGCGTTTAACACTGAAGATGCGGCGTTTAGTTTTCCACCGCAGAGCCTGACGCTGCGCTGGTTTAGCGTGGCAGCACAGCGTAGTGATATTCTTGATGCCGTGACACTGTCACTTAAAGTGGCGGCGCTGGCGACATTAAt                                                                                                                                                     <  2:299640/251‑1 (MQ=255)                                                                                                       cTGGGAGCGTTCGATGCACTCTGAACGCGCACCGTTTTTCCTCAAACTGGCGGCCTGGGGCGGCGTTGTTTTCCTACATTTTCCCATCCTGATAATCGCCGCCTATGCGTTTAACACTGAAGATGCGGCGTTTAGTTTTCCACCGCAGAGCCTGACGCTGCGCTGGTTTAGCGTGGCAGCACAGCGTAGTGATATTCTTGATGCCGTGACACTGTCACTTAAAGTGGCGGCGCTGGCGACCTTAATTtcgc                                                                                                                                                >  2:191553/1‑251 (MQ=255)                                                                                                           gAGCGTGCGATGCACGCGGAACGCGCACCGTGTTTCCGCAAACTGGCGGCCTGGGGCGGCTTTGTTTTCCTACATTTTCCCATCCTGATAAGCGCCGCCTATGCGTTTAACACTTAAGATGCTGCGTTTAGGTTTCCACCGCAGAGCCTGACGCTGCGCTGGTGTAGCGTGGCAGCACAGCGTAGTGATATTCTTGATGCCGGTACTCTGTCACTTAAAGTGGCGGCGCTTGCGACATTAATTGCGCttgt                                                                                                                                            <  2:318106/251‑4 (MQ=255)                                                                                                                 tCGATGCACTCTGAACGCGCACCGTTTTTCCTCAAACTGGCGGCCTGGGGCGGCGTTGTTTTCCTCCATTTTCCCATCCTGAGAATCGCCGCCTATGCGTTTAACACTGAAGATGCGGCGTTTAGTTTTCCACCGCAGAGCCTGACGCTGCGCTGGTTTAGCGTGGCAGCACAGCGTAGTGATATGCTTGATGCCGTGACACTGTCACTTAAAGTGGCGGCGCTGGCGACATTAATTGCGCTGGTGTTAg                                                                                                                                       <  2:40850/250‑1 (MQ=255)                                                                                                                     cgCACTCTTAACGCGCACCGTTTTTCCTCAACCTGGCGGCCTGGGGCGGCGTGGTTTTCCTACATTTTCCCAGCCTGATAATCGCCGCCTATGCGTTTAACACTGAAGATGCGGCGTTTAGTTTCCCACCGCAGTGCCTGCCGCTGCGCTGGTTTAGCGTGGCAGCACAGCGTAGTGATATTCTTGATGCCGTGACACTGTCACTTAAAGTGGCGGCGCTTGCGACATGAATGGCGCCGGTGTTAGGGACg                                                                                                                                  <  2:109294/250‑1 (MQ=255)                                                                                                                                        gTTTTTCCTCAAACTGGCGGCCGGGGGCGGCGTTGTTTTCCTACATTTTCCCATCCTGATAATCGCCGCCTATGCGTTTAACACTGAAGATGCGGCGTTTAGTTTGCCACCGCAGGGCCTGACGCTGCGCTGGTTTAGCGTGGCAGCACAGCGTAGTGATATTCTTGATGCCGTGACACTGTCACTTAAAGTGGCGGCGCTGGCGACATTAATTGCGCTGGGGTTAGGGACGCTGGCAGCTGCCGCGCtgt                                                                                                               <  2:327286/251‑1 (MQ=255)                                                                                                                                             tCCTCAAACTGGCGGCCTGGGGCGGCGTTGTTTTCCTACATTTTCCCATCCTGATAATCGCCGCCTATGCGTTTAACACTGAAGATGCGGCGTTTAGTTTTCCACCGCAGAGCCTGACGCTGCGCTGGTTTAGCGTGGCAGCACAGCGTAGTGATATTCTTGATGCCGTGACACTGTCACTTAAAGTGGCGGCGCTGGCGACATTAATTGCGCTGGTGTTAGGGACGCTGGCAGCTGCCGCGCTGTGGcg                                                                                                           <  2:151410/250‑1 (MQ=255)                                                                                                                                              ccTCAAACTGGCGGCCTGGGGCGGCGTTGTTTTCCTACATTTTCCCATCCTGATAATCGCCGCCTATGCGTTTAACACTGAAGATGCGGCGTTTAGTTTTCCACCGCAGGGCCTGACGCTGCGCTGGTTTAGCGTGGCAGCACAGCGTAGTGATATTCTTGATGCCGTGACACTGTCACTTAAAGTGGCGGCGCTGGCGACATTAATTGCGCTGGTGTTAGGGACGCTGGCAGCTGCCGCGCTGTGGcgac                                                                                                         >  1:479057/1‑251 (MQ=255)                                                                                                                                                tCAAACTGGCGGCCTGGGGCGGCGTTGTTTTCCTACATTTTCCCATCCTGATAATCGCCGCCTATGCGTTTAACACTGAAGATGCGGCGTTTAGTTTTCCACCGCAGGGCCTGACGCTGCGCTGGTTTAGCGTGGCAGCACAGCGTAGTGATATTCTTGATGCCGTGACACTGTCACTTACAGTGGCGGCGCTGGCGACCTTTATTCGGCTTGGGTTAGGGAAGCTGGCAGATTCCGCGCTGTGGcgcgg                                                                                                        >  2:586547/1‑247 (MQ=255)                                                                                                                                                  aaaCTGGCGGCCTGGGGCGGCGTTGTTTTCCTACATTTTCCCATCCTGATAATCGCCGCCTATGCGTTTAACACTGAAGATGCGGCGTTTAGTTTTCCACCGCAGAGCCTGACGCTGCGCTGGTTTAGCGTGGCAGCACAGCGTAGTGATATTCTTGATGCCGTGACACTGTCACTTAAAGTGGCGGCGCTGGCGACATTAATTGCGCTGGTGTTAGGGACGCTGGCAGCTGCCGCGCTGTGGCGACgaga                                                                                                     <  1:575748/251‑1 (MQ=255)                                                                                                                                                    aCTGGCGGCCTGGGGCGGCGTTGTTTTCCTACATTTTCCCATCCTGATAATCGCCGCCTATGCGTTTAACACTGAAGATGCGGCGTTTAGTTTTCCCCCGCAGAGCCTGACGCTGCGCGGGTTTAGCGTGGCAGCACAGCGTAGTGATATTCTTGAGGCCGTGACACTGTCACTTAAAGTGGCGGCGCTGGCGACATTAATTGCGCTGGTGTTAGGGACGCTGGCAGCTGCCTCGCTGTGGCGACGAGACt                                                                                                   <  2:602642/251‑1 (MQ=255)                                                                                                                                                             cTGGGGCGGCGTTGTTTTCCTACATTTTCCCATCCTGATAATCGCCGCCTATGCGTTTAACACTGAAGATGCGGCGTTTAGTTTTCCACCGCAGGGCCTGACGCTGCGCTGGTTTAGCGTGGCAGCACAGCGTAGTGATATTCTTGATGCCGTGACACTGTCACTTAAAGTGGCGGCGCTGGCGACATTAATTGCGCTGGTGTTAGGGACGCTGGCAGCTGCCGCGCTGTGGCGAAGAGACTTTTTCGGc                                                                                           >  1:222849/1‑250 (MQ=255)                                                                                                                                                              tGGGGCGGCGTTGTTTTCCTACATTTTCCCATCCTGATAATCGCCGCCTATGCGTTTAACACTGAAGATGCGGAGTTTAGTTTTCCACCGCAGAGCCTGACGCTGCGCTGGTTTAGCGTGGAAGCACAGCGTAGTGATATTCTTGATGCCGTGACACTGTCACTTAAAGTGGCGGCGCTGGCGGCATTAATTGCGCTGGTGTTTGGGACGCTGGCAGCTGCCGCGCTGTGGCGCCGGGACTTTTTTGGCa                                                                                          >  1:309490/1‑250 (MQ=255)                                                                                                                                                                   cggcgTTGTTTTCCTACATTTTCCCATCCTGATAATCGCCGCCTATGCGTTTAACACTGAAGATGCGGCGTTTAGTTTTCCACCGCAGAGCCTGACGCTGCGCTGGTTTAGCGTGGCAGCACAGCGTAGTGATATTCTTGATGCCGTGACACTGTCACTTAAAGTGGCGGCGCTGGCGAAATTAATTGCGCTGGTGTTAGGGACGCTGGCAGCTGCCGCGCCGGGGCGACGAGACTTTTTCGGCAAAAACg                                                                                    >  2:384184/1‑251 (MQ=255)                                                                                                                                                                     gcgTTGTTTTCCTACATTTTCCCATCCTGATAATCGCCGCCTATGCGTTTAACACTGAAGATGCGGCGTTTAGTTTTCCACCGCAGAGCCTGACGCTGCGCTGGTTTAGCGTGGCAGCACAGCGTAGTGATATTCTTGATGCCGTGACACTGTCACTTAAAGTGGCGGCGCTGGCGACATTAATTGCGCTGGTGTTAGGGACGCTGGCAGCTGCCGCGCTGTGGCGGCGAGACTTTTTTGGCAAAAACGcc                                                                                  >  2:65765/1‑251 (MQ=255)                                                                                                                                                                         tgttTTCATACATTATCCCCTCCTGATAATCGCCCCCTCTGCGTTTAACACTGAAGATCCGGCGTTTAGTTTTCCACCGAAGGGCCTGTCGCTGCGCTGGTTTAGCGTGGCAGCACAGCGTAGTGATATTCTTGATGCCGTGACACTGTCACTTAAAGTGGCGGCGCTGGCGACATTAATTGCGCTGGTGTTAGGGACGCCGGCAGCTGCCGCCCTGTGGCGACTGGACTTTTTCGGCACAAACGCCAttt                                                                              >  2:59631/1‑251 (MQ=255)                                                                                                                                                                            tttCCTACATTTTCCCATCCTGATAATCGCCGCCTATGCGTTTAACACTGAAGATGCGGCGTTTAGTTTTCCACCGCAGGGCCTGACGCTGCGCTGGTTTAGCGTGGCAGCACAGCGTAGTGATATTCTTGATGCCGTGACACTGTCACTTAAAGTGGCGGCGCTGGCGACATTAATTGCGCTGGTGTTAGGGACGCTGGCAGCTGCCGCGCTGTGGCGACGAGACTTTTTCGGCAAAAACGCCATTTCg                                                                            <  1:311750/250‑1 (MQ=255)                                                                                                                                                                            tttCCTACATTTTCCCATCCTGATAATCGCCGCCTATGCGTTTAACACTGAAGATGCGGCGTTTAGTTTTCCACCGCAGAGCCTGACGCTGCGCTGGTTTAGCGTGGCAGCACAGCGTAGTGATATTCTTGATGCCGTGACACTGTCACTTAAAGTGGCGGCGCTGGCGACATTAATTGCGCTGGTGTTAGGGACGCTGGCAGCTGCCGCGCTGTGGCGACGAGACTTTTTCGGCAAAAACGCCATTTCg                                                                            <  1:589011/250‑1 (MQ=255)                                                                                                                                                                              tCCTACATTTTCCCATCCTGATAATCGCCGCCTATGCGTTTAACACTGAAGATGCGGCGTTTAGTTTTCCACCGCAGAGCCTGACGCTGCGCTGGTTTAGCGTGGCAGCACAGCGTAGTGATATTCTTGATGCCGTGACACTGTCACTTAAAGTGGCGGCGCTGGCGACATTAATTGCGCTGGTGTTAGGGACGCTGGCAGCTGCCGCGCTGTGGCGACGAGACTTTTTCGGCAAAAACGCCATTTCGCt                                                                          <  1:30304/250‑1 (MQ=255)                                                                                                                                                                              tCCTACATTTTCCCAGCCTGATAAGCGCCGCCTATGCGTTTAACACGGAAGATGCGGCGTGTAGTTTGCCACCGCAGGGCCTGACGCTGCGCTGGTTTAGCGTGGCAGCACAGCGTAGTGATATTCTTGATGCCGTGACACTGTCACTTAAAGTGGCGGCGCTGGCGACATTAATTGCGCTGGTGTTAGGGACGCTGGCAGCTGCCGCGCTGTGGCGACGAGACTTTTTCGGCAAAAACGCCATTTCGCt                                                                          <  2:626237/250‑1 (MQ=255)                                                                                                                                                                                    aTTTTCCCATCCTGATAATCGCCGCCTATGCGTTTAACACTGAAGATGCGGCGTTTAGTTTTCCACCGCAGGGCCTGACGCTGCGCTGGTTTAGCGTGGCAGCACAGCGTAGTGATATTCTTGATGCCGTGACACTGTCACTTAAAGTGGCGGCGCTGGCGACATTAATTGCGCTGGTGTTAGGGACGCTGGCAGCTTCCGCGCTGTGGCGACGAGACTTTTTCGGCAAAAACGCCATTTCGCTGTTActg                                                                   <  2:222849/251‑1 (MQ=255)                                                                                                                                                                                          ccATCCTGATAATCGCCGCCTATGCGTTTAACACTGAAGATGCGGCGTTTAGTTTTCCACCGCAGAGCCTGACGCTGCGCTGGTTTAGCGTGGCAGCACAGCGTAGTGATATTCTTGATGCCGTGACACTGTCACTTAAAGTGGCGGCGCTGGCGACATTAATTGCGCTGGTGTTAGGGACGCTGGCAGCTGCCGCGCTGTGGCGACGAGACTTTTTCGGCAAAAACGCCATTTCGCTGTTActgctgctg                                                             >  1:195526/1‑251 (MQ=255)                                                                                                                                                                                                  aTAATCGCCGCCTATGCAGTTAACACTGAAGATGCGGCGTTTAGTTTTCCACCGCAGAGCCTGACGCTGCGCTGGTTTAGCGTGGCAGCACAGCGAAGTGATATTCTTGATGCCGTGACACTGTCACTTAAAGAGGCGGCGCTTGCGACATTAATTTCGCTGGTGTTTGGGACGCTGGCAGCTGCCGCGCTGTGGCGACGGGACTTTTTCTGCAAAAACGCCATTTCGCTGTTACTGCTGCTGCCCCTTgc                                                     >  2:416817/1‑251 (MQ=255)                                                                                                                                                                                                               aTGCGTTTAACACTGAAGATGCGGCGTTTAGTTTTCCACCGCAGAGCCTGACGCTGCGCTGGTTTAGCGTGGCAGCACAGCGTAGTGATATTCTTGATGCCGTGACACTGTCACTTAAAGTGGCGGCGCTGGCGACATTAATTGCGCTGGTGTTAGGGACGCTGGCAGCTGCCGCGCTGTGGCGACGAGACTTTTTCGGCAAAAACGCCATTTCGCTGTTACTGCTGCTGCCCATTTCGCTGCCGGGcctt                                        >  1:181228/1‑248 (MQ=255)                                                                                                                                                                                                                      tAACACTGAATATGCGGCGTTTAGTTTTCCACCGCAGAGCCTGACGCTGCGCTGGTTTAGCGTGGCAGCACAGCGTAGTGATATTCTTGATGCCGTGACACTGTCACTTAAAGTGGCGGCGCTGGCGACATTAATTGCGCTGGTGTTAGGGACGCTGGCAGCTGCCGCGCTGTGGCGACGAGACTTTTTCGGCAAAAACGCCATTTCGCTGTTACTGCTGCTGCCCATTGCGCTGCCGGGCATTGTCACt                                  <  2:459294/250‑1 (MQ=255)                                                                                                                                                                                                                         cacTGAAGATGCGGCGTTTAGTTTTCCACCGCCGAGACTGACGCTGCGCTGGTTTAGCGTGGCAGCACAGCGTAGTGCTATTCTTGATGCCGTGACACTGTCACTTAAAGCGGCGGCGCTGGCGACATTTATTGCGCTGGTGTTAGGGACGCTGGCAGGTGCCGCGCTGTGGCGACGAGCCTTTTTCGGCAAAAACGCCATTTCGCTGCTACTGCTGCTGCCCCATGCGCTGCCGGGCCTTGTCACCGGTc                              >  2:34249/1‑251 (MQ=255)                                                                                                                                                                                                                             gAAGATGCGGCGTTTAGTTTTCCACCGCAGAGCCTGACGCTGCGCTGGTTTAGCGTGGCAGCACAGCGTAGTGATATTCTTGATGCCGTGACACTGTCACTTAAAGTGGCGGCGCTGGCGACATTAATTGCGCTGGTGTTAGGGACGCTGGCAGCTGCCGCGCTGTGGCGACGAGACTTTTTCGGCAAAAACGCCATTTCGCTGTTACTGCTGCTGCCCATTGCGCTGCCGGGCATTGTCACTGGTCTGGc                          <  2:597526/251‑1 (MQ=255)                                                                                                                                                                                                                                 aTGCGGCGTTTAGGTTTCCACCGCAGAGCCTGACGCTGCGCTGGTTTAGCGTGGCAGCACAGCGTAGTGATATTCTTGATGCCGTGACACTGTCACTTAAAGTGGCGGCGCTGGCGACATTAATTGCGCTGGTGTTAGGGACGCTGGCAGCTGCCGCGCTGTGGCGACGAGACGTTTTCGGCAAAAACGCCATTTCGCTGTTACTGCTGCTGCCCATGGCGGTGCCGGGCATTGTCACTGGTCTGGCGtta                      <  2:291913/251‑1 (MQ=255)                                                                                                                                                                                                                                  tGCGGCGTTTAGTTTTCCACCGCAGGGCCTGACGCTGCGCTGGTTTAGCGTGGCAGCACAGCGTAGTGATATTCTTGATGCCGTGACACTGTCACTTAAAGTGGCGGCGCTGGCGACATTAATTGCGCTGGTGTTAGGGACGCTGGCAGCTGCCGCGCTGTGGCGACGAGACTTTTTCGGCAAAAACGCCATTTCGCTGTTACTGCTGCTGCCCATTGCGCTGCCGGGCATTGTCACTGGTCTGGCGttat                     >  1:433622/1‑251 (MQ=255)                                                                                                                                                                                                                                                     ccGCAGAGCCTGACGCTGCGCTGGTTTAGCGTGGCAGCACAGCGTAGTGATATTCTTGATGCCGTGAAACTGTCACTTAAAGTGGCGGCGCTGGCGACATTAATTGCGCTGGTGTTAGGGACGCTGGCAGCTGCCGCGCTGTGGCGACGAGCCTTTTTCGGCAAAAACGCCATTTCGCTGTTACTGCTGCTGCCCCTTGCGCTGCCGGGCATTGTCACTGGTCTGGCGTTATTAACCGCCTTTAAAACCAt  >  1:253743/1‑251 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                     GTTTATTCCCAGCAGGGGGCGATTGGCAATATGCCGATGGCGGCGGCATTCACCCTGGTGCCGATTATTCTCATCGCACTGTACCTGGCGTTCGTGAAACGTCTGGGAGCGTTCGATGCACTCTGAACGCGCACCGTTTTTCCTCAAACTGGCGGCCTGGGGCGGCGTTGTTTTCCTACATTTTCCCATCCTGATAATCGCCGCCTATGCGTTTAACACTGAAGATGCGGCGTTTAGTTTTCCACCGCAGGGCCTGACGCTGCGCTGGTTTAGCGTGGCAGCACAGCGTAGTGATATTCTTGATGCCGTGACACTGTCACTTAAAGTGGCGGCGCTGGCGACATTAATTGCGCTGGTGTTAGGGACGCTGGCAGCTGCCGCGCTGTGGCGACGAGACTTTTTCGGCAAAAACGCCATTTCGCTGTTACTGCTGCTGCCCATTGCGCTGCCGGGCATTGTCACTGGTCTGGCGTTATTAACCGCCTTTAAAACCAT  >  NC_000913/1514647‑1515141

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