Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	40,522	0	T	G	42.2%	101.3 / 13.3	65	K470N (AAA→AAC)	caiT	putative transporter
Reads supporting (aligned to +/- strand):  ref base T (12/25);  new base G (27/0);  total (40/25)
Fisher's exact test for biased strand distribution p-value = 6.21e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑GAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAA‑TGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGT  >  NC_000913/40373‑40671                                                                                                                                                        |                                                                                                                                                        aaattaaaATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑GAAT‑CGGtt                                                                                                                                                        <  1:281599/151‑1 (MQ=255)  aattaaaATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑GAAT‑CGGttt                                                                                                                                                       <  1:392808/151‑1 (MQ=255)  aattaaaATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑GAAT‑CGGttt                                                                                                                                                       <  1:819571/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:144886/150‑1 (MQ=255)                            gATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑GAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑a                                                                                                                             <  2:108965/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                                                                                                   <  2:684352/150‑1 (MQ=255)                                                     tCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑GAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAA‑TGCCAACCAGAATTg                                                                                                   <  1:682444/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)                                                                                                                                                   t‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAA‑TGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAAcc       >  1:1137084/1‑149 (MQ=255)                                                                                                                                                     cGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCTATAA‑TGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGc    <  1:326932/151‑1 (MQ=255)                                                                                                                                                       gTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAA‑TGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGt  >  2:574615/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                        AAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTT‑GAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAAC‑AATACCGATAA‑TGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGT  >  NC_000913/40373‑40671

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