Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	40,525	0	C	T	37.2%	103.7 / 29.5	81	L469L (CTG→CTA)	caiT	putative transporter
Reads supporting (aligned to +/- strand):  ref base C (12/37);  new base T (29/0);  total (44/37)
Fisher's exact test for biased strand distribution p-value = 6.19e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAATCGGTTTCAGGCCGCCG‑AGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGC  >  NC_000913/40379‑40669                                                                                                                                                     |                                                                                                                                                  aaaTCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTTAAAAGGTTTTAgg                                                                                                                                                >  2:632339/1‑150 (MQ=255) aaaTCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑GTTTGATTGGTTTTAGgc                                                                                                                                               >  1:420653/1‑151 (MQ=255)  aaTCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTTAATCGGGTTTTGgcc                                                                                                                                              >  2:358694/1‑151 (MQ=255)    tCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑GTTGAATTGGGTTTAGgccgc                                                                                                                                            >  2:1349839/1‑151 (MQ=255)         cTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTTAATCGGGTTTAGGGCGCCC‑Ag                                                                                                                                       >  2:1377776/1‑150 (MQ=255)           ttCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAATCGGTTTCAGGCCGCCG‑Agc                                                                                                                                      <  2:495629/149‑1 (MQ=255)            tCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGTAACGGTTTTATGCCGCCG‑CACGcc                                                                                                                                   >  2:1331456/1‑151 (MQ=255)              aaCCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GGTTTTGAATGGGTTTCA‑GCCGCCG‑GGCGcc                                                                                                                                   >  2:86859/1‑149 (MQ=255)              aaCCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAATCGGTTTCAGGCCGCCG‑AGCGCcag                                                                                                                                 <  1:1377684/151‑1 (MQ=255)                ccTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTTGATTGGGTTTAGgccgcccggcggcc                                                                                                                                   >  2:529511/1‑140 (MQ=255)                 ctctTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGG‑TTTTATCCGGTTTAAGCcccccgccgcc                                                                                                                                     >  1:951251/1‑141 (MQ=255)                 ctctTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTTAATCGGGTTCAGGCCGCCG‑CGCGcc                                                                                                                                   >  1:1036232/1‑146 (MQ=255)                     tGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTTAATAGGGTTTAGGGCGCCC‑GGCGGCCGCAGaa                                                                                                                            >  2:547353/1‑149 (MQ=255)                     tGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAATCGGTTTCAGGCCGCCG‑AGCGCCAGCAGaa                                                                                                                            <  1:103897/149‑1 (MQ=255)                     tGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAATCGGTTTCAGGCCGCCG‑AGCGCCAGCAGaa                                                                                                                            <  1:444674/149‑1 (MQ=255)                      gATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTTAATCGGTTTCAGGCCGCCC‑CGCGGCCGCAGAAc                                                                                                                           >  2:1284643/1‑149 (MQ=255)                         atTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAATCGGTTTTAGGGCGCCG‑AGCGCCAGCCGAACaaaa                                                                                                                       >  1:914594/1‑148 (MQ=255)                           tttGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAATCGGTTTCAGGCCGCCG‑AGCGCCAGCAGAACAATAcc                                                                                                                     <  2:226120/150‑1 (MQ=255)                            ttGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑GTTTAATCGGGTTTTGGTCGCCG‑AGCGCCAGCAGAACaaaaacg                                                                                                                    >  1:1030431/1‑145 (MQ=255)                                gggTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTTGATCGGTTTTAGGCCGCCC‑CGCGCCAGCAGAAAAAAAcccaaaa                                                                                                                >  2:183845/1‑145 (MQ=255)                                gggTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTTAAATGGGTTCAGGCCGCCG‑AGGGCCCGCAGAACaaaaaacaaaa                                                                                                                >  2:810405/1‑141 (MQ=255)                                gggTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAATCGGTTTCAGGCCGCCG‑AGCGCTAGCAGAACAATACCGATAAt                                                                                                               <  1:287724/151‑1 (MQ=255)                                  gTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAATCGGTTTCAGGCCGCCC‑CGCGCCAGCAGAACAAAAACGAAAAAGc                                                                                                             >  2:164209/1‑151 (MQ=255)                                      attaatCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTTAATTGGGTTTAGGCCGCCG‑CGCGCCCGCAGAACAATACCGATaa                                                                                                                >  1:330107/1‑144 (MQ=255)                                      attaatCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAATCGGTTTCAGGCCGCCG‑AGCGCCAGCAGAACAATACCGATAAt                                                                                                               <  2:330119/145‑1 (MQ=255)                                            cTTTCCAGTTCTGGTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GGCTTGGAATGGGTTTAATG‑CGCCGCAGCGCCCGC‑GAAAAATACAaataatattcaacaaa                                                                                                     >  1:271642/1‑139 (MQ=255)                                             tttCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTTAATCGGGGTTAGGGCGCCG‑AGCGCCCGCAGAAAAAAAACGATAAAACGAACCaaa                                                                                                     >  1:582174/1‑146 (MQ=255)                                              ttCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTTAATCGGGTTTTAGCCGCCC‑AGCGCCCGCCGAACAAAAACGATAATGCCAAcc                                                                                                        >  2:1247600/1‑144 (MQ=255)                                              ttCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAATCGGTTTCAGGCCGCCG‑AGCGCCAGCAGAACAATACCGATAATGCCAACCAGaa                                                                                                    <  2:380627/148‑1 (MQ=255)                                              ttCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAATCGGTTTCAGGCCGCCG‑AGCGCCAGCAGAACAATACCGATAATGCCAACCAGAAtt                                                                                                  <  1:1363027/150‑1 (MQ=255)                                              ttCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAATCGGTTTCAGGCCGCCG‑AGCGCCAGCAGAACAATACCGATAATGCCAACCAGAAt                                                                                                   <  2:1329527/149‑1 (MQ=255)                                              ttCCAGTTCTGTTGCGCGTCTTTAATAAAGGAGGGCGTCACCATAATGTTGACGAGGAACTGCGCGCATCCTCCGTCGATAATGGC‑GG‑TTTGAATCGGTTACAGGCCGCCG‑AGCGCCAGCAGAACAATACCGATAATTCCAACCAGaa                                                                                                    <  1:1136429/148‑1 (MQ=255)                                                ccAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑GTTGAATCGGGTTCAGGCCGCCC‑AGCGCCCGCAGAAAAAAAACGCTAATGCCAACCAGaaattac                                                                                               >  2:492218/1‑146 (MQ=255)                                                 cAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAATCGGTTTCAGGCCGCCG‑AGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGAcc                                                                                              <  2:441005/151‑1 (MQ=255)                                                 cAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAATCGGGTTCCGGCCGCCC‑AGCGCCCGCCGAACAAAACCGCAAATGCCAACCAGaa                                                                                                    >  2:719854/1‑145 (MQ=255)                                                    ttCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTTAAACGGGTTTAGGCCCCCC‑CGCGCCAGCAGAAAAAAAACGATAAAGCCAACCAGAAATTACCaa                                                                                            >  2:1133806/1‑150 (MQ=255)                                                        gTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAATCGGTTTCAGGCCGCCG‑AGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGa                                                                                        <  2:906177/150‑1 (MQ=255)                                                            cgcgTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAATAGGGTTTCGGGCGCCG‑CGCGGCCGCAGAAAAAAAACGATAATGACAACCAGAATTTAACaaaccacac                                                                                     >  1:318265/1‑142 (MQ=255)                                                              cgTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTTAATCCGGTTTTGGGCGCCC‑CGCCGCCGCAGCACAAAACCCACAATTCCCACCAGAAATGAACAACCGATa                                                                                      >  1:799086/1‑146 (MQ=255)                                                               gTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTTAATCGGGTTTTGGCCGCCG‑AGCGCCCGCAGAACAAAACCGATAATGCCAAACAGAAATTACCAACCGATACACAcc                                                                                >  2:1093743/1‑151 (MQ=255)                                                                  tttAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAATCGGGGTTAAGCCCCCG‑AGCGCCAGCAGAAAAAAACCGATAAATACAAACAGAAATTAACAAACGACACccaccaac                                                                             >  1:1055182/1‑149 (MQ=255)                                                                   ttAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAATCGGTTTCAGGCCGCCG‑AGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACcag                                                                              <  1:86854/149‑1 (MQ=255)                                                                   ttAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAATCGGTTTCAGGCCGCCG‑AGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACca                                                                               <  2:1302217/148‑1 (MQ=255)                                                                    taataaAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAATCGGGTTTAGGCCGCCC‑AGCGCCAGCAGAAAAATAACCATAATGCCAACCAGAAATGACCCACCGATACGCACcaac                                                                             >  2:1223038/1‑147 (MQ=255)                                                                     aataaAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTTAATTTGTTTTTAGCCGCCC‑CGCGCCAGCCGAAAAATAACCATAAAGCCAACCACAAATGACCAACAGCTAAGCAcc                                                                                >  1:459002/1‑145 (MQ=255)                                                                     aataaAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAATCGGTTTCAGGCCGCCG‑AGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACcagcag                                                                           <  1:619707/150‑1 (MQ=255)                                                                     aataaAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAATCGGTTTCAGGCCGCCG‑AGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCaga                                                                          <  2:430557/151‑1 (MQ=255)                                                                          aGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAATCGGTTTCAGGCCGCCG‑AGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAggtgg                                                                     <  2:799123/151‑1 (MQ=255)                                                                             agagCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTTAATAGGGTTTAGGCCGCCC‑GGAGCCCGCAGAACAAAAACCATAAATGCAACCAGAAATTACCCACCGATACGCCCCAGCAGAggggggt                                                                   >  2:1143718/1‑148 (MQ=255)                                                                              gagCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAATCGGTTTCAGGCCGCCG‑AGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGttct                                                                 <  1:680978/151‑1 (MQ=255)                                                                              gagCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAATCGGTTTCAGGCCGCCG‑AGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGttct                                                                 <  2:860595/151‑1 (MQ=255)                                                                              gagCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAATCGGTTTCAGGCCGCCG‑AGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGttct                                                                 <  2:1020746/151‑1 (MQ=255)                                                                              gagCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAATCGGTTTCAGGCCGCCG‑AGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGttct                                                                 <  2:46956/151‑1 (MQ=255)                                                                               agCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGCATCGGGTTTAGGCCGCCG‑CGCGCCCGCAGAACAAAAACGCTAAAGTCAACCAGCAAATACCAAACCATACGCACCCCCAGAGGGGGgtctt                                                                >  1:1076744/1‑151 (MQ=255)                                                                               agCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAATCGGTTTCAGGCCGCCG‑AGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGttcgg                                                                <  2:968480/151‑3 (MQ=255)                                                                                    cACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAATCGGTTTCAGGCCGCCG‑AGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCa                                                           <  1:173109/151‑1 (MQ=255)                                                                                    cACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAATCGGGTTCTGGCCGCCG‑CGCGCCAGCAGAACAAAACCGATAATGCCAACCAGAAATGACCAACCGATAAGCACCAACAGAGGGGGGTCTTcccc                                                            >  2:185025/1‑150 (MQ=255)                                                                                     aCCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTTAATCGGTTTTAGGCCGCCG‑AGCGCCCGCAGAACAAAAACGAGAAAGCCAACCAAAAATTAACAACCCATACGCACCCACAGAGGGGGGTTTTCACCa                                                           >  1:1317194/1‑150 (MQ=255)                                                                                      ccATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAATCGGTTTCAGGCCGCCG‑AGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCAt                                                          <  2:276884/150‑1 (MQ=255)                                                                                      ccATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAATCGGTTTCAGGCCGCCG‑AGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATc                                                         <  2:818792/151‑1 (MQ=255)                                                                                       cATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAATCGGTTTCAGGCCGCCG‑AGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATcg                                                        <  2:331812/151‑1 (MQ=255)                                                                                          aaTGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTTAATCGGTTTTAAGCCGCCG‑AGCGCCCGCAGAAAAAAAACGATAAAGCCAACCAGAAATGACCAACCGATAAGCACCAGCAGAGGGGGTTTTTTACCCTcgc                                                       >  2:711731/1‑149 (MQ=255)                                                                                            tGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑CTTTAATGGGGGTCACGCCGCCC‑CGCGGCAGCAGAAAAAAAACCACAATGCCCAACACAACTGACCAACCCAAACGCCCCCGCAGAAGTGGTTTTTCTcccccacggac                                                   >  2:592684/1‑142 (MQ=255)                                                                                               tGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTTAATCGGGTTTAGGCCGCCG‑AGCGCCCCCAGAAAAAAACCCATAATAGCCACCAAAAATGACCAACCGAGACGCACCAGCAAAGGGGGGTCTTCTCCATCGCGGAc                                                   >  1:1278553/1‑148 (MQ=255)                                                                                                  cGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAATCGGTTTCAGGCCGCCG‑AGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTcgcg                                             <  2:840802/151‑1 (MQ=255)                                                                                                          aGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAATTGGGTTGAGGGCGCCC‑GGCCCCCGCCGAACAAAACCGAAAATGGCCACCAGAATAGAACAACCCAAACGCACCAGCAGAGGGGGGTTTTTATCACCCCGCAATTCGCGGCAAGAgg                                     >  1:1189645/1‑151 (MQ=255)                                                                                                             gggCATCCTCCGGCGATAATGGC‑GG‑TTTTGATCGGCTTTAAGCCCCCC‑GGCGCCCGCCGAACAAAAACCAAAAAACCAAACACAATTGAACAAACCATAACCACCAACAAAAGGGGGTCTTTTCCAAAGCCTACTTctccgccagggg                                     >  2:730674/1‑137 (MQ=255)                                                                                                              ggCATCCTCCGGCGATAATGGC‑GG‑TTTGAATCGGTTTCAGGCCGCCG‑AGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACAt                                 <  2:872534/151‑1 (MQ=255)                                                                                                                 aTCCTCCGGCGATAATGGC‑GG‑TTTTTATTGGTTTTCGGGCGCCC‑CGCGCCCCCCAAAAAAAACCCATAAAGACAAACACAAATGACCCACCCATACACACCCCCAAAGGTGGGTCTTTACCCTCCCGCAATTCCCGGCGAGGGGACa                                  >  1:1184321/1‑147 (MQ=255)                                                                                                                  tcctccGGCGATAATGGC‑GG‑TTTTAATCGGTTTCAGGCCGCCG‑AGCGCCAGCAGAACAATAACGATAATGCCAACCAGAAATTACCAACCGATACGCACCAGCAGAGGGGGGTCTTTACCATCGCGTACTTCGCGGCCAGTGGACATCGcc                             >  2:747376/1‑151 (MQ=255)                                                                                                                    ctccGGCGATAATGGC‑GG‑TTTTAATCGGGTTTAGGCCGCCG‑GGCGCCAGCAGAACAAAAACGATAAAGCCCAACAGAATTGACCAACCCATACGCACCAGCAGAGGGGGTGCTTTACCCACACGTACTTCGCGGGAAGGGGGCATCg                               >  1:747684/1‑147 (MQ=255)                                                                                                                       tggCGATAATTGC‑GG‑TTTGAATCGGTTTCAGGCCGCCG‑AGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGt                        >  1:956469/2‑151 (MQ=255)                                                                                                                              aaTGGC‑GG‑TTTGAATCGGTTTCAGGCCGCCG‑AGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAagag                 <  1:1022973/151‑1 (MQ=255)                                                                                                                              aaTGGC‑GG‑TTTGAATCGGTTTCAGGCCGCCG‑AGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAaga                  <  2:1352528/150‑1 (MQ=255)                                                                                                                               aTGTC‑GG‑TTTGAATAGGTTTCAGGCCGCCG‑AGCCCGAACAGAACAACACCGATAATGCCAGCCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCTCCATCGCGTTCTTCGCAGCAAGTGGACATCGCCAGGGTAGaaaa                  <  1:789046/149‑3 (MQ=255)                                                                                                                                 ggC‑GG‑TTTGAATCGGTTTCAGGCCGCCG‑AGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCa               <  1:605067/150‑1 (MQ=255)                                                                                                                                   c‑GG‑TTTGAATCGGTTTCAGGCCGCCG‑AGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGc            >  1:1280050/1‑151 (MQ=255)                                                                                                                                          tGAATCGGTTTCAGGCCGCCG‑AGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTa        <  2:382993/150‑1 (MQ=255)                                                                                                                                             aTCGGTTTCAGGCCGCCG‑AGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCa    <  1:632308/151‑1 (MQ=255)                                                                                                                                               cGGTTTCAGGCCGCCG‑AGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAg   <  1:471189/150‑1 (MQ=255)                                                                                                                                                ggTTTCAGGCCGCAG‑GGCGCCAGCAGACCAATTCCGATAATACCATCCAGAATTGACCATCCGATCCGCACCAGCAGAGGTGGTTCTTCGCCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGTGTATAAGAGCAGGCGTTAACCAGc  <  1:450763/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  AAATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAATCGGTTTCAGGCCGCCG‑AGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGC  >  NC_000913/40379‑40669

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