Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,538,109	0	A	C	42.4%	104.3 / 16.4	60	G219G (GGT→GGG)	nanM	N‑acetylneuraminic acid mutarotase
Reads supporting (aligned to +/- strand):  ref base A (31/3);  new base C (0/25);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 1.08e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCG‑CAATCCT‑GG‑TTTGGCTTCGCCATTAATAAG‑CCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTT  >  NC_000913/4537960‑4538258                                                                                                                                                           |                                                                                                                                                      aTCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCG‑CAATCCT‑GG‑TTTGGCTTCGCCATTAATAAG‑CCAGG‑TTTTATCAc                                                                                                                                                      <  1:857203/151‑1 (MQ=255) aTCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCG‑CAATCCT‑GG‑TTTGGCTTCGCCATTAATAAG‑CCAGG‑ATTTATCAc                                                                                                                                                      >  1:336846/1‑151 (MQ=255)          aaaCCGCCAGCTACGCCATCTGGTGATGAGAAGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCG‑CAATCCT‑GG‑TTTGGCTTCGCCATTAATAAG‑CCAGG‑TTTTATCACCTTTATTc                                                                                                                                              >  1:555872/1‑150 (MQ=255)            aCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCG‑CAATCCT‑GG‑TTTGGCTTCGCCATTAATAAG‑CCAGG‑TTTTATCACCTTTATTcaca                                                                                                                                           <  1:277339/151‑1 (MQ=255)             ccgccgGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCG‑CAATCCT‑GG‑TTTGGCTTCGCCATTAATAAG‑CCAGG‑TTTTATCACCTTTATTcaca                                                                                                                                           >  2:399356/1‑150 (MQ=255)             ccGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCG‑CAATCCT‑GG‑TTTGGCTTCGCCATTAATAAG‑CCAGG‑TTTTATCACCTTTATTCACaa                                                                                                                                          <  2:586174/151‑1 (MQ=255)                ccAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCG‑CAATCCT‑GG‑TTTGGCTTCGCCATTAATAAG‑CCAGG‑TTTTATCACCTTTATTCACAAccg                                                                                                                                       >  1:621102/1‑151 (MQ=255)                 cAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCG‑CAATCCT‑GG‑TTTGGCTTCGCCATTAATAAG‑CCAGG‑TTTTATCACCTTTATTCACAAccgc                                                                                                                                      >  2:760227/1‑151 (MQ=255)                      aCGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCG‑CAATCCT‑GG‑TTTGGCTTCGCCATTAATAAG‑CCAGG‑TTTTATCACCTTTATTCACAACCGCCGCAc                                                                                                                                 >  2:813034/1‑151 (MQ=255)                        gCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCG‑CAATCCT‑GG‑TTTGGCTTCGCCATTAATAAG‑CCAGG‑TTTTATCACCTTTATTCACAACCGCCGCAcc                                                                                                                                >  2:813234/1‑150 (MQ=255)                            tCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCG‑CAATCCT‑GG‑TTTGGCTTCGCCATTAATAAG‑CCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCg                                                                                                                           >  1:404915/1‑151 (MQ=255)                                gTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCG‑CAATCCT‑GG‑TTTGGCTTCGCCATTAATAAG‑CCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTcc                                                                                                                       >  1:965250/1‑151 (MQ=255)                                                aaGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCG‑CAATCCT‑GG‑TTTGGCTTCGCCATTAATAAG‑CCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCg                                                                                                       >  1:34641/1‑151 (MQ=255)                                                 aGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCG‑CAATCCT‑GG‑TTTGGCTTCGCCATTAATAAG‑CCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGc                                                                                                      >  1:151226/1‑151 (MQ=255)                                                      aTTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCG‑CAATCCT‑GG‑TTTGGCTTCGCCATTAATAAG‑CCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCg                                                                                                 >  1:196506/1‑151 (MQ=255)                                                      aTTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCG‑CAATCCT‑GG‑TTTGGCTTCGCCATTAATAAG‑CCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCg                                                                                                 >  2:817607/1‑151 (MQ=255)                                                            ttttAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCG‑CAATCCT‑GG‑TTTGGCTTCGCCATTAATAAG‑CCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTc                                                                                           >  2:806785/1‑151 (MQ=255)                                                            tttaAAATTTTTCCCGGT‑AAAACAAATTAAAAATCCGCCTCCCTTCC‑CCACTCT‑GG‑TTTTGCCTCCCCCATAAAAAG‑GCAGGTTTTTACCCCTTTTATCCCCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAAc                                                                                             <  2:671728/149‑1 (MQ=255)                                                                aaaTTTTTCCCGGTAAAATCAAGTTAAAATCCGGATTCCGTCCC‑AAATCTT‑GT‑TTTGGCTCCCCAATTAATAAC‑CCGGT‑TTTTATCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCa                                                                                         <  2:553401/149‑1 (MQ=255)                                                                aaaTTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCG‑CAATCCT‑GG‑TTTGGCTTCGCCATTAATAAG‑CCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCAtt                                                                                       >  2:972316/1‑151 (MQ=255)                                                                   ttattaCCGGTGAAATCAAGTTCAAATACGGCATCCGTTCG‑CAATCCT‑GG‑TTTGGCTTCGCCATTAATAAG‑CCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCt                                                                                    >  1:153658/1‑151 (MQ=255)                                                                     attaCCGGTGAAATCAAGTTCAAATACGGCATCCGTTCG‑CAATCCT‑GG‑TTTGGCTTCGCCATTAATAAG‑CCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtgt                                                                                  >  2:293574/1‑151 (MQ=255)                                                                        aCCGGTGAAATCAAGTTCAAATACGGCATCCGTTCG‑CAATCCT‑GG‑TTTGGCTTCGCCATTAATAAG‑CCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGtt                                                                               >  1:949451/1‑151 (MQ=255)                                                                         ccGGTTAAAACAAATTTAAAATCCGCCTCCCTTCCCCAACCCT‑GT‑TTTGTCTCCCCCATAAATAAG‑CCAGT‑TTTTACCCCTTTTATCCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGtt                                                                               <  2:779538/151‑1 (MQ=255)                                                                            gtaaaaTAAAGTTAAAATACGGCATCCGTTCG‑CAATCCT‑GG‑TTTGGCTTCCCCATTAATAAG‑CCAGT‑TTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAgg                                                                           <  1:919284/148‑1 (MQ=255)                                                                             taaaaTCAAGTTCAAATCCGGCTCCCTTCCG‑CAACCCG‑GT‑TTTGGTTCCGCTATTAAAAAG‑CCTGG‑TTTTATCTCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGa                                                                             <  1:944808/146‑1 (MQ=255)                                                                              aaaacaaatttaaaaTGCGCCACCCGTCCG‑AAACCCT‑GT‑TTTGGCTCCCCCATTAATAAC‑CCAGT‑TTTTACCCCCTTTATCCCCACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAgg                                                                           <  2:898882/137‑1 (MQ=255)                                                                               aaaCCAAGTTCAAATACGGCCTCCGTTCG‑CAACCCT‑GT‑TTTGTCTTCCCCTTTAAAAAC‑CCGGG‑TTTTATCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGa                                                                         <  1:387482/150‑1 (MQ=255)                                                                                aataaattacAAATCCGCCACCCGTTCG‑AACTCCG‑GT‑TTTGGTTCCCCCATAAAAAAC‑CCAGT‑TTTTACCCCTTTAATCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATc                                                                       <  2:604745/142‑1 (MQ=255)                                                                                aagaaatttCAAATCCGGCACCGGTTCC‑CAATCCG‑GG‑TTGGGCTTCGCC‑TTAAAAAGCCCAGG‑TTTTATCCCCTTTATCCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATc                                                                       <  2:11534/144‑1 (MQ=255)                                                                                  tCAAGTTCAAATACGGCATCCGTTCG‑CAATCCT‑GG‑TTTGGCTTCGCCATTAATAAG‑CCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCa                                                                      >  2:417934/1‑150 (MQ=255)                                                                                       ttaaaaTACGCCCTCC‑TTCC‑CCACCCT‑GGTTTTGGCTTCCCATTTAAAAAC‑CCAGT‑TTTTATCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGa                                                                 <  2:288987/147‑1 (MQ=255)                                                                                       ttCAAATACGGCATCCGTTCG‑CAATCCT‑GG‑TTTGGCTTCGCCATTAATAAG‑CCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGac                                                                >  2:391568/1‑151 (MQ=255)                                                                                          aaaTAGGGCTTCCTTTCC‑CAATCCT‑GT‑TTTGGTTTCCCAATTAATAAG‑CCAGG‑TTTTATCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaac                                                             <  2:936067/151‑1 (MQ=255)                                                                                          aaaTACGGCATCCGTTCG‑CAATCCT‑GG‑TTTGGCTTCGCCATTAATAAG‑CCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaac                                                             >  1:694157/1‑151 (MQ=255)                                                                                               cgccATCCTTTCG‑C‑ATCCTGGG‑TTTGGCTCCGCCATTAATAAG‑CCGGT‑TTTTTCCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGaaa                                                        <  2:715246/148‑1 (MQ=255)                                                                                                ggatccctttCG‑CAACCTT‑GG‑TTTGGCTCCGCAATTAATAAG‑CCAGT‑TTTTACCCCCTTTATTCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGaaa                                                        <  1:223122/142‑1 (MQ=255)                                                                                                  cATCCGTTCG‑CAATCCT‑GG‑TTTGGCTTCGCCATTAATAAG‑CCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACtt                                                     >  1:596956/1‑151 (MQ=255)                                                                                                              aacccT‑GT‑TTGGCCTTCGCATTAAATAAG‑CCAGT‑TTTTTCCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaag                                              <  2:540367/144‑1 (MQ=255)                                                                                                                   cttg‑TTTGGCTTCCCAATTAATAAC‑CCAGT‑TTTTAACCCCTTTCTCCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaaga                                             <  2:476699/141‑1 (MQ=255)                                                                                                                     gg‑TTGGGTTTCCCCATTAAAAAG‑CCAGT‑TTTTACCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaat                                          <  2:351545/145‑1 (MQ=255)                                                                                                                      g‑TTTGGCTTCGCCATTAATAAG‑CCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtcttc                                   >  2:897498/1‑151 (MQ=255)                                                                                                                        tttGGCTTCGCCATAAATAAG‑CCAGT‑TTTTACCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtc                                      <  1:314850/147‑1 (MQ=255)                                                                                                                         ttGGCTTCGCCATTAATAAG‑CCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTg                                 >  2:957607/1‑151 (MQ=255)                                                                                                                          tGGCTTCGCCATTAATAAG‑CCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGc                                >  2:497334/1‑151 (MQ=255)                                                                                                                           ggccccGCCATAAAAAAG‑CCAGT‑TTTTACCCCCTTAATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTg                                 <  1:837956/144‑1 (MQ=255)                                                                                                                                cGCAATAAATAAG‑CCGGT‑TTTTATCCCCTTTATACACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCttt                             <  2:152431/148‑1 (MQ=255)                                                                                                                                 gCCATTAATAAG‑CCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTg                         >  1:286275/1‑151 (MQ=255)                                                                                                                                 ccctttAATAAC‑CATGT‑TTTTATCCCCTTTATTCCAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTg                         <  2:865100/147‑1 (MQ=255)                                                                                                                                     ttAATAAG‑CCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCa                      >  1:853482/1‑150 (MQ=255)                                                                                                                                        aatactccAGT‑TTTTTCCCCCTTTATCCCAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaa                     <  1:138301/143‑1 (MQ=255)                                                                                                                                          aaG‑CCAGT‑TTTTATCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAg                <  1:377631/151‑1 (MQ=255)                                                                                                                                            aatcgtgtttttATCCCTTTTATTCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAg                <  2:778872/143‑1 (MQ=255)                                                                                                                                                    ttttATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCAtt        >  1:19397/1‑151 (MQ=255)                                                                                                                                                    ttttACCCCCTTTTTCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCAt         <  2:317532/150‑1 (MQ=255)                                                                                                                                                    ttttACCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCAt         <  1:896781/150‑1 (MQ=255)                                                                                                                                                      ttATCCCCTTTATTCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTg       <  1:78034/150‑1 (MQ=255)                                                                                                                                                      ttATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTg       >  1:238486/1‑150 (MQ=255)                                                                                                                                                      ttATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTg       >  1:601275/1‑150 (MQ=255)                                                                                                                                                          cACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGAtttt  >  1:512770/1‑151 (MQ=255)                                                                                                                                                           |                                                                                                                                                      ATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCG‑CAATCCT‑GG‑TTTGGCTTCGCCATTAATAAG‑CCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTT  >  NC_000913/4537960‑4538258

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