Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	40,515	0	G	T	37.8%	102.3 / 36.5	82	Q473K (CAA→AAA)	caiT	putative transporter
Reads supporting (aligned to +/- strand):  ref base G (16/35);  new base T (31/0);  total (47/35)
Fisher's exact test for biased strand distribution p-value = 4.33e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

CATCATTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAATCGGTTTCAGGCCGCCG‑AGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTA  >  NC_000913/40365‑40663                                                                                                                                                         |                                                                                                                                                      catcatTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTg                                                                                                                                                       <  2:433523/151‑1 (MQ=255)     atTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAATc                                                                                                                                                   <  2:506263/151‑1 (MQ=255)       ttaaattaaaATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTTAATTgg                                                                                                                                                 >  1:233850/1‑151 (MQ=255)             taaaATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTTAATAGGgttt                                                                                                                                             >  1:1287869/1‑148 (MQ=255)               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)                        tttCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTTGATCCGttttcagaccccgg                                                                                                                                   >  2:723545/1‑137 (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                                                           <  2:860595/151‑1 (MQ=255)                                                                                            gagCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAATCGGTTTCAGGCCGCCG‑AGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGttct                                                           <  1:680978/151‑1 (MQ=255)                                                                                            gagCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAATCGGTTTCAGGCCGCCG‑AGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGttct                                                           <  2:46956/151‑1 (MQ=255)                                                                                            gagCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAATCGGTTTCAGGCCGCCG‑AGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGttct                                                           <  2:1020746/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)                                                                                                                                                         |                                                                                                                                                      CATCATTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAATCGGTTTCAGGCCGCCG‑AGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTA  >  NC_000913/40365‑40663

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: ‑