Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	40,515	0	G	T	41.6%	136.1 / 47.1	89	Q473K (CAA→AAA)	caiT	putative transporter
Reads supporting (aligned to +/- strand):  ref base G (11/41);  new base T (37/0);  total (48/41)
Fisher's exact test for biased strand distribution p-value = 1.86e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGC‑AGAACAATACCGATAATGCCAACC‑AGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAAC  >  NC_000913/40368‑40665                                                                                                                                                    |                                                                                                                                                          catTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGaa                                                                                                                                                         <  1:1503468/150‑1 (MQ=255)  atTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAAT‑c                                                                                                                                                      <  1:963118/151‑1 (MQ=255)  atTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAAT‑c                                                                                                                                                      <  1:685456/151‑1 (MQ=255)     taaattaaaATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTACTAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTAAT‑Cggg                                                                                                                                                   >  2:1308880/1‑150 (MQ=255)      aaattaaaATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAACAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAAT‑CGGtt                                                                                                                                                  <  1:312986/151‑1 (MQ=255)        attaaaATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAAT‑CGGTTTc                                                                                                                                                <  2:1318125/151‑1 (MQ=255)        attaaaATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAAT‑CGGTTTc                                                                                                                                                <  1:421559/151‑1 (MQ=255)          taaaATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTAAT‑CGGgtt                                                                                                                                                 >  2:693922/1‑148 (MQ=255)              aTCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTAAT‑TGGTTTTAGgcc                                                                                                                                           >  2:1551191/1‑150 (MQ=255)              aTCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAAT‑CGGTTTCAGgccg                                                                                                                                          <  1:63929/151‑1 (MQ=255)                ccATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTAAT‑CGGGTTCTGgccgcc                                                                                                                                        >  2:1003852/1‑151 (MQ=255)                    cTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTAAT‑TGGgtt                                                                                                                                                 >  1:276959/1‑138 (MQ=255)                      ttCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTAAT‑TAGGTTCAGGCCCCCGacc                                                                                                                                    >  1:904062/1‑147 (MQ=255)                      ttCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAAT‑CGGTTTCAGGCCGCCGAgc                                                                                                                                    <  2:464069/149‑1 (MQ=255)                       tCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTAAT‑TGGTTTCAGGGCGCCCGGCGcc                                                                                                                                 >  1:1466323/1‑151 (MQ=255)                       tCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGATTGTAT‑CGGGTTCACGGCGCCGAGCGcc                                                                                                                                 >  1:901011/1‑151 (MQ=255)                         aaCCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTTAT‑CGGGTTCAGGGCGCCGAgcgc                                                                                                                                  >  1:954363/1‑148 (MQ=255)                         aaCCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTAAT‑TGGGTTTAgggcgccgggcgcc                                                                                                                                 >  2:859549/1‑149 (MQ=255)                           ccTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAAT‑CGGTTTCAGGCCGCCGAGCGCcagc‑a                                                                                                                            <  2:58092/151‑1 (MQ=255)                           ccTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAAT‑CGGTTTCAGGCCGCCGAGCGCcagc‑a                                                                                                                            <  2:669544/151‑1 (MQ=255)                           ccTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAAT‑CGGTTTCAGGCCGCCGAGCGCcagc‑a                                                                                                                            <  1:778976/151‑1 (MQ=255)                           ccTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAAT‑CGGATTCAGGCCGCCGAGCGCcagc‑a                                                                                                                            <  2:153824/151‑1 (MQ=255)                              ctTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAAT‑CGGTTTCAGGCCGCCGAGCGCcagc‑ag                                                                                                                           <  1:1111906/149‑1 (MQ=255)                               ttGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGGTTTAAT‑TGGGTTTATGGCGCCGCGCGCccgc                                                                                                                              >  1:1421667/1‑146 (MQ=255)                                tGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGC‑AGAAc                                                                                                                        <  2:1667221/150‑1 (MQ=255)                                 gATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTGAT‑CGGGTTTAGGCCGCCGGGCGCCAGC‑AGaaa                                                                                                                        >  2:1680061/1‑148 (MQ=255)                                 gATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGC‑AGAACaa                                                                                                                      <  1:318029/151‑1 (MQ=255)                                  atatTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGC‑AGAACAAt                                                                                                                     <  1:1413720/151‑1 (MQ=255)                                    atTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGC‑AGAACAATAc                                                                                                                   <  2:1018497/151‑1 (MQ=255)                                    atTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAAT‑CGGGTTCAGGCCGCCGAGCGCCAGC‑AGAACaaaa                                                                                                                    >  2:420069/1‑148 (MQ=255)                                      tttGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTTAT‑CGGTTTTTAGCCGCCGCGCGCCCGC‑AGAaaaaaaac                                                                                                                  >  1:689799/1‑146 (MQ=255)                                       ttGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGCTTTAAT‑CGGTTTCAGGCCCCCGACGGGCAGC‑AGAACAATaac                                                                                                                  >  2:1276122/1‑147 (MQ=255)                                                aattaatCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTAAT‑CGGGTTTAGGCCGCCGGGCGCCAGC‑AGAACAAAACCCATAATGCCa                                                                                                        >  2:40628/1‑150 (MQ=255)                                                aattaatCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTAAA‑CGGGTTTACGCCGCCCAGCGCCCGC‑AGAACAATAACGATaaaggc                                                                                                         >  2:1277765/1‑145 (MQ=255)                                                  ttaatCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTAAA‑TGGGTTCAGGGCGCCGAGCGGCCGG‑AGAAAAAAAACGATAAAGCCaa                                                                                                       >  2:65335/1‑149 (MQ=255)                                                    aatCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGC‑AGAACAATACCGATAATGCCAACC‑Ag                                                                                                  <  1:572856/151‑1 (MQ=255)                                                     atCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTAAT‑CGGGTTTTGGGCGCCGCGCGCCCGC‑AGAAAAATATCAATAATgacca                                                                                                       >  1:1442887/1‑142 (MQ=255)                                                      tCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTAAT‑CGGGTTTAGGCCCCCCAGCGCCCGC‑AGAAAAAAACCGATAATGTCaaac                                                                                                     >  2:20840/1‑145 (MQ=255)                                                        tttCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTAAA‑AGGGTTTAGGGCGCCGGGCGCCCGC‑GAAACAAAACCCATAAAGGCCACC‑AGaa                                                                                                >  2:531279/1‑149 (MQ=255)                                                         ttCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGTAT‑AGGTTTTAGGGCGCCGAGGGCCACG‑AGAACAAAAACCATAAAGCCAAAC‑AGaa                                                                                                >  1:1206690/1‑148 (MQ=255)                                                          tCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTAAT‑CGGTTTTATGCCGCCGGGCGCCCGC‑AGAACAATAACGATAATTcccaaccga                                                                                                  >  1:1005274/1‑139 (MQ=255)                                                           ccAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTGAT‑AGGGTTTTAGCCGCCCCGCGCCCGC‑AGAAAAAAAACGATAATACCAACC‑aaaa                                                                                                >  1:826306/1‑146 (MQ=255)                                                           ccAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGC‑AGAACAATACCGATAATGCCAACC‑AGAATTGAc                                                                                           <  2:748815/151‑1 (MQ=255)                                                            cAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGC‑AGAACAATACCGATAATGCCAACC‑AGAATTGAcc                                                                                          <  2:383670/151‑1 (MQ=255)                                                              gTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTTAT‑CGGGTTCAGGCCGCCCAGCGCCAGC‑AGAAAAAAAACCATAAAGCCAAAC‑AAAAATGaac                                                                                          >  2:1055459/1‑147 (MQ=255)                                                                tCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGC‑AGAACAATACCGATAATGCCAACC‑AGAATTGACCAAc                                                                                       <  2:837297/150‑1 (MQ=255)                                                                   gTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGC‑AGAACAATACCGATAATGCCAACC‑AGAATTGACCAAcc                                                                                      <  1:1558435/148‑1 (MQ=255)                                                                    tttCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGGTTGAAA‑TGGGGTTAGGCCGGCGAACCCCCGG‑GGAGAAAAAACGATAATGcccaaacgaagatac                                                                                            >  2:278900/1‑129 (MQ=255)                                                                        gcgTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTAAT‑CGGGTTTAGGCCGCCGGGCGCCAGC‑AGAAAAAAAACCATAATGACAACC‑AGAAATGACCCACCCATAccc                                                                               >  2:54065/1‑148 (MQ=255)                                                                         cgTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGC‑AGAACAATACCGATAATGCCAACC‑AGAATTGACCAACCGATACGCAc                                                                             <  1:717332/151‑1 (MQ=255)                                                                              ttAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTTAT‑AGGGTTTAAGCCGCCGAGCGCCCGC‑AAAACAATAACGAAAATTCCCACC‑AGAATTTACCCACCGATACGCaacccc                                                                         >  1:1548471/1‑145 (MQ=255)                                                                              ttAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTAAA‑CGGGTTCAGGCCGCCCAGGGCCCGC‑AGAACAATACCGATAATGCCAACC‑AGAATTGACCCACCGATACGCACccgc                                                                         >  2:1500460/1‑150 (MQ=255)                                                                              ttAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGGAT‑AGGGTTTCGGGCGCCGGGGGCCAGG‑AGAACAAAACCGATAATTACAAAC‑CGAATTGAACAAACCATAAGCAc                                                                             >  1:746351/1‑146 (MQ=255)                                                                              ttAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGC‑AGAACAATACCGATAATGCCAACC‑AGAATTGACCAACCGATACGCACca                                                                           <  2:1095709/148‑1 (MQ=255)                                                                              ttAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGGTTTAATGGGGTTTACGGCCG‑CGCGCGCCCGC‑AGAAAAAAACCGAAAATACCCACC‑AGAAATGAACAAACCATACGCAcccac                                                                         >  1:1007436/1‑147 (MQ=255)                                                                                aataaAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGC‑AGAACAATACCGATAATGCCAACC‑AGAATTGACCAACCGATACGCACcagcag                                                                       <  1:709446/150‑1 (MQ=255)                                                                                aataaAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGC‑AGAACAATACCGATAATGCCAACC‑AGAATTGACCAACCGATACGCACCAGCaga                                                                      <  2:1560114/151‑1 (MQ=255)                                                                                   aaaGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTATT‑ACAGTTCTCGCCGCCGCGCGCCACC‑AGGACAAAAATCATAATGTCCACC‑CGAAATTACCAAACGAAACGACCCCCCAGAggg                                                                   >  1:497713/1‑150 (MQ=255)                                                                                   aaaGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGC‑AGAACAATACCGATAATGCCAACC‑AGAATTGACCAACCGATACGCACCAGCAGAgg                                                                    <  1:1500914/150‑1 (MQ=255)                                                                                    aaGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTAAT‑TGGTTTTAGGCCGCCGCGCGCCCAC‑CGAACAAAACCGATAAAGCCACCC‑CGAAATGAACAACCCACACGCCCccgcaaaaggg                                                                  >  1:1079728/1‑145 (MQ=255)                                                                                    aaGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGC‑AGAACAATACCGATAATGCCAACC‑AGAATTGACCAACCGATACGCACCAGCAGAggtg                                                                  <  2:1465027/151‑1 (MQ=255)                                                                                     aGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGC‑AGAACAATACCGATAATGCCAACC‑AGAATTGACCAACCGATACGCACCAGCAGAggtgg                                                                 <  1:568738/151‑1 (MQ=255)                                                                                     aGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAAT‑CGGGTTTAGGCCGCCCAGCGCCCGC‑AGAACAAAAACGATAATGCCAACC‑AGAAATGACCAACCCATACGCACCAGCAGAggggg                                                                 >  1:1694747/1‑151 (MQ=255)                                                                                        agagCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTAAT‑CGGGTTTAGGCCGCCCAGCGCCCGC‑AGAAAAAAAACGATAATGCCAAAC‑AGAAATGACCCACCCATACGCACCCGCAGAggggggt                                                               >  2:139430/1‑148 (MQ=255)                                                                                         gagCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGC‑AGAACAATACCGATAATGCCAACC‑AGAATTGACCAACCGATACGCACCAGCAGAGGTGGttct                                                             <  1:513541/151‑1 (MQ=255)                                                                                            cGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTTAT‑CGGGTTTTGGCCGCCGAGGGCCCGC‑CGAAAAAAAACCATAATGCCCACC‑AAAAATGACCCACCCAAACGCCCCACCAAAGGGGGgtctt                                                            >  1:154081/1‑149 (MQ=255)                                                                                            cGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGC‑AGAACAATACCGATAATGCCAACC‑AGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCa                                                          <  2:1150987/151‑1 (MQ=255)                                                                                               cACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTTAT‑CGGGTTTAGGGCGCCGAGCGCCCGC‑AGAACAAAAACCAAAAAGCCCAAC‑AGAAATGAACAAACCCTACGCACCCCCAGAGGGGGGTCTTTAcc                                                        >  1:1401747/1‑150 (MQ=255)                                                                                               cACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGC‑AGAACAATACCGATAATGCCAACC‑AGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCa                                                       <  2:1277839/151‑1 (MQ=255)                                                                                                aCCATAATGTTGATGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGC‑AGAACAATACCGATAATGCCAACC‑AGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCa                                                       <  1:1592353/150‑1 (MQ=255)                                                                                                aCCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTAAT‑CGGGTTTAAGTCGCCGAGCCCCCAC‑AGAACAAAACCCATAATGCACAAC‑AGAAATGATCCACCGATCAGCCCCACCAGCggggggt                                                               >  1:718705/1‑140 (MQ=255)                                                                                                 ccATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGC‑AGAACAATACCGATAATGCCAACC‑AGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATc                                                     <  2:237663/151‑1 (MQ=255)                                                                                                  cATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTAAT‑CGGGTTCCGGCCGCCCCCGGCCAGC‑ACGACAAAACCCATAAAGGCAACCAAAAATTAACAAACCGAAAACCACCA‑CAGAggggggatttta                                                          >  1:1406884/1‑138 (MQ=255)                                                                                                  cATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGC‑AGAACAATACCGATAATGCCAACC‑AGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATcg                                                    <  2:381652/151‑1 (MQ=255)                                                                                                   aTAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTAAA‑CGGGTTTAGGCCGCCGAGCGCCCGC‑AGAAAAATACCCATAAAGCCAACC‑AGAATTTACCAAACGCGAAGCACCCGCAGAGGGGGGTCTTTAcc                                                        >  2:744482/1‑146 (MQ=255)                                                                                                    tAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTAAA‑CGGTTTTCGGCCGCCCGGCGCCCGC‑AGAAAAAAACCGAAAATTCCCACC‑CAAATTTACCAACCGAAACGCACCAGCAGAGGGGGGTCTTTAcccccgc                                                   >  2:1026370/1‑145 (MQ=255)                                                                                                         ttGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTAAA‑CGGGTTCAGGCCGCCCGGCGCCAGCAAGAAAAAAACC‑ATAAATCCAACC‑AGAAATGACCAACCGAAACGCACCAGCAGAGGGGGGTCTTTACCAACGCGCACtt                                             >  2:1692048/1‑151 (MQ=255)                                                                                                           gACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGTAA‑TGGGTTCAGGGCGCCCGGCGCCCAC‑CGAAAAAAACCCATAAAGGCAACC‑AGAAATGAACAAACCAAAAGCCCCCGCAGAGGGGGGTTTTCTCCACCGCCTAcct                                             >  1:699846/1‑147 (MQ=255)                                                                                                                  aaCAGCGGGCATCCTCCGGCGATAATGGCGGTTTTAAT‑TGGTTTTCGGGCGGCGGGCGGCCGC‑AGAACAAAAACGAAAAATCCAAAC‑AGAATTGGCCAACCCATACGCACCCGAAGAGGGGGTTCTTCTCCATcgccacct                                              >  2:324330/1‑136 (MQ=255)                                                                                                                          gCATCCTCCGGCGATAATGGCGGTTTGAAT‑TGGTTTTAGGCCGCCGCGCGCCCGC‑CGAAAAATAACGATAATACCAACC‑CAAATTTACCAACCCATAACCACCACCAGAAGGGGGTTTTTACCATCCCGTACCTCGCGGCCAGTgg                                 >  1:733287/1‑146 (MQ=255)                                                                                                                          gCATCCTCCGGCGATAATGGCGGTTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGC‑AGAACAATACCGATAATGCCAACC‑AGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATc                            <  2:838770/151‑1 (MQ=255)                                                                                                                                 ccGGCGATAATGGCGGTTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGC‑AGAACAATACCGATAATGCCAACC‑AGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAggg                     <  1:1230544/151‑1 (MQ=255)                                                                                                                                 ccGGCGATAATGGCGGTTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGC‑AGAACAATACCGATAATGCCAACC‑AGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAgg                      <  2:824489/150‑1 (MQ=255)                                                                                                                                           tGGCGGTTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGC‑AGCACAATACCGATAATGCCAACC‑AGAATTGACCAACCGATACGCACCAGCAGAGGGGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGc            >  1:1074216/1‑150 (MQ=255)                                                                                                                                            ggCGGTTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGC‑AGAACAATACCGATAATGCCAACC‑AGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCa           <  2:1555487/150‑1 (MQ=255)                                                                                                                                              cGGTTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGC‑AGAACAATACCGATAATGCCAACC‑AGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGc        <  1:5202/151‑1 (MQ=255)                                                                                                                                               ggTTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGC‑AGAACAATACCGATAATGCCAACC‑AGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCg       >  1:287761/1‑151 (MQ=255)                                                                                                                                               ggTTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGC‑AGAACAATACCGATAATGCCAACC‑AGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCg       <  2:1696280/151‑1 (MQ=255)                                                                                                                                                   tGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGC‑AGAACAATACCGATAATGCCAACC‑AGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTa    >  2:48016/1‑150 (MQ=255)                                                                                                                                                    gAAT‑CGGTTTCAGGCCGCCGAGCGCCAGC‑AGAACAATACCGATAATGCCAACC‑AGAATTGACCAACCGATACGCACCAGCAGGGGTGGTTCTTCACCATCGCGTACTTCGCGGCATGGGGACATCGCCAGGGGATAAGAGCAGGCGTTAAc  >  2:1436485/1‑151 (MQ=255)                                                                                                                                                    gAAT‑CGGTTTCAGGCCGCCGAGCGCCAGC‑AGAACAATACCGATAATGCCAACC‑AGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTaa   <  2:1071176/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                          CATTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGC‑AGAACAATACCGATAATGCCAACC‑AGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAAC  >  NC_000913/40368‑40665

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