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

ATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCCAGG‑TTTT‑AT‑CACCTTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTT  >  NC_000913/4537960‑4538258                                                                                                                                                           |                                                                                                                                                       attcccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCCAGG‑TTTT‑AT‑CAc                                                                                                                                                       >  2:987923/4‑151 (MQ=255)    cccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCCAGG‑TTTT‑AT‑CACCtt                                                                                                                                                    >  1:784224/1‑151 (MQ=255)      cGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCCAGG‑TTTT‑AT‑CACCTTT‑a                                                                                                                                                 <  1:985925/151‑1 (MQ=255)       gCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCCAGG‑TTTT‑AT‑CACCTTT‑At                                                                                                                                                <  2:689282/151‑1 (MQ=255)        cAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCCAGG‑TTTT‑AT‑CACCTTT‑Att                                                                                                                                               >  2:942607/1‑151 (MQ=255)         aaaaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCCAGG‑TTTT‑AT‑CACCTTT‑ATTc                                                                                                                                              <  2:967026/151‑1 (MQ=255)             ccGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCCAGG‑TTTT‑AT‑CACCTTT‑ATTcaca                                                                                                                                           >  2:682145/1‑150 (MQ=255)             ccGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCCAGG‑TTTT‑AT‑CACCTTT‑ATTCACaa                                                                                                                                          <  2:677906/151‑1 (MQ=255)              cGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCCAGG‑TTTT‑AT‑CACCTTT‑ATTCACAAc                                                                                                                                         <  2:700297/151‑1 (MQ=255)               gCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCCAGG‑TTTT‑AT‑CACCTTT‑ATTCACAAcc                                                                                                                                        <  2:604517/151‑1 (MQ=255)                  aGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCCAGG‑TTTT‑AT‑CACCTTT‑ATTCACAAccgcc                                                                                                                                     >  2:891521/1‑151 (MQ=255)                      aCGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCCAGG‑TTTT‑AT‑CACCTTT‑ATTCACAACCGCCGCAc                                                                                                                                 <  2:97689/151‑1 (MQ=255)                         ccATCTGGTGATGAGACGGGACCAAGTTTATTCCATTTTAATTTATTCCCGGTAAAATCAATTTCAAATACGGCATCCTTTCGCAATCCTGT‑TTTGGCTT‑CGCCATAAATAACCCAGG‑TTTT‑AT‑CCCCTTT‑ATTCACAACCGCCGCACCa                                                                                                                               <  1:219727/150‑1 (MQ=255)                               ggTGATGAGACGGGAGCAAGTTTATTCCTTTTTAAATTATTCCCGGTAAAATCAAGTTCAAATCCGGCTTCCGTTCGCACTCCTGT‑TTTGGCTT‑CGCCATTAATAAGCCGGG‑TTTT‑AT‑CCCCTTT‑ATCCACAACCGCCGCACCAGc                                                                                                                             <  1:796755/146‑1 (MQ=255)                                      agaCGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCCAGG‑TTTT‑AT‑CACCTTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTAcc                                                                                                                  >  1:672400/1‑150 (MQ=255)                                          gggacaaatcTTTTTCCATTTTAAATTATACCCGGTAAAATAAAGTTAAAATACGGCACCCTTCCCCAACCCTGG‑TTGGGCTC‑CGCAATAAATAAGCCAGT‑TTTT‑AT‑CCCCTTT‑ATCCAAAACCGCCGCACCAGCCGTTCCGTACCAg                                                                                                                <  2:816975/139‑1 (MQ=255)                                                aaGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCCAGG‑TTTT‑AT‑CACCTTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCg                                                                                                       >  1:733341/1‑151 (MQ=255)                                                   cTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCCAGG‑TTTT‑AT‑CACCTTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc                                                                                                     >  1:118228/1‑150 (MQ=255)                                                    ttATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCCAGG‑TTTT‑AT‑CACCTTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc                                                                                                     >  2:167214/1‑149 (MQ=255)                                                        tCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCCAGG‑TTTT‑AT‑CACCTTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGACAGCGt                                                                                                >  1:202807/1‑150 (MQ=255)                                                         ccATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCCAGG‑TTTT‑AT‑CACCTTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTa                                                                                               >  2:401785/1‑150 (MQ=255)                                                           aTTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCCAGG‑TTTT‑AT‑CACCTTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACt                                                                                            >  1:12079/1‑151 (MQ=255)                                                            ttttAAATTATTCCCGGTAAAACCAAGTTAAAATACGGCTTCTTTCCGCAATCCTGT‑TTGGGCTT‑CGCCATAAATAAGCAGGT‑TTTT‑AT‑CCCTTTT‑ATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACt                                                                                            <  2:800502/150‑1 (MQ=255)                                                            tttaAAATTATTCCCGGTAAAATCAAGTTAAAACCCGGCTCCCTTTCGCAATCCTGG‑TTTGGCTT‑CCCCTTTAAAAAGCCAGT‑TTTT‑CC‑CCCCTTT‑ATCCCCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGc                                                                                                  <  1:345479/144‑1 (MQ=255)                                                            ttaaaaaTTATTCCCGGAAAAATAAATTAAAATTCCGGCACCCGTCCCAATTCCTGT‑TTTGCTTT‑CCCCTTTAAAAACCCGGT‑TTTT‑AC‑CCCCTTT‑ATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTaa                                                                                              <  1:8423/144‑1 (MQ=255)                                                                aaaTTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCCAGG‑TTTT‑AT‑CACCTTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCAtt                                                                                       >  1:226450/1‑151 (MQ=255)                                                                   ttattaCCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCCAGG‑TTTT‑AT‑CACCTTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCt                                                                                    >  1:583116/1‑151 (MQ=255)                                                                     attaCCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCCAGG‑TTTT‑AT‑CACCTTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtgt                                                                                  >  1:657830/1‑151 (MQ=255)                                                                         cggttaaaaaaaaGTACAAATCCGCCACCCTTTCGAAATCCGGG‑TTTGCCTT‑CCCCATAAATACGCCGGT‑TTTT‑AC‑CCCCTTT‑ATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGtt                                                                               <  2:787933/139‑1 (MQ=255)                                                                         ccGGTAAAATAAAGTTAAAATCCGGCCCCCGTTCGCAATCCTGT‑TTTGGCTC‑CGCAATTAATAAGCCAGT‑TTTT‑AC‑CCCCTTT‑ATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtgtg                                                                                 <  2:1037765/148‑1 (MQ=255)                                                                               aaataaatTTAAAATACGACATCCGTTCGCAATCCTGT‑TTTGGTTT‑CGCCATAAATAACCCAGT‑TTTT‑AC‑CCCTTTT‑ATTCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAggg                                                                          <  2:117834/149‑1 (MQ=255)                                                                               aaaTCAAGTTCAAATACGGCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCCAGG‑TTTT‑AT‑CACCTTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGAt                                                                        >  1:969559/1‑151 (MQ=255)                                                                                 aTCAAGTTCAAATACGGCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCCAGG‑TTTT‑AT‑CACCTTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATc                                                                       >  1:497738/1‑150 (MQ=255)                                                                                   cAAGTTCAAATACGGCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCCAGG‑TTTT‑AT‑CACCTTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaa                                                                     >  1:749643/1‑150 (MQ=255)                                                                                    aatttaaaaTGCGGCACCCGTCCGCAATCCTGG‑TTTGCCTC‑CCCCATAAAAAACCCAGT‑TTTT‑AT‑CCCTTTT‑ATCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaaaa                                                                   <  2:339386/145‑1 (MQ=255)                                                                                    aaGTTCAAATACGGCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCCAGG‑TTTT‑AT‑CACCTTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaaa                                                                    >  1:555803/1‑150 (MQ=255)                                                                                          aaaCACGGCCCCCGTCCCCATTCCTTG‑TTTGGCTC‑CGCATTAAATAACCCAGT‑TTTT‑AC‑CCCCTTT‑TTACACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAg                                                                  <  1:767363/146‑1 (MQ=255)                                                                                             tACGGCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCCAGG‑TTTT‑AT‑CACCTTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaac                                                             >  2:881752/1‑148 (MQ=255)                                                                                               cGGCATCCTTTCGCAATCCGGG‑TTTGGCTT‑CGCCATAAAAAAGCCAGT‑TTTT‑AC‑CCCCTTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGaa                                                         <  2:77111/150‑1 (MQ=255)                                                                                                  cacccGTTCGCACCCCTGT‑TTAGGCTT‑CCCCTTTAAAAAGCCAGT‑TTTT‑TT‑CCCTTTT‑ATCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACtt                                                     <  2:234383/148‑1 (MQ=255)                                                                                                     cgtttCGCACCCCTGG‑TTTGTCTT‑CCCCTTAAAAAACCCGGT‑TTTT‑AC‑CCCTTTT‑ATCCCCATCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACt                                                      <  2:404722/144‑1 (MQ=255)                                                                                                     ccttccccatcccTGT‑TTTGGCTC‑CCCCATAAATAACCCGGT‑TTTT‑AC‑CCCTTTT‑ATCCAAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTAtt                                                  <  1:223652/140‑1 (MQ=255)                                                                                                     ccGTCCGCAACCCTGT‑TTTGCCTC‑CCCCATAAAAAAGCCAGT‑TTTT‑AC‑CCCTTTT‑ATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACt                                                      <  2:161449/147‑1 (MQ=255)                                                                                                      cGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCCAGG‑TTTT‑AT‑CACCTTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTg                                                 >  1:715970/1‑151 (MQ=255)                                                                                                      cGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCCAGG‑TTTT‑AT‑CACCTTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTg                                                 >  2:760839/1‑151 (MQ=255)                                                                                                        ttCGCAATCCTGT‑TTTGCCTT‑CGCCATTAAAAAGCCAGT‑TTTT‑AT‑CCCTTTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTga                                                <  2:727844/150‑1 (MQ=255)                                                                                                        ttCGCAATCCTGG‑TTTGGTTT‑CCCCATTAATAAGCCAGT‑TTTT‑AT‑CCCTTTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTAtt                                                  <  2:360879/148‑1 (MQ=255)                                                                                                          cGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCCAGG‑TTTT‑AT‑CACCTTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTga                                                >  2:130683/1‑148 (MQ=255)                                                                                                           cgaaTTCTTG‑TTTGGGTTCCCCCATTAA‑AAGCCCGGTTTTT‑AT‑CCCTTTT‑ATCCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTg                                                 <  1:343475/145‑1 (MQ=255)                                                                                                             aaTCCTGG‑TTTGGCTT‑CGCCATTAATAAGCCAGG‑TTTT‑AT‑CACCTTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaat                                          >  1:291534/1‑151 (MQ=255)                                                                                                              aTCCTGGTTTTGGTTC‑CGCAATAAAAAACCCAGT‑TTTT‑AC‑CCCTTTT‑ATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaa                                               <  1:520836/146‑1 (MQ=255)                                                                                                                 cTGG‑TTTGGCTT‑CGCCATTAATAAGCCAGG‑TTTT‑AT‑CACCTTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaataat                                       >  2:40045/1‑150 (MQ=255)                                                                                                                    gttttgcttt‑CGCAATTAATAAGCAAGT‑TTTT‑AT‑CCCCTTT‑ATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtct                                     <  2:705560/142‑1 (MQ=255)                                                                                                                    g‑TTTGGCTT‑CGCCATTAATAAGCCAGG‑TTTT‑AT‑CACCTTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtcttc                                   >  2:1014948/1‑151 (MQ=255)                                                                                                                      tttGCCTC‑CGCCA‑TAATAAACCAAG‑TTTT‑AT‑CCCCTTTAATTCCCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaata                                         <  1:507407/144‑1 (MQ=255)                                                                                                                       ttGGCTT‑CGCAATTAATAGGCCAGT‑TTTT‑AT‑CACCTTT‑ATCCACAACCCCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtcttct                                  <  1:600311/150‑1 (MQ=255)                                                                                                                        tGGCTT‑CGCCATTAATAAGCCAGC‑TTTT‑AT‑CACCTTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGc                                >  1:63287/1‑151 (MQ=255)                                                                                                                          gtctcccccATAAATAAGCCAGT‑TTTT‑AC‑CCCTTTT‑ATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTg                                 <  2:11420/142‑1 (MQ=255)                                                                                                                          gcctcccccATAAATAAGCCAGT‑TTTT‑AT‑CCCTTTT‑ATCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtcttct                                  <  1:88170/141‑1 (MQ=255)                                                                                                                                 cctataaaaaagcAGG‑TTTTAAT‑CCCCTTT‑ATTCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCttttt                           <  1:646909/137‑1 (MQ=255)                                                                                                                                    ttAATAAGCCGGT‑TTTT‑AC‑CCCCTTT‑ATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTc                       <  2:813189/149‑1 (MQ=255)                                                                                                                                    ttAATAAGCCAGG‑TTTT‑AT‑CACCTTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCa                      >  2:561325/1‑150 (MQ=255)                                                                                                                                    ttAAAAAGCCAGT‑TTTT‑AC‑CCCTTTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCa                      <  2:571071/150‑1 (MQ=255)                                                                                                                                      aataaCCCAGG‑TTTT‑AC‑CCCTTTT‑TTTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaa                     <  2:297884/149‑1 (MQ=255)                                                                                                                                       aaaaacccggttttt‑AC‑CCCTTTT‑TTCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaaaa                   <  1:359927/140‑1 (MQ=255)                                                                                                                                          aGCCAGG‑ATTT‑AT‑CACCTGT‑ATTCCCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGCAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATa                 >  1:403010/1‑149 (MQ=255)                                                                                                                                            ccAGT‑TTTT‑AC‑CCCCTTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGa             <  1:800558/151‑1 (MQ=255)                                                                                                                                             ccggttttt‑AC‑CCCTTTT‑ATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGt               <  1:367793/144‑1 (MQ=255)                                                                                                                                             cAGG‑TTTT‑AT‑CACCTTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTg              >  1:924944/1‑149 (MQ=255)                                                                                                                                               gg‑TTTT‑AACCACTTTT‑ATCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGa             <  2:462325/149‑1 (MQ=255)                                                                                                                                                  tttt‑AT‑CCCTTTT‑ATCCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGc           <  2:909256/148‑1 (MQ=255)                                                                                                                                                          cACCTTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGAtttt  >  2:81238/1‑151 (MQ=255)                                                                                                                                                           |                                                                                                                                                       ATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCCAGG‑TTTT‑AT‑CACCTTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTT  >  NC_000913/4537960‑4538258

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