Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	40,515	0	G	T	41.8%	72.7 / 32.3	68	Q473K (CAA→AAA)	caiT	putative transporter
Reads supporting (aligned to +/- strand):  ref base G (9/30);  new base T (28/0);  total (38/30)
Fisher's exact test for biased strand distribution p-value = 2.44e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

CATTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCG  >  NC_000913/40368‑40660                                                                                                                                                      |                                                                                                                                                   catTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAAt                                                                                                                                                 <  1:459199/151‑1 (MQ=255)     taaattaaaATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCAAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑CTTTaagacgg                                                                                                                                             >  1:929585/1‑146 (MQ=255)       aattaaaATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTTAAT‑CGGttt                                                                                                                                          >  1:810928/1‑151 (MQ=255)         ttaaaATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTTAAT‑Tgggt                                                                                                                                           >  2:345164/1‑146 (MQ=255)          taaaATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGTAT‑CGGgtt                                                                                                                                          >  2:708822/1‑148 (MQ=255)          taaaATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAAT‑CGGTTTCa                                                                                                                                        <  1:633060/150‑1 (MQ=255)           aaaaTCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑GTTGAAT‑CGGTTTTAgg                                                                                                                                      >  1:114437/1‑151 (MQ=255)            aaaTCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAAT‑GGGattta                                                                                                                                        >  1:887287/1‑146 (MQ=255)             aaTCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAAT‑CGGTTTCAGgcc                                                                                                                                    <  1:534728/151‑1 (MQ=255)              aTCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAAT‑CGGTTTCAGgccg                                                                                                                                   <  1:1037798/151‑1 (MQ=255)                ccATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAAT‑CGGTTTCAGgccgcc                                                                                                                                 <  2:513639/151‑1 (MQ=255)                       tCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAAT‑CGGGTTCAGGCCGCCCAGCGcc                                                                                                                          >  1:29051/1‑151 (MQ=255)                        cAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTTTAT‑TGGGTTTAGGGCGCCGAGCGcc                                                                                                                          >  2:594799/1‑150 (MQ=255)                             tctTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTTAAT‑AGGGTTTAGggcgcccggcggc                                                                                                                          >  2:446316/1‑138 (MQ=255)                              ctTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGG‑TTTGAA‑‑CGGGTTCCGGCCGCCCCGCGcccccc                                                                                                                      >  2:295014/1‑144 (MQ=255)                              ctTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTTTAA‑CGGTTTCATGGCGCCCGGCGCcagca                                                                                                                      >  1:247175/1‑148 (MQ=255)                              ctTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTTAAT‑AGGGTTCAgggggcc                                                                                                                                 >  2:1091576/1‑132 (MQ=255)                              ctTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCAGGCGATAATGGC‑GG‑TTTGAAT‑CGGTTTCAGGCCGCCGAGCGCcagcag                                                                                                                     <  1:391707/149‑1 (MQ=255)                                 gATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAACa                                                                                                                 <  2:116358/150‑1 (MQ=255)                                      tttGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTTAAT‑TGGGTTTAAGCCGCCGAGCGCcagccggaaaa                                                                                                                >  1:271971/1‑139 (MQ=255)                                      tttGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATAcc                                                                                                            <  2:610399/150‑1 (MQ=255)                                       ttGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTTGAT‑AGGGTTTTGGCCGCCGAGCCCCCGCCCAACAATaacc                                                                                                           >  2:994220/1‑147 (MQ=255)                                        tGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTTGAT‑TGGGTTTAGGCCGCCGCGCGCCCGCCGAaaaaa                                                                                                               >  2:251971/1‑144 (MQ=255)                                        tGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑GTTTAAT‑TGGGTTCAAGGCGCCCCGCGCCCGCAGAaaaaaaac                                                                                                            >  1:927793/1‑144 (MQ=255)                                           gtgTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTTAAT‑TGGGTTTAGGCCGCCCGGCGCCCGCCGAACAAAAACGATaa                                                                                                       >  1:854480/3‑150 (MQ=255)                                                      tCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GGTTTTGAAT‑‑GGTTTTTGGCCGCCGAGCGCCCCCAGAAAAAAAACCATAATGccca                                                                                                 >  2:581602/1‑143 (MQ=255)                                                      tCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTTTAT‑CGGGTTTACGCCCCCCCGCGCCCGCAGAACAATACCGATAATGCCAAcc                                                                                               >  1:994937/1‑147 (MQ=255)                                                         ttCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGaa                                                                                           <  1:929445/148‑1 (MQ=255)                                                           ccAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTTAAT‑TGGGTTTCTGCCGGCCCGCGCCCGCAGAAAAATAACCATAAtatca                                                                                                  >  2:239945/1‑135 (MQ=255)                                                           ccAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAAA‑CGGTTTTAGGGCGCCGAGCGCCCGCAGAACAAAACCGATAATGCCAACCAGaaat                                                                                         >  1:1007410/1‑146 (MQ=255)                                                             aGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCa                                                                                    <  2:100344/151‑1 (MQ=255)                                                              gTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTTGAT‑TGGGTTCAGGCCGCCCAGGGCCAGCAGAACAATACCCACAATGCCCAACCGAAATGaac                                                                                     >  1:242601/1‑147 (MQ=255)                                                               ttCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTTAAT‑CGGTTTTACGGCGCCCCGCGGCCGCCGAAAAAAACCGCAAATGCCAAACAAAATTGaac                                                                                     >  2:630431/1‑146 (MQ=255)                                                                  tGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GGTTTTAAAT‑GGGTTTAGGGCCCCCG‑GCGCCCGCAGAAAAAAAACGAAAATGCCAACCaaaaattaacaaac                                                                                 >  1:600159/1‑139 (MQ=255)                                                                   gTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAAcc                                                                                 <  1:461246/148‑1 (MQ=255)                                                                   gTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAAcc                                                                                 <  1:355507/148‑1 (MQ=255)                                                                   gTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGa                                                                               <  1:1092226/150‑1 (MQ=255)                                                                    tttCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTTAAT‑CGGGTTTAGGCCGCCGCGCGCCCGCAGAAAAAAAACGATAAAGCCAACCAGAAATTAcccaac                                                                                 >  1:661194/1‑143 (MQ=255)                                                                     ttCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAATACGGTTTAAGG‑CGGCGCGCGCCCGCCGGAAAAAACCGAAAATGCCAACCAGAAATGACCAACCGaca                                                                             >  1:127641/1‑148 (MQ=255)                                                                       cgcgTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTTAAT‑CGGGGTTAGGCCGCCCCGCGCCCGCCGAAAAAAACCGATAATATCAACCAGAAATGACCAAACGATa                                                                             >  1:1041855/1‑148 (MQ=255)                                                                            cTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTTGAA‑CCGTTTTTGGCCGGCGGGAGCCCGCCAAAAAAAAACAAAAAAACCAACCACAAATGACCAAACGAAACCCAccc                                                                      >  1:389963/1‑149 (MQ=255)                                                                              ttAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCGCcagc                                                                    <  1:1018306/150‑1 (MQ=255)                                                                              ttAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACcagc                                                                    <  2:181603/150‑1 (MQ=255)                                                                                aataaAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTTAAT‑CGGTTTTATGCCGCCCAGCGCCCGCAGAAAAAAAACCATAATTCCAACCAAAAATGACCAACCCATaacc                                                                          >  1:774216/1‑139 (MQ=255)                                                                                aataaAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTTAAT‑AGGGTTTAAGCCGCCCAGCGCCAACAGAAAAAAACCCATAATGCCAAACAGAAATTAACAAACGAAAAGCAcccac                                                                    >  1:569238/1‑145 (MQ=255)                                                                                  taaAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTTAAT‑CGGGTTCAGGCCGCCCGGCGCCCGCCGAAAAAAAACGATAAATCCAACCCAAATTGACCAACCGATACGCAcccacaaagg                                                               >  1:430263/1‑143 (MQ=255)                                                                                     aGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTTTAT‑TGGGTTTAGGCCGCCCCGCGCCCGCAGAACAAAACCGACAAAGGCAACCAAAAATGTCCAACCGGAACGCACCAGCCGAggggg                                                            >  2:500353/1‑151 (MQ=255)                                                                                        agagCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGttc                                                         <  2:950236/151‑1 (MQ=255)                                                                                            cGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCATCCGATACGCACCAGCAGAGGTGGTTCTTCa                                                     <  1:583841/151‑1 (MQ=255)                                                                                                aCCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTTAAT‑CGGTTTTTCGCCCCCCAGGGCCCGCAGAACAAAAACGATAAAGACAACCAGAAATTACCCACCCAAACGCCCCCGCAAAggggggtttttccc                                                   >  1:295797/1‑140 (MQ=255)                                                                                                 ccATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTTAAT‑AGGTTTCAAGCCGCCGAGCGCCCGCCGAACAAAAACCAAAAAGCCAAACAGAAATTACCAACCCATACCCACCCGCAGAGGGGGGTTTTCAcc                                                   >  1:523072/1‑148 (MQ=255)                                                                                                 ccATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATc                                                <  2:5725/151‑1 (MQ=255)                                                                                                  cATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTTAAT‑CGGGTTTTGGCCGCCGGGCGCCCGCCGAACAAAACCGATAATGCCCACCAGAAATTACCAAACGATACGCACCAGCAGAGGGGGGTTTTTAcc                                                   >  1:208937/1‑147 (MQ=255)                                                                                                    tAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTTGAT‑CGGTTTTAGGGCGCCGGGCGCCCGCCGAAAAATAACGATTATACCAAACAGAAATGAACCACCCATACGCACCCCCAGAGGGGGgtcttcctccaccc                                              >  2:19609/1‑142 (MQ=255)                                                                                                          tGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACtt                                        <  1:937466/150‑1 (MQ=255)                                                                                                           gACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTcg                                      <  1:389059/151‑1 (MQ=255)                                                                                                                  aaCAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTTTAA‑AGGTTTCCGGCCGCCGGGGGCCCGCAGAAAAAAAACGATAATATCAACACGAAATGACACACCGCTACGGGCCCCCAGCGGGGGGTCTTCACCCCCGCGCCCTTcgcg                                    >  2:58573/1‑146 (MQ=255)                                                                                                                         ggCATCCTCCGGCGATAATGGC‑GG‑TTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACAt                        <  2:494985/151‑1 (MQ=255)                                                                                                                          gCATCCTCCGGCGATAATGGC‑GG‑TTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATc                       <  1:251976/151‑1 (MQ=255)                                                                                                                           cATCCTCCGGCGATAATGGC‑GG‑TTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCg                      <  1:446329/151‑1 (MQ=255)                                                                                                                           cATCCTCCGGCGATAATGGC‑GG‑TTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCg                      <  1:522236/151‑1 (MQ=255)                                                                                                                                 ccGGCGATAATGGC‑GG‑TTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAggg                <  1:848242/151‑1 (MQ=255)                                                                                                                                 ccGGCGATAATGGC‑GG‑TTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAggg                <  2:764860/151‑1 (MQ=255)                                                                                                                                 ccGGCGATAATGGC‑GG‑TTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAgg                 <  1:41244/150‑1 (MQ=255)                                                                                                                                   ggCGATAATGGC‑GG‑TTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGta              <  2:723488/151‑1 (MQ=255)                                                                                                                                   ggCGATAATGGC‑GG‑TTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGta              <  2:814368/151‑1 (MQ=255)                                                                                                                                    gCGATAATGGC‑GG‑TTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGtat             <  2:1029454/151‑1 (MQ=255)                                                                                                                                                gg‑TTTGACT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCg  >  2:11847/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                   CATTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGC‑GG‑TTTGAAT‑CGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCG  >  NC_000913/40368‑40660

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