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.9%	115.7 / 20.8	72	G219G (GGT→GGG)	nanM	N‑acetylneuraminic acid mutarotase
Reads supporting (aligned to +/- strand):  ref base A (36/4);  new base C (0/30);  total (37/35)
Fisher's exact test for biased strand distribution p-value = 8.38e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGG‑TGAAATCAAGTTC‑AAATA‑CGGCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCC‑A‑GGTTTTATCACCTTTA‑TTCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATT  >  NC_000913/4537960‑4538256                                                                                                                                                             |                                                                                                                                                      aTCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGG‑TGAAATCAAGTTC‑AAATA‑CGGCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCC‑A‑GGTTTTATCa                                                                                                                                                       <  2:546737/150‑1 (MQ=255)        cAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGG‑TGAAATCAAGTTC‑AAATA‑CGGCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCGTTAATAAGCC‑A‑GGTTTTATCACCTTTA‑t                                                                                                                                               <  2:497245/150‑1 (MQ=255)        cAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGG‑TGAAATCAAGTTC‑AAATA‑CGGCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCC‑A‑GGTTTTATCACCTTTA‑tt                                                                                                                                              >  2:18882/1‑151 (MQ=255)          aaaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTCCCGG‑TGAAATCAAGTTC‑AAATA‑CGGCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCC‑A‑GGTTTTATCACCTTTA‑TTca                                                                                                                                            <  1:566890/151‑1 (MQ=255)            aCCGCCAGCTACGCCATCTGGAGATGAGACGGGAGCAAGCTTAATCCATTTTAAATTATTACCGG‑AGAAATCAAGTTC‑AAATA‑CGGCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCC‑A‑GGTTTTATCACCTTTA‑TTcac                                                                                                                                           >  2:759356/1‑150 (MQ=255)                 cAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGG‑TGAAATCAAGTTC‑AAATA‑CGGCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCC‑A‑GGTTTTATCACCTTTA‑TTCACAA‑ccgc                                                                                                                                    >  2:364503/1‑151 (MQ=255)                     tACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGG‑TGAAATCAAGTTC‑AAATA‑CGGCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCC‑A‑GGTTTTATCACCTTTA‑TTCACAA‑Ccgccgc                                                                                                                                 >  1:106664/1‑150 (MQ=255)                         ccATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGG‑TGAAATCAAGTTC‑AAATA‑CGGCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCC‑A‑GGTTTTATCACCTTTA‑TTCACAA‑CCGCCGCAc                                                                                                                               >  1:830904/1‑148 (MQ=255)                              tGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGG‑TGAAATCAAGTTC‑AAATA‑CGGCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCC‑A‑GGTTTTATCACCTTTA‑TTCACAA‑CCGCCGCACCAGCCGtt                                                                                                                       >  1:646194/1‑151 (MQ=255)                                      agaCGGGAGCAAGCTTATTCCATTTTAAATTATTACCGG‑TGAAATCAAGTTC‑AAATA‑CGGCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCC‑A‑GGTTTTATCACCTTTA‑TTCACAA‑CCGCCGCACCAGCCGTTCCGTAcc                                                                                                                >  2:528710/1‑150 (MQ=255)                                          gggAGCAAGTTTATCCCATTTAAAATTATTCCGGG‑TAAAATCAATTTA‑AAATC‑CGGCATCCTTCCGCAACCCTGT‑TTTGGCTT‑CGCCATTAATAAGCC‑A‑GTTTTTACCCCCTTTA‑TTCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGc                                                                                                           <  2:418214/151‑1 (MQ=255)                                          gggAGCAAGCTTATTCCATTTTAAATTATTACCGG‑TGAAATCAAGTTC‑AAATA‑CGGCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCC‑A‑GGTTTTATCACCTTTA‑TTCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGc                                                                                                           >  1:405605/1‑151 (MQ=255)                                            gAGCAAGCTTATTCCATTTTAAATTATTACCGG‑TGAAATCAAGTTC‑AAATA‑CGGCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCC‑A‑GGTTTTATCACCTTTA‑TTCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCg                                                                                                          >  2:126145/1‑150 (MQ=255)                                                aaGCTT‑TTCCATTTTAAATTTTTTCCGGTTAAAATCAATTTA‑AAATA‑GGCCATCCTTTCGCATCCCTGT‑TTTGGTTT‑CCCCTTAAATAAGCA‑A‑GTTTTTATCCCTTTTT‑TTCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTc                                                                                                      <  2:344403/150‑1 (MQ=255)                                                 aGCTTATTCCATTTTAAATTATTACCGG‑TGAAATCAAGTTC‑AAATA‑CGGCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCC‑A‑GGTTTTATCACCTTTA‑TTCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGc                                                                                                    >  2:184090/1‑151 (MQ=255)                                                  gttgcctCCTTTTTAATTTATTCCTGG‑TAAATTCAAG‑TA‑AAATACCGTCACCCCTTCACAACCCTGG‑TAT‑GTTT‑CCTCATTAAAAAACC‑AGGGTTTTACCCCTTTTA‑TCAAAAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGa                                                                                                         <  1:42081/139‑1 (MQ=255)                                                  gCTTATTCCATTTTAAATTATTACCGG‑TGAAATCAAGTTC‑AAATA‑CGGCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCC‑A‑GGTTTTATCACCTTTA‑TTCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc                                                                                                   >  1:376028/1‑151 (MQ=255)                                                         ccATTTTAAATTATTACCGG‑TGAAATCAAGTTC‑AAATA‑CGGCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCC‑A‑GGTTTTATCACCTTTA‑TTCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTa                                                                                             >  1:659410/1‑150 (MQ=255)                                                         catttttAAATTATTACCGG‑AAAAATCAAGTAC‑AAATC‑CGCCACCCTTTCACATCCCGGG‑TTTGGTTC‑CCCCTTTAATAAGCC‑G‑GTTTTTATC‑CCTTTATTTCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTaa                                                                                            <  2:675141/148‑1 (MQ=255)                                                          cATTTTAAATTATTACCGG‑TGAAATCAAGTTC‑AAATA‑CGGCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCC‑A‑GGTTTTATCACCTTTA‑TTCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAAc                                                                                           >  1:51704/1‑151 (MQ=255)                                                            ttttAATTTATT‑CCGG‑TTAAAATAAATTCAAAATA‑CGGCACCCTTCCGAAACCCTGT‑TTGGGCTC‑CGCAATTAATAGCCC‑A‑GTTTTTATCCCTTTTA‑TCCAAAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAAc                                                                                           <  2:173136/149‑1 (MQ=255)                                                            tttaAAATTTTTCCGGG‑AAAAATCAATTTA‑AAAAC‑CGCCCCCCGTCCGC‑ATCCTTG‑TTTGGGTTCCCCCATAAAAAAGCC‑G‑GTTTTTATCCCTTTTA‑TCCACAC‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTa                                                                                             <  1:671561/147‑1 (MQ=255)                                                            tttaAAATTATTCCGGG‑TAAAACCAATTTA‑AAATC‑CGGCACCCTTTCCCAATCCTGT‑TTGGGTTT‑CCCCATAAAAAAGCC‑G‑GGTTTTATCCCCTTTA‑TCCAAAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTa                                                                                             <  2:829966/147‑1 (MQ=255)                                                             tttAAATTATTACCGG‑TGAAATCAAGTTC‑AAATA‑CGGCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCC‑A‑GGTTTTATCACCTTTA‑TTCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTcc                                                                                        >  1:384859/1‑151 (MQ=255)                                                               tAAATTATTACCGG‑TGAAATCAAGTTC‑AAATA‑CGGCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCC‑A‑GGTTTTATCACCTTTA‑TTCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCAt                                                                                      >  2:848315/1‑151 (MQ=255)                                                               tAAATTATTACCGG‑TGAAATCAAGTTC‑AAATA‑CGGCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCC‑A‑GGTTTTATCACCTTTA‑TTCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCAt                                                                                      >  1:243506/1‑151 (MQ=255)                                                                  attattACCGG‑TGAAATCAAGTTC‑AAATA‑CGGCATCCGTTCGCAATCCTGG‑TTTGGCAT‑CGCCATTAATAAGCC‑A‑GGTTTTATCACCTTTA‑TTCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGc                                                                                   >  1:197375/1‑151 (MQ=255)                                                                        aCCGG‑TGAAATCAAGTTC‑AAATA‑CGGCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCC‑A‑GGTTTTATCACCTTTA‑TTCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGtt                                                                             >  1:361227/1‑151 (MQ=255)                                                                          cGG‑TGAAATCAAGTTC‑AAATA‑CGGCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCC‑A‑GGTTTTATCACCTTTA‑TTCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTg                                                                            >  2:182696/1‑150 (MQ=255)                                                                           gg‑TAAAATCAATTTA‑AAATA‑CGCCATCCTTCCGCAATCTTGT‑TTTGGCTT‑CCCCTTAAATAAGCC‑A‑GTTTTTATCCCCTTTA‑TCCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGa                                                                           <  2:348581/150‑1 (MQ=255)                                                                                 aaTCAATTTA‑AAATA‑CGGCATCCTTCCGCAATCCTGG‑TTTGGCTT‑CCCCATTAATAAGCC‑A‑GTTTTTACCCCCTTTA‑TCCAAAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATc                                                                     <  2:608362/151‑1 (MQ=255)                                                                                 aaTCAAGTTC‑AAATA‑CGGCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCC‑A‑GGTTTTATCACCTTTA‑TTCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATc                                                                     >  2:674928/1‑151 (MQ=255)                                                                                      aatttaaaaTA‑CGGCACCCGTCCGCAATCCTGT‑TTTGGTTC‑CCCCATTAATAAGCC‑A‑GTTTTTATCCCCTTTA‑TCCAAAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaa                                                                   <  2:715828/143‑1 (MQ=255)                                                                                             aaTC‑CGGCACCCTTTCGAAATCCTGT‑TTTGGCTT‑CCCCATTAATAAGCC‑A‑GGTTTTCTCCCCTTTA‑TCCCAAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaaca                                                          <  2:71876/151‑1 (MQ=255)                                                                                               tA‑CGGCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCC‑A‑GGTTTTATCACCTTTA‑TTCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaac                                                           >  1:585391/1‑148 (MQ=255)                                                                                                   ggCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCC‑A‑GGTTTTATCTCCTTTA‑TTCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAAc                                                     >  2:493657/1‑151 (MQ=255)                                                                                                   ggCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCC‑A‑GGTTTTATCACCTTTA‑TTCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAAc                                                     >  2:460918/1‑151 (MQ=255)                                                                                                   gccctccttccGCACTCCTTG‑TTTGGCTC‑CCCAATAAAAAACCC‑G‑GTTTTTACCCCTTTTA‑TCCACAC‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGaaa                                                      <  1:173195/140‑1 (MQ=255)                                                                                                       tCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCC‑A‑GGTTTTATCACCTTTA‑TTCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTAt                                                 >  1:752330/1‑151 (MQ=255)                                                                                                        ccGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCC‑A‑GGTTTTATCACCTTTA‑TTCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTAtt                                                >  2:381156/1‑151 (MQ=255)                                                                                                         ccgtccgaaTCCTGGTTTTGGCTG‑CCCCATTAAAAAGCG‑A‑GGTTTTTCCCCTTTTA‑TCCCCAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTAtt                                                <  1:416998/144‑1 (MQ=255)                                                                                                           ttCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCC‑A‑GGTTTTATCACCTTTA‑TTCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTga                                              >  2:702360/1‑150 (MQ=255)                                                                                                           ttCGCAATCCTGG‑TTGGGTTT‑CCCCATTAATAAGCC‑G‑GGTTTTATCCCCTTTA‑TCCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTa                                                  <  2:149754/146‑1 (MQ=255)                                                                                                           ttCGCAACCTTGT‑TTTGGCTC‑CCCCATAAATAAGCC‑A‑GTTTTTATCCCCTTTA‑TCCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTg                                               <  2:76285/149‑1 (MQ=255)                                                                                                            tCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCC‑A‑GGTTTTATCACCTTTA‑TTCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTga                                              >  2:22145/1‑149 (MQ=255)                                                                                                             cGCAATCCTGT‑TTTGGCTC‑CGCCTTTAAAAAGCC‑A‑GTTTTTATCCCTTTTA‑TTCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaaga                                           <  2:310202/151‑1 (MQ=255)                                                                                                             cGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCC‑A‑GGTTTTATCACCTTTA‑TTCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTga                                              >  2:761035/1‑148 (MQ=255)                                                                                                                     tttttttGGTTA‑CCCCATAAAAAACCA‑A‑GGTTTTACCCCCTTTA‑CTCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGaaa                                         <  1:44020/142‑1 (MQ=255)                                                                                                                      gg‑TTTGGCTT‑CGCCATTAATAAGCC‑A‑GGTTTTATCACCTTTA‑TTCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtctt                                  >  2:145213/1‑151 (MQ=255)                                                                                                                         tttgcttt‑CGCTATAAATATGCC‑A‑GTTTTAAACACATTTA‑TCAACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtct                                   <  1:759261/142‑1 (MQ=255)                                                                                                                         tttGGCTT‑CGCCATTAATAAGCC‑A‑GTTTTTATCCCCTTTA‑TCCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaataat                                     <  2:657784/146‑1 (MQ=255)                                                                                                                         tttGGCTT‑CGCAATTAA‑AAGCCTA‑GGTTTTATCCCTTTTT‑TCCAAAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaataat                                     <  2:690359/146‑1 (MQ=255)                                                                                                                          ttggcctcccccATAAAAAACCC‑A‑GTTTTTATCCCTTTTA‑TCCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtctt                                  <  1:218912/139‑1 (MQ=255)                                                                                                                          ttGGCTT‑CGCCATTAATAAGCC‑A‑GGTTTTATCACCTTTA‑TTCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTg                               >  1:735675/1‑151 (MQ=255)                                                                                                                           tGGCCT‑CCCCAATAAAAAGCC‑A‑GGTTTTTTCCCCTTTT‑TTCCCAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtcttct                                <  1:789971/150‑1 (MQ=255)                                                                                                                                  cGCAATTAAAAAGCC‑A‑GTTTTTATCCCTTTTA‑TTCACAC‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttt                          <  1:817432/149‑1 (MQ=255)                                                                                                                                   cccATAAATAACCC‑A‑GTTTTTCTCCCTTTTA‑CTCACAC‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtcttct                                <  1:380193/141‑1 (MQ=255)                                                                                                                                      aTTAATAAGCC‑A‑GGTTTTATCACCTTTA‑TTCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTc                     >  2:475141/1‑150 (MQ=255)                                                                                                                                      aTTAATAAGCC‑A‑GGTTTTATCACCTTTA‑TTCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTc                     >  2:299414/1‑150 (MQ=255)                                                                                                                                         aataaGCC‑A‑GTTTTTATCCCCTTTA‑TCCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaaa                  <  2:266216/150‑1 (MQ=255)                                                                                                                                         aataaGCC‑A‑GTTTTTACCCCTTTTA‑TCCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTg                       <  1:364455/145‑1 (MQ=255)                                                                                                                                         aataaGCA‑G‑GTTTTTACCCCTTTTA‑TCCCAAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaaa                  <  1:191713/150‑1 (MQ=255)                                                                                                                                         aaaaggCC‑A‑GTTTTTATCCCTTTTC‑TCCAAAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCa                    <  1:568939/143‑1 (MQ=255)                                                                                                                                         aaaaaccc‑A‑GTTTTTACCCCTTTTA‑TCCACAC‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaaa                  <  2:106685/144‑1 (MQ=255)                                                                                                                                         aaaaaGCA‑A‑TGTTTTACCCCTTTTT‑TCCCCAC‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaaaa                 <  1:848203/148‑1 (MQ=255)                                                                                                                                          ataaGCC‑A‑GGTTTTATCACCTTTA‑TTCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAAt                >  2:149480/1‑151 (MQ=255)                                                                                                                                           taaGCC‑A‑GGTTTTATCACCTTTA‑TTCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAAt                >  2:329364/1‑150 (MQ=255)                                                                                                                                                 cggggTTTTATCCCTTTTA‑TTCAAAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGt             <  2:51711/146‑1 (MQ=255)                                                                                                                                                  cgggTTTTATCCCCTTTA‑TTCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAg          <  2:405654/149‑1 (MQ=255)                                                                                                                                                  a‑GGTTTTATCACCTTTA‑TTCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGc         >  1:546344/1‑151 (MQ=255)                                                                                                                                                    gtttttATCACCTTTA‑TCCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCa        <  2:254155/149‑1 (MQ=255)                                                                                                                                                        ttATCACCTTTA‑TTCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTg     >  1:256968/1‑150 (MQ=255)                                                                                                                                                          aTCACCTTTA‑TTCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGAtt  >  2:669911/1‑151 (MQ=255)                                                                                                                                                             |                                                                                                                                                      ATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGG‑TGAAATCAAGTTC‑AAATA‑CGGCATCCGTTCGCAATCCTGG‑TTTGGCTT‑CGCCATTAATAAGCC‑A‑GGTTTTATCACCTTTA‑TTCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATT  >  NC_000913/4537960‑4538256

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