Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	40,515	0	G	T	44.3%	93.2 / 33.6	80	Q473K (CAA→AAA)	caiT	putative transporter
Reads supporting (aligned to +/- strand):  ref base G (13/31);  new base T (35/0);  total (49/31)
Fisher's exact test for biased strand distribution p-value = 6.10e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCATTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAAC  >  NC_000913/40366‑40665                                                                                                                                                       |                                                                                                                                                        atcatTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTGa                                                                                                                                                        <  1:1303911/151‑1 (MQ=255)    atTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTGCAt                                                                                                                                                      >  2:1134586/1‑150 (MQ=255)      ttaaattaaaATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTTAATAgg                                                                                                                                                   >  1:1069042/1‑151 (MQ=255)       taaattaaaATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTGAATCggg                                                                                                                                                  >  2:733822/1‑150 (MQ=255)        aaattaaaATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTTAATTgggt                                                                                                                                                 >  2:777570/1‑149 (MQ=255)             aaaaTCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTGAATCGGTTTCAgg                                                                                                                                            >  1:923781/1‑151 (MQ=255)              aaaTCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTTGATCGGGTTCAGgc                                                                                                                                           >  2:281461/1‑151 (MQ=255)              aaaTCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTTAAACGGGTTctagg                                                                                                                                           >  2:879510/1‑147 (MQ=255)               aaTCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTGAATCGGTTTCAGgcc                                                                                                                                          <  2:112743/151‑1 (MQ=255)                  ccATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTGAATCGGTTTCAGgccgcc                                                                                                                                       <  2:529693/151‑1 (MQ=255)                   cATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTTTATAGGTTTTAAgccgccc                                                                                                                                      >  2:579744/1‑150 (MQ=255)                                ctTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTTAATTGGGTTTAGGCCGCCGAGCGCCCGCCGa                                                                                                                          >  1:1122769/1‑150 (MQ=255)                                  tGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGaa                                                                                                                         <  1:1350622/149‑1 (MQ=255)                                      atTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATa                                                                                                                   <  1:239002/150‑1 (MQ=255)                                      atTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATa                                                                                                                   <  1:471382/150‑1 (MQ=255)                                      atTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATAc                                                                                                                  <  2:473059/151‑1 (MQ=255)                                       ttttGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTGTATTGGTTTTAGGACGCCGCGCGCCCGCAGAaa‑aaaaac                                                                                                                 >  2:1247563/1‑147 (MQ=255)                                         ttGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTGTATCGGGTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATAACg                                                                                                                >  1:232746/1‑150 (MQ=255)                                           gggggTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTTGATCGGGTTCAGGCCGCCCCGCGGcagccaaaacaaaacgca                                                                                                               >  2:1226965/1‑136 (MQ=255)                                             gggTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTTAATCGGTTTTAAGCCCCCCCGCGCCCGCCGAAC‑AAAAACCATAAt                                                                                                           >  1:61409/1‑151 (MQ=255)                                              ggTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTTAATCGGGTTCAGGCCGCCGCGCGCCAGCAGAAC‑AAAACCGaaaaat                                                                                                          >  1:420872/1‑149 (MQ=255)                                               gTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTTAATCGGGTTCAGGCCCCCGAGCGCCCGCAGAAC‑AAAACCGATaaagg                                                                                                         >  1:58036/1‑148 (MQ=255)                                                ttaattaaTCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTTAATAGGGTTTAGGGCGCCGAGCGCCAGCCGAAA‑AAAAACGATaaattc                                                                                                        >  2:892415/1‑147 (MQ=255)                                                 taattaatCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTTAATAGGGTTTAGGCCGCCCCGCGCCCGCCGAAA‑AAAAcccaaaaaagcc                                                                                                       >  2:833295/1‑141 (MQ=255)                                                  aattaatCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTTTATTGGGTTTAGGCCGCCGAGCGCCCGCAGAAA‑AATAACCATAAttgca                                                                                                       >  2:398961/1‑146 (MQ=255)                                                   attaatCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTTAATCGGGTTTAGGCCCCCGAGCGCCCGCAGAAC‑AATAACGATAATTCCAAc                                                                                                     >  1:188779/1‑151 (MQ=255)                                                    ttaatCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTTAATCGGGTTTAGGCCGCCCAGCGCCCGCAGAAC‑AATAACGAAAATGCCAAcc                                                                                                    >  1:720888/1‑151 (MQ=255)                                                      aatCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTTAATAGGTATAAGGCCGCCGCGCGCCCGCAGAAC‑AAAACCGATAATGCCAACCa                                                                                                   >  2:1654/1‑150 (MQ=255)                                                       atCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTGAAATGGTTTTAAGCCGCCGAGCGCCAGCAGAAA‑AATAACGATAAAACCAAca                                                                                                    >  2:948015/1‑147 (MQ=255)                                                         cTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTGAATCGGTTTTAGGCCGCCCCGCGCCCCCAGAAC‑AAAACCGATAAAACCAACCaaaa                                                                                                >  1:1301130/1‑150 (MQ=255)                                                            tCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTTGATCGGGTTTAGGGCGCCCCGCGCCCGCAGAAC‑AATAACGGTAAAGCCaaaccaaaatg                                                                                             >  1:1321530/1‑141 (MQ=255)                                                             ccAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAATGCCAACCAGAATTGAc                                                                                           <  1:1013461/151‑1 (MQ=255)                                                              cAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTTTATTGGGTTTAGGGCGCCCCGCGCCACCCGAACAAAAACCAAAAATGCAAACCAAAAttacc                                                                                           >  2:1322597/1‑148 (MQ=255)                                                              cAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAATGCCAACCAGAATTGAcc                                                                                          <  1:38379/151‑1 (MQ=255)                                                               aGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAATGCCAACCAGAATTGACCa                                                                                         <  1:1189934/151‑1 (MQ=255)                                                                gTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAATGCCAACCAGAATTGACCaa                                                                                        <  1:928144/151‑1 (MQ=255)                                                                    tGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTTAATCGGTTTTCGGGCGCCCAGCGCCCGGAGAAA‑AAAACCCATAAAACCAACCAGAAATGAcccaacgc                                                                                    >  2:266945/1‑145 (MQ=255)                                                                     gTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTAAATGGGTTTAAGG‑CGCCGAGGGCCCGCAGGAC‑AAAAACCAAAAAGCCAAACAAAAATGAACAAcccac                                                                                   >  2:994465/1‑148 (MQ=255)                                                                     gTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAATGCCAACCAGAATTGACCAAcc                                                                                      <  1:31515/148‑1 (MQ=255)                                                                     gTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAATGCCAACCAGAATTGACCAACCGa                                                                                    <  1:501399/150‑1 (MQ=255)                                                                       ttCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAATGCCAACCAGAATTGACCAACCGATa                                                                                  <  1:979582/150‑1 (MQ=255)                                                                          gcgTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAATGCCAACCAGAATTGACCAACCGATACGCa                                                                              <  1:1142866/151‑1 (MQ=255)                                                                             tCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑GTTTAATCGGGTTTAGGCCCCCGAGCGCCCGCCGAAC‑AAAACCGAAAATGCCAACCAGAAATGAACAACCCAAACGCAccc                                                                           >  2:283328/1‑150 (MQ=255)                                                                              cTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTTATTGGGTTTTAGG‑CGCCCCGGGGCCGCAGAAA‑AAAAACCAAAAAGCCAACCAAAAATGACCAACCGATAAGCAccc                                                                           >  2:1081689/1‑149 (MQ=255)                                                                                ttAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTTAATCGGTTTTAGGCCGCCGCGCGCCCGCAGAAA‑AAAAACGATAAAGCCAACCAGAAATTACCAACCGATACGCACcagc                                                                         >  1:253876/1‑150 (MQ=255)                                                                                ttAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACcagc                                                                         <  2:1255701/150‑1 (MQ=255)                                                                                ttAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACca                                                                           <  1:486903/148‑1 (MQ=255)                                                                                 taataaAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACcagca                                                                        <  2:235350/150‑1 (MQ=255)                                                                                  aataaAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACcagcag                                                                       <  2:644276/150‑1 (MQ=255)                                                                                      aaGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTTGATCGGTTTCAGGCCGCCGAGCGCCCGCAGAAA‑AAAAACGATAATGCCCACCAGAATTGACCAACCGATACGCACCAGCAGAgggg                                                                  >  1:457896/1‑149 (MQ=255)                                                                                       aGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTTAATTGG‑TTCAAGCCCCCCAGCGCCCGCAGAAA‑AAAAACCATAAAACCAAACCGAAATGACCAACCCAAAAGCAcccaccgaggggg                                                                 >  2:523677/1‑140 (MQ=255)                                                                                       aGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTGAATCGGTTTTAGGCCGCCGGGCGCCCGCAGAAC‑AAAACCGATAATGCCAACCAGAAATGACCAACCGATACGCACCACCAGAggggg                                                                 >  2:856836/1‑151 (MQ=255)                                                                                          agagCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTTTATAGGGTTTTTGGCGCCGAGCGCCCGCAGAAA‑AAAAACGCTAAAGCCAAACACAATTGAACAACCGATAAACACCCACAGAggggggtt                                                              >  2:451236/1‑148 (MQ=255)                                                                                          agagCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTTAAACGGTTTTAGGGCCCCGAGCGCCCGCAGGAA‑AAAAACCATAATGGCCACCAGAATTGACCAACCGAAACGCACCAGCAGAggggggtt                                                              >  1:932146/1‑148 (MQ=255)                                                                                            agCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTTTATCGGTTTTAGGGCGCCCAGCGCCCGCCGAAC‑AAAAACGAAAATGCCCACCACAAAATACCAACCCACACCCACCAACAGAGGGGGgtctt                                                            >  2:182660/1‑151 (MQ=255)                                                                                            agCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGttctt                                                            <  2:999895/151‑1 (MQ=255)                                                                                               gTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTGAATCGGGTTTAGGCCCCCGCGCGCCCGCAGAAA‑AAAACCGATAATGCCAACCCGAAATGACCAACCGATAAGCACCAGCAGAGGGGGGTCTTCAc                                                         >  2:760784/1‑151 (MQ=255)                                                                                                 cACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTTAAAAGGGTTTCGGCCGCCGCGCGGCAGCAGAAC‑AAAAACCAAAAAGACAACCCGAAATTACCCACCCACACGCACCCGCAGGGGGGGGTCTTCAccc                                                       >  2:52622/1‑150 (MQ=255)                                                                                                  aCCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTTAATCGGGTTTATGGCGCCCGGCGCCCGCCGAAC‑AAAAACCATAATGCCAAACAGAAATGACCAAACCATACGCACCAGCAGAGGGGGGTCTTTACCa                                                       >  2:1044925/1‑150 (MQ=255)                                                                                                    cATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATcg                                                    <  2:1243110/151‑1 (MQ=255)                                                                                                    cATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATcg                                                    <  1:700744/151‑1 (MQ=255)                                                                                                     aTAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTTAATCGGGTTTAGGCCGCCCAGGGCCAGCAGAAC‑AAAAACCATAATGCCCACCAAAATTGACCAACCCATACGCACCAGAAAAGGGGGGTCTTTACCATccc                                                   >  1:1010186/1‑149 (MQ=255)                                                                                                        aTGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTTTATTGGGTTTTGGGCGCCCCGCGCCCGCAGAAC‑AAAACCGATAAAGCCAAACAGAAATTACCAACCGATACGCCCCAGCAGAGGGGGGTCTTTTCCATcgcgg                                                 >  2:781362/1‑149 (MQ=255)                                                                                                         tGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTGAATCGGTTTTAGGCCGCCCAGCGCCCGCAGAAC‑AAAACCGATAAAGCCCACCAAAATTGAACAACCGATACGCACCAGCCGAGGGGGGTCTTCACCCTACCGTAc                                               >  2:1200207/1‑151 (MQ=255)                                                                                                          gTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTGAATAGGTTTTAGGCCGCCCAGCGCCCGCAGAAC‑AAAAACGATAATGTCAAACAGAAATGACCAAACGATACGCACCAGCAGAGGGGGGTGTTTCACATcgcggaat                                              >  2:1200163/1‑147 (MQ=255)                                                                                                             gACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTTTATCGGGTTCAGGCCCCCCGGCGCCCGCCGAAC‑AATACCGCTAATGCCAACCACAAATGACCCACCGCTACGCCCCCGCAGAGGTGGTTTTTTACCATAGCCTAccttg                                           >  1:866644/1‑147 (MQ=255)                                                                                                             gACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTTAATCGGTTTTAGGCCGCCGCGCGCCAGCAGAAC‑AAAAACGATAATGCCAACCAGAAATGAACAACCCATACGCACCAACAGAGGGGGGTTTTTACCCACGCGTACTTcc                                           >  1:760204/1‑150 (MQ=255)                                                                                                                   gaaCAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCaa                                     <  1:728049/151‑1 (MQ=255)                                                                                                                    aaCAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTTTATCGGGTTTAGGCCGCCGAGCGCCAGCCGAAC‑AAAAACGATAATGCCCACCAGAAATTACCAACCCATAAGCCCCCGCAGAGGGGGGTCTTTACCATCCCGTACTTCGCGGCaa                                     >  1:1245285/1‑150 (MQ=255)                                                                                                                          gggCATCCTCCGGCGATAATGGCGG‑TTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACa                              <  1:164556/151‑1 (MQ=255)                                                                                                                             cATCCTCCGGCGATAATGGCGG‑TTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTACTTCACCATCGCGTACTTCGCGGCAAGTGGACATCg                           <  1:1021614/151‑1 (MQ=255)                                                                                                                              aTCCTCCGGCGATAATGGCGG‑TTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGc                          <  1:1213568/151‑1 (MQ=255)                                                                                                                              aTCCTCCGGCGATAATGGCGG‑GTTTAATCGGTTTCAGGCCGCCCCGCGCCCGCAGAAC‑AAAAACGAAAAAGCCAACCAGAATTTACCCACCGATACGCACCAGCAGAGGGGGGTCTTTACCATCGCCGACTTCGCGGCAAGGGGGCAACGc                          >  2:216082/1‑151 (MQ=255)                                                                                                                                cctccGGCGATAATGGCGG‑TTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCa                        <  1:1335006/151‑1 (MQ=255)                                                                                                                                   ccGGCGATAATGGCGG‑TTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAgg                      <  1:26558/150‑1 (MQ=255)                                                                                                                                   ccGGCGATAATGGCGG‑TTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAgg                      <  1:812219/150‑1 (MQ=255)                                                                                                                                     ggCGATAATGGCGG‑TTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGta                   <  2:685177/151‑1 (MQ=255)                                                                                                                                             tGGCGG‑TTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGc            >  1:1150205/1‑150 (MQ=255)                                                                                                                                                     ttGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGtt     >  1:399171/1‑150 (MQ=255)                                                                                                                                                     ttGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGtt     <  2:814307/150‑1 (MQ=255)                                                                                                                                                       gAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTaa   <  2:372909/150‑1 (MQ=255)                                                                                                                                                       gAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAAc  >  1:614415/1‑151 (MQ=255)                                                                                                                                                       gAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAAc  <  1:869296/151‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                        ATCATTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAAC  >  NC_000913/40366‑40665

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