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

CCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTT‑AAATTATTACCGGTGAAATCAAGTTC‑AAA‑TACGGCATCCGTTCGCAATCCTGG‑TTT‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTT  >  NC_000913/4537964‑4538258                                                                                                                                                       |                                                                                                                                                       ccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTT‑AAATTATTACCGGTGAAATCAAGTTC‑AAA‑TACGGCATCCGTTCGCAATCCTGG‑TTT‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCttt                                                                                                                                                   >  2:418509/1‑151 (MQ=255) ccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTT‑AAATTATTACCGGTGAAATCAAGTTC‑AAA‑TACGGCATCCGTTCGCAATCCTGG‑TTT‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCttt                                                                                                                                                   <  1:392161/151‑1 (MQ=255)   gCAAAACCGCCAGCTACGCCATCTGGTGATGAGCCGGGAGCAAGCTTATTCCATTTT‑AAATTATTCCCGGTGAAATCAAGTTC‑AAA‑TACCGCATCCGTTCGCAATCCTGG‑TTT‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTa                                                                                                                                                  <  1:283478/150‑1 (MQ=255)       aaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTT‑AAATTATTACCGGTGAAATCAAGTTC‑AAA‑TACGGCATCCGTTCGCAATCCTGG‑TTT‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATT‑cac                                                                                                                                            >  1:450093/1‑151 (MQ=255)             cAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTT‑AAATTATTACCGGTGAAATCAAGTTC‑AAA‑TACGGCATCCGTTCGCAATCCTGG‑TTT‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATT‑CACAAccgc                                                                                                                                      <  2:144450/151‑1 (MQ=255)                 tACGCCATCTGGCGATGAGACGGGAGCAAGCTTATTCCATTTT‑AAATTATTACCGGTGAAATCAAGTTC‑AAA‑TACGGCATCCGTTCGCAATCCTGG‑TTT‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATT‑CACAACcgccgc                                                                                                                                   >  1:690481/1‑150 (MQ=255)                  aCGCCATCGGGTGATGAGACGGGAGCAAGCTTATTCCTTTTA‑AAATTATTCCCGGTAAAATCAAGTTA‑AAA‑TCCGGCATCCTTTCGCAACCCTGG‑TTG‑GGCTTCGCCATTAATAAGCCGGGTTTTACCCCCTTAATT‑CACAACcgccgc                                                                                                                                   <  1:795061/149‑1 (MQ=255)                      cATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTT‑AAATTATTACCGGTGAAATCAAGTTC‑AAA‑TACGGCATCCGTTCGCAATCCTGG‑TTT‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATT‑CACAACCGCCGCACCAGc                                                                                                                             >  2:450414/1‑151 (MQ=255)                      cATCTGGTGAGGAGAGGGGAGCAAGTTTTTTCCTTTTT‑AAATTATTACCGGTAAAATCAAGTCA‑AAA‑TCGGCCATCCGTCCGCAATCCTGT‑TTT‑GGCTTCGCCATAAATAAGCCGGGTTTTATCACCTTTATT‑CACAACCGCCGCACCAGc                                                                                                                             <  1:850797/151‑1 (MQ=255)                         cTGGTGATGAGACGGGAGCAAGCTTATTCCATTTT‑AAATTATTACCGGTGAAATCAAGTTC‑AAA‑TACGGCATCCGTTCGCAATCCTGG‑TTT‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGt                                                                                                                          >  1:77672/1‑151 (MQ=255)                                   gaCGGGAGCAAGCTTATTCCATTTT‑AAATTATTACCGGTGAAATCAAGTTC‑AAA‑TACGGCATCCGTTCGCAATCCTGG‑TTT‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTAcc                                                                                                                  >  1:672866/1‑149 (MQ=255)                                      gggAGCAATCTTTTTCCATTTA‑AAATTATTCCGGGTAAAATCAATTTA‑AAA‑TACGGCCTCCGTCCCCATCCCGGG‑TTG‑GCCTCCGCCATAAAAAACCCAGTTTTTACCCCTTTTATA‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGc                                                                                                             <  1:184581/151‑1 (MQ=255)                                      gggAGCAAGTTTTTTCCTTTTT‑AAATTATTCCCGGTAAAATCAAGTTC‑AAA‑TACGCCATCCTTCCGCAATCCTGT‑TTT‑GGTTTCCCCATAAATAAGCCAGTTTTTATCCCTTTTATT‑CAAAACCGCCGCACCAGCCGTTCCGTACCAg                                                                                                                <  2:153225/148‑1 (MQ=255)                                      gggAGCAACCTT‑TTCCTTTTTAAAATTATT‑CCGGGAAAAACAAATTCAAAA‑TCCGGCACCCTTTCAC‑ATCCTGG‑TTTGGGCTTCCCCATTAAAACGCCAGTTTTTACCCCTTTTATT‑CAAACCCGCCGCACCAGCCGTTCCGTACCAGGGc                                                                                                             <  1:460567/151‑1 (MQ=255)                                       ggAGCAAGCTTATTCCATTTT‑AAATTATTACCGGTGAAATCAAGTTC‑AAA‑TACGGCATCCGTTCGCAATCCTGG‑TTT‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCg                                                                                                            >  1:1034495/1‑151 (MQ=255)                                            aaGCTTATTCCATTTT‑AAATTATTACCGGTGAAATCAAGTTC‑AAA‑TACGGCATCCGTTCGCAATCCTGG‑TTT‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCg                                                                                                       >  1:442717/1‑151 (MQ=255)                                            aaGCTTATTCCATTTT‑AAATTATTACCGGTGAAATCAAGTTC‑AAA‑TACGGCATCCGTTCGCAATCCTGG‑TTT‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCg                                                                                                       >  2:473945/1‑151 (MQ=255)                                            aaGCTTATTCCATTTT‑AAATTATTACCGGTGAAATCAAGTTC‑AAA‑TACGGCATCCGTTCGCAATCCTGG‑TTT‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCg                                                                                                       >  1:825965/1‑151 (MQ=255)                                             aGCTTATTCCATTTT‑AAATTATTACCGGTGAAATCAAGTTC‑AAA‑TACGGCATCCGTTCGCAATCCTGG‑TTT‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGc                                                                                                      >  2:340503/1‑151 (MQ=255)                                              gCTTATTCCATTTT‑AAATTATTACCGGTGAAATCAAGTTC‑AAA‑TACGGCATCCGTTCGCAATCCTGG‑TTT‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc                                                                                                     >  2:1055763/1‑151 (MQ=255)                                              agtttttCCAATTT‑AAATTTTTACCGGTGAAAATAAATTA‑AAA‑AGCCGCCCCCCTTCCGAAACCTGG‑TTTGGGCTTCCCCTTAAATAA‑CCAGGTTTTAACCCCTTTTTTCCAAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGc                                                                                                      <  1:824862/146‑1 (MQ=255)                                                         ttt‑AAATTATTACCGGTGAAATCAAGTTC‑AAA‑TACGGCATCCGTTCGCAATCCTGG‑TTT‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTcc                                                                                          >  1:285024/1‑151 (MQ=255)                                                          tt‑AAATTATTACCGGTGAAATCAAGTTC‑AAA‑TACGGCATCCGTTCGCAATCCTGG‑TTT‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTcc                                                                                          >  2:242383/1‑150 (MQ=255)                                                           taaaaTTATTCCCGGTAAAAACAAGTAA‑AAA‑TAGGGCCCCCGTCCCAATTCCTGT‑TTT‑GTCTCCGCAATAAAAACCCCGGTTTTTACCCCTTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTcc                                                                                          <  2:606314/148‑1 (MQ=255)                                                             aaaTTATTACCGGTGAAATCAAGTTC‑AAA‑TACGGCATCCGTTCGCAATCCTGG‑TTT‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCAtt                                                                                       >  2:211806/1‑151 (MQ=255)                                                                ttattaCCGGTGAAATCAAGTTC‑AAA‑TACGGCATCCGTTCGCAATCCTGG‑TTT‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCt                                                                                    >  1:923652/1‑151 (MQ=255)                                                                  attaCCGGTGAAATCAAGTTC‑AAA‑TACGGCATCCGTTCGCAATCCTGG‑TTT‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtgt                                                                                  >  1:1029614/1‑151 (MQ=255)                                                                   ttcccGGTAAAATCAATTTC‑AAA‑TACGGCATCCGTTCGCAATCCTGT‑TTT‑GCCTTCGCCATTAATAAGCCAGGTTTTACCCCCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtgt                                                                                  <  2:813790/147‑1 (MQ=255)                                                                      cgggTAAAATCAAGTAC‑AAA‑TACGGCATCCGTTCGCAATCCTGG‑TTT‑GGCTTCCCCATTAATAAGCCGGTTTTTATCCCTTTTATC‑CAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtgtg                                                                                 <  2:676068/146‑1 (MQ=255)                                                                      ccGGTAAAATCAATTTA‑AAA‑TCCGGCACCCGTTCGCAACCCTGT‑TTT‑GGTTTCGCCATTAAAAAGCCAGTTTTTATCCCCTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTg                                                                              <  2:843130/151‑1 (MQ=255)                                                                         ggtaaaacaaatttaaaa‑TCCGGCCCCCGTTCGCAATCCTGG‑TTG‑GGTTTCGCCTTTAATAAGCCAGTTTTTATCCCCTTTATC‑CACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGa                                                                             <  1:279757/135‑1 (MQ=255)                                                                         ggtaaaaCAAATTC‑AAATTCCGCCATCCTTTCGCAATCCTGT‑TTT‑GTTTTCCCCATAAAAACGCGGGTTTTTATCCCCTTTTTT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGa                                                                             <  2:349200/147‑1 (MQ=255)                                                                          tGAAATCAAGTTC‑AAA‑TACGGCATCCGTTCGCAATCCTGG‑TTT‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAggg                                                                          >  1:573692/1‑151 (MQ=255)                                                                             attCAATTTA‑AAA‑TCCGGCATCCTTTCGCAATCCTGG‑TTT‑GGCTCCCCCATTAATAACCCAGTTTTTACCCCCTTTATC‑CACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATc                                                                       <  1:750738/149‑1 (MQ=255)                                                                                cAAGTTC‑AAA‑TACGGCATCCGTTCGCAATCCTGG‑TTT‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaa                                                                     >  1:387423/1‑150 (MQ=255)                                                                                cAAGTTC‑AAA‑TACGGCATCCGTTCGCAATCCTGG‑TTT‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaa                                                                     >  1:445081/1‑150 (MQ=255)                                                                                cAAGTTC‑AAA‑TACGGCATCCGTTCGCAATCCTGG‑TTT‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaa                                                                     >  1:39890/1‑150 (MQ=255)                                                                                 aaGTTC‑AAA‑TACGGCATCCGTTCGCAATCCTGG‑TTT‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaaa                                                                    >  2:831900/1‑150 (MQ=255)                                                                                 aaGTTC‑AAA‑TACGGCATCCGTTCGCAATCCTGG‑TTT‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaaa                                                                    >  2:384999/1‑150 (MQ=255)                                                                                  aatttaaaa‑TCCGGCTTCCGTTCGCACTCCTGT‑TTT‑GGCTCCGCCATTAAAAAGCCAGTTTTTATCCCTTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCa                                                                      <  1:450326/142‑1 (MQ=255)                                                                                  aatttaaaa‑TCCGCCATCCTTTCGAAATCCTGT‑TTT‑GGTTTCCCTATTAATAGGCCAGTTTTTACCCCTTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCa                                                                      <  2:688442/142‑1 (MQ=255)                                                                                    ttC‑AAA‑TACGGCATCGTTTCGAAATCCGGT‑TTG‑GGCTTCCCCATTAATAGGCCAGTTTTTACCCCTTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGa                                                                 <  2:61861/150‑1 (MQ=255)                                                                                        aaa‑TCCGCCTTCCTTTCCCAACCCTGT‑TTT‑GGCTTCCCCTTTAATAACCCAGTTTTTATCCCCTTTTTC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGac                                                                <  2:149117/148‑1 (MQ=255)                                                                                        aaa‑TACGGCATCCGTTCGCAATCCTGG‑TTT‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaac                                                             >  1:1763/1‑151 (MQ=255)                                                                                          a‑TACGCCATCCTTCCGCAACCCTGT‑TTT‑GGTTCCGCTTTTAAAAACCAAGTTTTTACCACTTTTATT‑CAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaaca                                                            <  2:144180/150‑1 (MQ=255)                                                                                               gcctccttttCCCAATCCTGG‑TTG‑GGCTCCCCCTTTAAAAAACCCGGTTTAATCCCTTTTATC‑CCAAGCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGa                                                          <  1:297411/140‑1 (MQ=255)                                                                                               gccATCCTTTCCCATTCCTGT‑TTT‑GGCTTCCCAATAAATAGGCCGGTTTTTATCCCTTTTATC‑AAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGaaa                                                        <  1:142687/148‑1 (MQ=255)                                                                                                  acccGTTCCCATTCCTGT‑TTT‑GCCTTCGCATTAAATAAGCCGGGTTTTACCCCCTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACt                                                      <  2:118953/147‑1 (MQ=255)                                                                                                  aTCCGTTCGCAATCCTGG‑TTT‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACtt                                                     >  2:811126/1‑150 (MQ=255)                                                                                                    cctttCGCAATCCTGG‑TTT‑GGCTTCCCCATAAAAAAGCCGGTTTTTACCCCTTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACtt                                                     <  2:672981/145‑1 (MQ=255)                                                                                                    ccGTCCGCAATCCTGG‑TTT‑GGTTCCCCCATTAAAAACCCGGGTTTTACCCCCTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTa                                                    <  2:366740/149‑1 (MQ=255)                                                                                                       ttCGCAATCCTGG‑TTG‑GGCTTCCCCATTAAAAAGCCAGTTTTTACCCCTTTTATC‑CAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTAtt                                                  <  2:20761/148‑1 (MQ=255)                                                                                                       ttCGCAA‑CCTGGTTTT‑GGTTTCGCCATAAATAAGCCAGGTTTTATCCCCTTTATT‑CACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTga                                                <  2:982355/150‑1 (MQ=255)                                                                                                       tcccaaaTCCTGT‑TTT‑GGCTCCCCCATAAATAAGCCAGGTTTTATCACTTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAg                                                           <  1:328433/134‑1 (MQ=255)                                                                                                        tCGCAATCCTGG‑TTT‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTga                                                >  1:287453/1‑149 (MQ=255)                                                                                                           cAATCCTGG‑TTT‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaagaa                                            >  2:710806/1‑150 (MQ=255)                                                                                                                 tGG‑TTT‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCAATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtct                                     >  1:853119/1‑151 (MQ=255)                                                                                                                  gg‑TTT‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtctt                                    >  1:3136/1‑151 (MQ=255)                                                                                                                     ttt‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtcttct                                  >  2:712490/1‑151 (MQ=255)                                                                                                                      tt‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTg                                 >  2:66753/1‑151 (MQ=255)                                                                                                                           cTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCttt                             >  2:1056409/1‑151 (MQ=255)                                                                                                                            ttCGCCATAAAAAACCCGGTTTTTATCCCCTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGc                                <  2:1024751/147‑1 (MQ=255)                                                                                                                              cGCCATTAATAAGCCGGGTTTTATCCCCTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCttttt                           <  2:1012575/150‑1 (MQ=255)                                                                                                                                   ttAAAAACCCAGTTTTTATCCCTTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTc                       <  2:880906/149‑1 (MQ=255)                                                                                                                                     aaaaacccggtttttACCCCTTTTATC‑CACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGt                        <  1:897822/135‑1 (MQ=255)                                                                                                                                      ataaGCCGGTTTTTACCCCTTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAAt                  <  2:754722/151‑1 (MQ=255)                                                                                                                                           ccAGTTTTTACCCCCTTTATC‑CACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGt               <  1:1055588/149‑1 (MQ=255)                                                                                                                                             aGGTTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGc           >  1:801703/1‑151 (MQ=255)                                                                                                                                              gttttaATCCCTTTAATC‑CAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCa          <  1:537917/149‑1 (MQ=255)                                                                                                                                                tttaATCCCCTTTTTC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGa             <  1:1052729/146‑1 (MQ=255)                                                                                                                                                      cACTTTTATC‑CACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGAttt   <  2:578134/150‑1 (MQ=255)                                                                                                                                                      cACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGAtttt  >  2:909162/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                       CCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTT‑AAATTATTACCGGTGAAATCAAGTTC‑AAA‑TACGGCATCCGTTCGCAATCCTGG‑TTT‑GGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTT  >  NC_000913/4537964‑4538258

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