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.2%	114.1 / 22.0	71	G219G (GGT→GGG)	nanM	N‑acetylneuraminic acid mutarotase
Reads supporting (aligned to +/- strand):  ref base A (31/10);  new base C (0/30);  total (31/40)
Fisher's exact test for biased strand distribution p-value = 1.21e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTC‑AAA‑TACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTG  >  NC_000913/4537961‑4538253                                                                                                                                                         |                                                                                                                                                 tCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTC‑AAA‑TACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGGTTTTATCAc                                                                                                                                                 <  2:1016239/150‑1 (MQ=255)   cccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTC‑AAA‑TACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGGTTTTATCACCtt                                                                                                                                              >  1:377830/1‑151 (MQ=255)           aCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTC‑AAA‑TACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGGTTTTATCACCTTTATTcaca                                                                                                                                      <  2:736807/151‑1 (MQ=255)           aCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTC‑AAA‑TACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGGTTTTATCACCTTTATTcaca                                                                                                                                      <  2:305000/151‑1 (MQ=255)            ccGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTC‑AAA‑TACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGGTTTTATCACCTTTATTcac                                                                                                                                       <  2:935737/149‑1 (MQ=255)            ccGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTC‑AAA‑TACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGGTTTTATCACCTTTATTCACaa                                                                                                                                     <  2:820614/151‑1 (MQ=255)              gCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTC‑AAA‑TACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGGTTTTATCACCTTTATTCACAAcc                                                                                                                                   <  1:252688/151‑1 (MQ=255)                   cTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTC‑AAA‑TACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGGTTTTATCACCTTTATTCACAACcgccgc                                                                                                                              >  1:658111/1‑151 (MQ=255)                    tACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTC‑AAA‑TACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGGTTTTATCACCTTTATTCACAACcgccgc                                                                                                                              >  2:351676/1‑150 (MQ=255)                        ccttCTGGTGATGAGCCGGGAGCAAGCTTATTCCATTTTAAT‑TTATTACCGGTGAAATAAAGTTA‑AAA‑TACGGCACCCGTTCCCAATCCTGGTTGGCCTTCCCCATT‑AATAAGCCAGGTTTTACCCCCTTTATTCCCAAccgccg                                                                                                                               <  2:568106/142‑1 (MQ=255)                             tGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTC‑AAA‑TACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGtt                                                                                                                    >  1:596484/1‑151 (MQ=255)                               gTGATGAGGCGGGAGCAATTTTATTCCATTTAAAT‑TTATTCCCGGTAAAATAAAGTCA‑AAA‑TACGCCACCCTTCCGAAACCCTGTTTGGCCTCCCCCATA‑AAAAACCCGGTTTTTATCCCTTTTATCCACAACCGCCGCACCAGCCGTTcc                                                                                                                  <  2:794652/151‑1 (MQ=255)                                 gatgaGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTC‑AAA‑TACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGt                                                                                                                >  1:657751/1‑151 (MQ=255)                                  atgaGAGGGGCGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTC‑AAA‑TACGACATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGGTTTTATCACCTTTATTCACAACCGACGCACCAGCCGTTCCGt                                                                                                                >  2:755604/1‑150 (MQ=255)                                               aaGCTTATTCCATTTTAAA‑TTATTACCGGGGAAATCAAGTTC‑AAA‑TACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTc                                                                                                   >  1:1086183/1‑150 (MQ=255)                                                 gCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTC‑AAA‑TACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc                                                                                                >  1:427388/1‑151 (MQ=255)                                                 gCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTC‑AAA‑TACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc                                                                                                >  1:628204/1‑151 (MQ=255)                                                    tttcccttttaAAA‑TTATTCCGGGGGAAATCAAGTTA‑AAA‑CACGGCACCGGTTCGCAACCCTGTGTGGCCTTCGCCATT‑AATAAGCCAGGTTTAACCCCCTTTATTCACAACCCCCGCACCAGCCGTTCCGTACCAGGGCGATTCGc                                                                                                 <  1:1109901/140‑1 (MQ=255)                                                      ttCCTTTTTAAA‑TTTTTCCCGGTAAAATCAAGTTA‑AAA‑TCCGGCTCCCGTCCCCAATCCTGGTTTGGCTTCGCCTTT‑AAAAACCAAGTTTTTATCCCCTTTATTCCCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGc                                                                                             <  1:1103755/149‑1 (MQ=255)                                                       tCCCTTTTAAATTTATTCCCGGTAAAATCAATTAC‑AAA‑TACGGCATCCTTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCGGTTTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCg                                                                                            <  1:562801/150‑1 (MQ=255)                                                         cATTTTAAA‑TTATTACCGGTGAAATCAAGTTC‑AAA‑TACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAAc                                                                                        >  1:297331/1‑151 (MQ=255)                                                           ttttAAA‑TTATTCCCGGT‑AAAACAAATTCAAAA‑TACGCCTCCCTTCCGAAACCTTGTTTTGTCTCCCCCATA‑AAAAACCCAGTTTTTATCCCCTTTATCCCCACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTa                                                                                          <  2:455212/147‑1 (MQ=255)                                                           tttaAAA‑TTATTCCCGGTAAAATAAAGTTA‑AAA‑TACGGCATCCGTCCGCAACCCTGGTTTGGCTTCGCCATT‑AATAAGCCGGGTTTTATCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACt                                                                                       <  1:927176/150‑1 (MQ=255)                                                                 a‑ttattACCGGTGAAATCAAGTTC‑AAA‑TACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGc                                                                                >  2:211624/1‑151 (MQ=255)                                                                      ttccgggTAAAATCAATTTC‑AAA‑TACGGCTCCTTTTCCCACCCCTGTTTTGTTTCCCCAATA‑AAAAAGCCAGTTTTTACCCCTTTTTTCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCt                                                                               <  2:619462/143‑1 (MQ=255)                                                                      ttaCCGGTGAAATCAAGTTC‑AAA‑TACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtgtg                                                                            >  1:545252/1‑151 (MQ=255)                                                                      ttaCCGGTAAAATAAAGTTA‑AAA‑TGCGGCATCCGTTCGCACTCCGGTTTGGGCTTCGCCATA‑AATAAGCCAGGTTTTATCACCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtgt                                                                             <  1:521396/150‑1 (MQ=255)                                                                      ttaCCGGGAAAACCAAGTTA‑AAA‑TCCGCCTCCCTTTCGCACTCCTGTTTGGCCTTCGCCATT‑AAAAAGCCA‑GTTTTACCCCCTTTATTCCCACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtgt                                                                             <  2:619685/149‑1 (MQ=255)                                                                       ttcccggtaaaaCAAATTCAAAA‑TCCGCCATCCTTTCGCAACCCTGTTTTGCTTTCGCCTTT‑AATAACCCGGTTTTTACCCCTTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtg                                                                              <  1:485719/141‑1 (MQ=255)                                                                       taCCGGTGAAATCAAGTTC‑AAA‑TACGGCATCCGTTCGCAATCNNGNTTTGGCTTCGCCATT‑AATAAGCCAGGTTTTATCACNNTTANNnnc                                                                                                                                       >  2:741162/1‑91 (MQ=255)                                                                         cgggTAAAATAAATTTA‑AAA‑TACGGCCCCCGTTCGAATCCCTGTTTTGGCTTCGCCATT‑AATAACCCAGGTTTTACCCCCTTTATCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTgtgt                                                                           <  2:225117/147‑1 (MQ=255)                                                                         ccGGTGAAATCAAGTTC‑AAA‑TACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTg                                                                         >  1:1076132/1‑151 (MQ=255)                                                                               aaaTCAATTTA‑AAA‑TCCGGCACCCGTCCGCATCCGTGGTTGGCCTTCGCCATT‑AATAGCCCGGTTTTTACCCCTTTTATCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGAt                                                                   <  2:692931/151‑1 (MQ=255)                                                                               aaaTCAA‑TTC‑AAATTCCGCCTTCCGTTGGCAATCTTGTTTTGGCTTCCCCATA‑AAAAGGCCAGGTTTTATCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGAt                                                                   <  1:197322/151‑1 (MQ=255)                                                                                aaTCAAGTTC‑AAA‑TACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATc                                                                  >  1:271933/1‑151 (MQ=255)                                                                                 aTCAAGTTC‑AAA‑TACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATc                                                                  >  2:226231/1‑150 (MQ=255)                                                                                  tCAAGTTC‑AAA‑TACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCa                                                                 >  1:144829/1‑150 (MQ=255)                                                                                   cAAGTTC‑AAA‑TACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaa                                                                >  2:306656/1‑150 (MQ=255)                                                                                   cAAGTTC‑AAA‑TACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaa                                                                >  2:249020/1‑150 (MQ=255)                                                                                    atttac‑AAA‑TCCGCCATCGTTTCGCATTCTTGGTTTGTCTCCGCTAAT‑AAAAACCCGGTTTTTATCCCCTTTATTCACAGCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaaa                                                               <  1:675447/145‑1 (MQ=255)                                                                                    aggTTA‑AAA‑TCCGGCATCCGTTCGCAATCCTGTTTGGGCTTCGCCATT‑AATAAGCCAGTTTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCa                                                                 <  1:998911/146‑1 (MQ=255)                                                                                           aaa‑TACGGCACCCGTTCGAAACCCTGGTTGGGCTTCCCATTT‑AAAAAGCCAGTTTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaa                                                         <  1:28555/150‑1 (MQ=255)                                                                                             a‑TACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaac                                                        >  1:228960/1‑149 (MQ=255)                                                                                             a‑TACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaac                                                        >  1:927407/1‑149 (MQ=255)                                                                                               tACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaac                                                        >  2:991544/1‑148 (MQ=255)                                                                                                  ggCTTCCGTTCGAATTCCGGTTTTGGTTTCCCTTTT‑AATAACCCGTTTTTTTTCACTTTTTTCCCAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAg                                                      <  1:255185/147‑1 (MQ=255)                                                                                                  ggCATGCGTTGGCAATGCTGGGTTGGCTTCGCCATG‑AATAAGCCAGGTTTTATCACCTTTATTCACGACCGCCGCACAAGCCGTTGCGTATCAGGGCGATTCGCCAGCGTAACTCCATTGCTGGGTTGAGGGATCAAAAGACAACAGAAAc                                                  >  2:326777/1‑151 (MQ=255)                                                                                                    cATCCGTTCGCAACCTTGTTTTGGTTTCGCCATT‑AATAAGCCAGGTTTATTCCCCTTTATTCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAAc                                                  <  1:825510/149‑1 (MQ=255)                                                                                                       ccGTTCGCAACCCTGTTTTGGTTTCCCCATA‑AAAAAGCCAGTTTTTACCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTAtt                                             <  1:385964/151‑1 (MQ=255)                                                                                                            cgaaccccggttttGGCTTC‑CCATTAAATAAGCCGGTTTTTATCACCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaa                                          <  1:55563/138‑1 (MQ=255)                                                                                                               aacccTGGTTTGGTTCCCCCATT‑AATAAGCCAGGTTTTATCCCTTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGaa                                                    <  2:1053851/133‑1 (MQ=255)                                                                                                                     gtttgggCTCCGCAATT‑AATAAGCCAGTTTTTACCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtc                                 <  1:787433/144‑1 (MQ=255)                                                                                                                         tGGTTTGGCCATA‑AAAAGGCCAGGTTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtcttct                             <  1:781864/149‑1 (MQ=255)                                                                                                                           ctttCGCAATA‑AATAACCCGGTTTTTATCCCCTTTTTCCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTg                            <  2:491336/146‑1 (MQ=255)                                                                                                                                gCCATT‑AATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttttt                     >  1:377456/1‑150 (MQ=255)                                                                                                                                 ccATT‑AATAAGCCGGTTTTTATCCCTTTTATCCCCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGt                   <  1:789037/151‑1 (MQ=255)                                                                                                                                 ccATT‑AATAAGCCAGTTTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTg                    <  1:1135702/150‑1 (MQ=255)                                                                                                                                 ccATT‑AATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGt                   >  1:905365/1‑151 (MQ=255)                                                                                                                                 ccATT‑AATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGt                   >  2:733951/1‑151 (MQ=255)                                                                                                                                    tt‑AATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCa                 >  2:715499/1‑150 (MQ=255)                                                                                                                                       aataaGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaaa               >  1:163309/1‑150 (MQ=255)                                                                                                                                       aaaaacccAGTTTTTATCCCTTTTATCCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaaa               <  2:1054236/144‑1 (MQ=255)                                                                                                                                             ccggtttttACCACTTTTATCCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGt          <  2:467666/144‑1 (MQ=255)                                                                                                                                             ccAGTTTTTACCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGt          <  1:438884/149‑1 (MQ=255)                                                                                                                                             ccAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTg         >  1:1493/1‑150 (MQ=255)                                                                                                                                                 tttttATCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTg         <  2:657644/145‑1 (MQ=255)                                                                                                                                                 tttttACCCCTTTTATCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCa     <  2:1036866/149‑1 (MQ=255)                                                                                                                                                  ttttATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCAtt   >  2:1137862/1‑151 (MQ=255)                                                                                                                                                   tttACCCCCTTTATACACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCAtt   <  2:165232/150‑1 (MQ=255)                                                                                                                                                    ttATCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTg  <  1:229277/150‑1 (MQ=255)                                                                                                                                                    atATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTg  >  1:48095/2‑150 (MQ=255)                                                                                                                                                         |                                                                                                                                                 TCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTC‑AAA‑TACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTG  >  NC_000913/4537961‑4538253

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