Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	40,515	0	G	T	42.4%	94.1 / 32.9	66	Q473K (CAA→AAA)	caiT	putative transporter
Reads supporting (aligned to +/- strand):  ref base G (11/27);  new base T (28/0);  total (39/27)
Fisher's exact test for biased strand distribution p-value = 5.62e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑GAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAA‑TGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTA  >  NC_000913/40372‑40663                                                                                                                                                 |                                                                                                                                                        taaattaaaATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑TAAT‑Agg                                                                                                                                                  >  2:54806/1‑150 (MQ=255)  aaattaaaATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑GAAT‑CGGtt                                                                                                                                                <  1:281599/151‑1 (MQ=255)   aattaaaATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑GAAT‑CGGttt                                                                                                                                               <  1:819571/151‑1 (MQ=255)   aattaaaATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑GAAT‑CGGttt                                                                                                                                               <  1:392808/151‑1 (MQ=255)           tCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCGTAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑GAAT‑CGGGTTTAGgcccc                                                                                                                                       >  2:842409/1‑149 (MQ=255)             cATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑TAAT‑CGGGTTTAGggcgcc                                                                                                                                      >  1:187245/1‑150 (MQ=255)             cATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑GAAT‑CGGGTTCAAgccgcc                                                                                                                                      >  2:484433/1‑150 (MQ=255)                cTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑GAAT‑CGGTTTTAGgccgcc                                                                                                                                      >  2:905063/1‑147 (MQ=255)                  ttCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑GAAT‑TGGGTTTAGgccgcc                                                                                                                                      >  2:967164/1‑145 (MQ=255)                  ttCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑GAAT‑CGGTTTCAGGCCGCCGa                                                                                                                                    <  2:562822/147‑1 (MQ=255)                   tCAAACTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑TAAT‑TCGGTTCTGGCCGCcgcgcg                                                                                                                                 >  2:770944/1‑149 (MQ=255)                     aaCCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑TAAT‑TGGGTTCAGGGCGCCGGgcggc                                                                                                                               >  1:531205/1‑147 (MQ=255)                             gATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑GAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑a                                                                                                                     <  2:108965/150‑1 (MQ=255)                             gATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑GAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑a                                                                                                                     <  2:144886/150‑1 (MQ=255)                                atTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑GTAT‑GGGGTTCAGGCCGCCCGGCGCCGACAGAAC‑aaa                                                                                                                   >  2:528133/1‑148 (MQ=255)                                atTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑GAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATa                                                                                                                  <  1:1095059/150‑1 (MQ=255)                                       gggTTACTCAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑GAAC‑CGGGTTACGGGCGCCGGGCGCCACCCCGAC‑AATAGCGAtga                                                                                                           >  1:426116/1‑148 (MQ=255)                                        ggTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑TAAT‑CGGGTTTAGGCCGCCGAGCGCCCGCCGAAA‑AAAACCGATaa                                                                                                           >  2:394755/1‑149 (MQ=255)                                            aattaatCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑GAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAA‑TGCCaa                                                                                                    <  2:357047/151‑1 (MQ=255)                                             attaatCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑TAAT‑CGGTTTTAGGCCGCCGAGCGCCCGCAGAAC‑AAAAACGATAA‑AGcc                                                                                                      >  2:195463/1‑148 (MQ=255)                                              ttaatCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑TGAT‑TGGGTTTTGGGCGCCGAGCGCCCGCCGAAA‑AAAAACGATAA‑AGCCAAcc                                                                                                  >  1:902282/1‑151 (MQ=255)                                                aatCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑TTAT‑CGGTTTCAAGGCGCCGAGCCCCCGCAGAAA‑AAAAcccgtaaatcccacc                                                                                                   >  1:633362/1‑136 (MQ=255)                                                aatCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑TAAT‑TTGGTTTAGGCCGCCCAGCGCCCGCAGAACAAAAACCGA‑AA‑TGCCaa                                                                                                    >  1:191045/1‑147 (MQ=255)                                                  tCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑TAAT‑CGGTTTCAGGCCGCCCAGCGCCCGCAGAAC‑AATACCGATAA‑TGCCAACCAGa                                                                                               >  2:50269/1‑150 (MQ=255)                                                   cTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑TAAT‑TGGGTTTAGGCCGCCGGGCGCCCGCAGAAC‑AAAAACGATAA‑TGCCAAACAGa                                                                                               >  1:833776/1‑149 (MQ=255)                                                   cTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑TAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AAAACCGATAA‑TGCCAACCAGaa                                                                                              >  2:350053/1‑150 (MQ=255)                                                     ttCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑TTAT‑TGGTTTCAGGCCGCCGGGCGCCCGCAGAAC‑AATACCGATAA‑TGCCAACCAGAAtt                                                                                            >  1:268948/1‑150 (MQ=255)                                                     ttCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑GGAT‑AGGTTTCAGGCCGCCGAGCCCCGCCAGAAC‑AAAACCGATAATTGCCaaacgaatt                                                                                             >  1:294697/1‑143 (MQ=255)                                                      tCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑GAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAA‑TGCCAACCAGAATTg                                                                                           <  1:682444/150‑1 (MQ=255)                                                      tCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑GAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAA‑TGCCAACCAGAATTg                                                                                           <  2:684352/150‑1 (MQ=255)                                                              tGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑TAAT‑CGGGTTCAAGCCGCCGAGCGCCCGCAGAAC‑AATAACGATAA‑TACCAACCAGAATTTAACaaac                                                                                    >  2:612483/1‑147 (MQ=255)                                                               gTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑GAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAA‑TGCCAACCAGAATTGACCAACCGAt                                                                                 <  2:737788/151‑1 (MQ=255)                                                                   cgcgTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑GAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAA‑TGCCAACCAGAATTGACCAACCGATACg                                                                              <  1:1183826/150‑1 (MQ=255)                                                                    gcgTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGGATT‑AGGTTTC‑GGGCGGCGGGGGGCAGGCGGAC‑AAAACCGTTAA‑TGCCGAACAGaaattaacaaaccata                                                                                >  1:1085068/1‑133 (MQ=255)                                                                    gcgTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑TAAA‑CGGGTTTAAGCCGCCGGGCGCCAGCAGAAA‑AAAAACGATAA‑AGCCCAACAGAAATGACCAACCGATACGc                                                                             >  1:1092392/1‑150 (MQ=255)                                                                    gcgTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑GAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAA‑TGCCAACCAGAATTGACCAACCGATACGCa                                                                            <  1:879987/151‑1 (MQ=255)                                                                     cgTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑GAAT‑CGGGTTCAAGCCGCCGAGCGCCCGCAGAAA‑AAAACCGATAA‑TTCCCAACCGAAATGAACaaacccaacgc                                                                             >  2:322371/1‑140 (MQ=255)                                                                     cgTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGATT‑TAAA‑CGGGTTTAGGCCGCCGAGCGCCAGCAGAAA‑AAAAACGATAA‑TGTCAACCAGAAATTAACAACCCATa                                                                                >  2:405533/1‑146 (MQ=255)                                                                          ttAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑TAAT‑CGGGTTCAGGCCGCCGAGCGCCCGCAGAAC‑AATACCGATAA‑TGCCAACCAGAAATTAACAAACAATACGCAcc                                                                          >  2:297332/1‑147 (MQ=255)                                                                          ttAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑TAAT‑CGGGTTCAGGCCGCCCAGCGCCAGCAGAAC‑AATAACCATAA‑ATCCAACCAAAAATGACCAAACGATACGCaac                                                                          >  1:414238/1‑145 (MQ=255)                                                                          ttAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑GAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAA‑TGCCAACCAGAATTGACCAACCGATACGCACca                                                                         <  2:19937/148‑1 (MQ=255)                                                                            aataaAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑TAATAGGGCTTCAGG‑GGCCGAGGGCGAGCAGAAC‑AAAACAGTGAA‑AGCCAACCAGAAGTGGACAAACGAGACg                                                                              >  1:389828/1‑141 (MQ=255)                                                                              taaAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑TGAT‑AGGGTTCTGGCCGCCGAGCGCCCGCCAAAA‑AAAACCGATAA‑AGCCCACCCAAAATTAACAACCCACACGAACCAGCaaag                                                                   >  2:1166461/1‑150 (MQ=255)                                                                                   gagaGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑GAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAA‑TGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGtt                                                             <  2:283191/151‑1 (MQ=255)                                                                                         gTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑TAAA‑CGGGTTCAGGCCCCCGAGCGCCCCCCGAAC‑AAAAACCGTAA‑AGCCCACCACAAATGACCCACCGAAACGCCCCAACAGAggggggtttttac                                                       >  2:517176/1‑143 (MQ=255)                                                                                            aCCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑TAAT‑CGGGTTTAAGCCGCCCGGCGCCCGCAGAAC‑AAAAACCATAA‑TGCCAACCACAAATGACCAACCGAAACGCCCCACCAAAGGGGGGTCTTcccc                                                      >  1:5186/1‑149 (MQ=255)                                                                                               aTAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑TGAT‑AGGGTTTATGACGCCGGGAGCCAGCCGAAA‑AAAAACGCTAA‑TGCCCAACCCAAATGACCAAACCATACGCACCACCACAGGGGGGTTTTTACCATccc                                                 >  2:709287/1‑149 (MQ=255)                                                                                                      tGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑TGAT‑AGGGTTGAGGGCGCCGGGCGCCCGCAAAAC‑AAAAACGCGAA‑AGACAACCCACATTGAACAAACGACACCCCCCACCCAAGGAGGttctttccc                                                      >  1:193499/1‑135 (MQ=255)                                                                                                       gACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑GAAT‑CGGTTTCAGGCCTCCGAGCGCCAGCAGAAC‑AATACCGATAA‑TGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTc                                          <  1:669779/150‑1 (MQ=255)                                                                                                              aaCAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑TAAA‑CGGTTTTAAGCCGCCCAGGGCCAGCAGGAC‑AATACCGATAA‑AGCCAACCAAAAATTACCAAACGATACCCCCCAACAGAAGGGGGTCTTCCCCCACCCGGACTTCGCgg                                      >  1:182206/1‑147 (MQ=255)                                                                                                               aCAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑TGAT‑AGGGTTCAAGCCGCCGGGGGGCAGCAGAAA‑AAAAACGAAAA‑AGACAAACAAAAATGAACAAACGATAAGCACCAGCAGAGGGGGGTTTTTACCAACGCGGACTTTGCGGGaa                                   >  2:293702/1‑149 (MQ=255)                                                                                                               aCAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑GAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAA‑TGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGt                                 <  2:199385/151‑1 (MQ=255)                                                                                                                cAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑TAAT‑CGGGTTTAGGCCGCCGAGCGCCAGCAGAAC‑AAAACCGATAA‑TGCCAACCAGAATTGACCCAACCATACGCACCAGCAGAGGTGGGTCTTCACCATCCCGTACTTCGCgg                                      >  2:870935/1‑145 (MQ=255)                                                                                                                 aGCGGGCATCCTCCGGCGATAATGGCGGTTT‑GAAA‑CGGTTTTAGGCCGCCCAGCGCCAGCAGAAC‑AATACCGATAA‑TGCCAACCAGAATTGACCAAACGATACGCACCAGCAGAGGGGGGTTTTCACCATCGCGTACTTCGCGGGAAgggg                               >  2:713683/1‑151 (MQ=255)                                                                                                                      gCATCCTCCGGCGATAATGGCGGTTT‑GAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAA‑TGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATc                          <  2:818611/151‑1 (MQ=255)                                                                                                                       cATCCTCCGGCGATAATGGCGGTTT‑TAAT‑CGGGTTCAGGCCCCCCAGCGCCAGCAGAAA‑AAAACCGATAA‑TGACAAACAGAATTGACCAACCCATACGCACCCGCAGAGGGGGGTCTTCCCCATCGCGTACTTCCCGGCAAgggggcctc                          >  1:414239/1‑145 (MQ=255)                                                                                                                       cATCCTCCGGCGATAATGGCGGTTT‑GAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAA‑TGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCg                         <  1:1016498/151‑1 (MQ=255)                                                                                                                        aTCCTCCGGCGATAATGGCGGTTT‑GAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAA‑TGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGc                        <  2:44270/151‑1 (MQ=255)                                                                                                                          cctccGGCGATAATGGCGGTTT‑GAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAA‑TGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGc                        <  1:71406/149‑1 (MQ=255)                                                                                                                               ggCGATAATGGCGGTTT‑GAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAA‑TGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGta                 <  2:493415/151‑1 (MQ=255)                                                                                                                                 cGATAATGGCGGTTT‑GAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAA‑TGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGtata               <  1:142624/151‑1 (MQ=255)                                                                                                                                      aTGGCGGGTT‑GAAT‑CGGGTTACAGCCGCCGAGCGCCAGCAGAAA‑AATACCGATAA‑AGCCAACCAGAAATGACCAACCAAGACGCACCAACAGAGGGGGTTCTTCACCAACGCGGACTTTGCGGCAACCGGACATCGCCAGGGTATAaga            >  1:886349/1‑149 (MQ=255)                                                                                                                                        ggCGGTTT‑GAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAA‑TGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCa         <  1:139330/150‑1 (MQ=255)                                                                                                                                        ggCGGTTT‑GAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAA‑TGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCa         <  1:1109039/150‑1 (MQ=255)                                                                                                                                        ggCGGTTT‑GAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAA‑TGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAg        <  1:353370/151‑1 (MQ=255)                                                                                                                                               t‑GAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAA‑TGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTa  <  1:604978/150‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                        TAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑GAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAA‑TGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTA  >  NC_000913/40372‑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: ‑