Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,538,109	0	A	C	39.4%	121.4 / 14.1	69	G219G (GGT→GGG)	nanM	N‑acetylneuraminic acid mutarotase
Reads supporting (aligned to +/- strand):  ref base A (39/1);  new base C (0/26);  total (39/30)
Fisher's exact test for biased strand distribution p-value = 1.63e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAA‑TACGGCA‑TCCGTTCGCAATCCTGGTTTGGCTT‑CGCC‑A‑TT‑AATAAGCCAGGTTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTT  >  NC_000913/4537960‑4538258                                                                                                                                                             |                                                                                                                                                        aTCCCCGCAAAAACGCCAGCTACGCCATCTGGTGATGAGACGGGCGCAAGATTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAA‑TACGGCA‑TCCGTTCGCAATCCTGGTTTGGCTT‑CGCC‑A‑TT‑AATAAGCCAGGTTTTATCa                                                                                                                                                         >  2:829209/1‑150 (MQ=255)  tCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAA‑TACGGCA‑TCCGTTCGCAATCCTGGTTTGGCTT‑CGCC‑A‑TT‑AATAAGCCAGGTTTTATCAcc                                                                                                                                                       >  1:505223/1‑151 (MQ=255)    cccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAA‑TACGGCA‑TCCGTTCGCAATCCTGGTTTGGCTT‑CGCC‑A‑TT‑AATAAGCCAGGTTTTATCACCtt                                                                                                                                                     >  2:992654/1‑151 (MQ=255)            aCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAA‑TACGGCA‑TCCGTTCGCAATCCTGGTTTGGCTT‑CGCC‑A‑TT‑AATAAGCCAGGTTTTATCACCTTTATT‑cac                                                                                                                                             >  1:441882/1‑150 (MQ=255)              cGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAA‑TACGGCA‑TCCGTTCGCAATCCTGGTTTGGCTT‑CGCC‑A‑TT‑AATAAGCCAGGTTTTATCACCTTTATT‑CACAA‑c                                                                                                                                         >  1:1000115/1‑151 (MQ=255)               gCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAA‑TACGGCA‑TCCGTTCGCAATCCTGGTTTGGCTT‑CGCC‑A‑TT‑AATAAGCCAGGTTTTATCACCTTTATT‑CACAA‑cc                                                                                                                                        <  1:732884/151‑1 (MQ=255)               gCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAA‑TACGGCA‑TCCGTTCGCAATCCTGGTTTGGCTT‑CGCC‑A‑TT‑AATAAGCCAGGTTTTATCACCTTTATT‑CACAA‑cc                                                                                                                                        >  2:242999/1‑151 (MQ=255)               gCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAA‑TACGGCA‑TCCGTTCGCAATCCTGGTTTGGCTT‑CGCC‑A‑TT‑AATAAGCCAGGTTTTATCACCTTTATT‑CACAA‑cc                                                                                                                                        >  1:970158/1‑151 (MQ=255)                  aGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAA‑TACGGCA‑TCCGTTCGCAATCCTGGTTTGGCTT‑CGCC‑A‑TT‑AATAAGCCAGGTTTTATCACCTTTATT‑CACAA‑ccgcc                                                                                                                                     >  1:220934/1‑151 (MQ=255)                   gCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAA‑TACGGCA‑TCCGTTCGCAATCCTGGTTTGGCTT‑CGCC‑A‑TT‑AATAAGCCAGGTTTTATCACCTTTATT‑CACAA‑ccgccg                                                                                                                                    >  2:201849/1‑151 (MQ=255)                   gCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAA‑TACGGCA‑TCCGTTCGCAATCCTGGTTTGGCTT‑CGCC‑A‑TT‑AATAAGCCAGGTTTTATCACCTTTATT‑CACAA‑ccgccg                                                                                                                                    >  1:559441/1‑151 (MQ=255)                       ctccATGTGGGGATGAAGCGGAAGCAATCTTTTTCCATTTTAAA‑TTATTACCGGTGAAATCAATTTAAAT‑TCCGGCA‑TCCTTCCCCAAGCTTGGTTTGGTTC‑CGCC‑A‑TA‑AATACGCCAGTTTTTACCCCTTTAATC‑CAAAA‑CCGCCGCa                                                                                                                                  <  2:792543/147‑1 (MQ=255)                       cGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAA‑TACGGCA‑TCCGTTCGCAATCCTGGTTTGGCTT‑CGCC‑A‑TT‑AATAAGCCAGGTTTTATCACCTTTATT‑CACAA‑CCGCCGCAcc                                                                                                                                >  2:385274/1‑151 (MQ=255)                             cTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAA‑TACGGCA‑TCCGTTCGCAATCCTGGTTTGGCTT‑CGCC‑A‑TT‑AATAAGCCAGGTTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGt                                                                                                                          >  1:593887/1‑151 (MQ=255)                               ggTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAA‑TACGGCA‑TCCGTTCGCAATCCTGGTTTGGCTT‑CGCC‑A‑TT‑AATAAGCCAGGTTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTc                                                                                                                        >  1:387523/1‑151 (MQ=255)                                 tgatgaGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAA‑TACGGCA‑TCCGTTCGCAATCCTGGTTTGGCTT‑CGCC‑A‑TT‑AATAAGCCAGGTTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCg                                                                                                                      >  2:813603/1‑151 (MQ=255)                                   aggagaCGGGAGAAAGTTTTTCCTTTTTTAAT‑TTATTCCGGGAAAAATAAATTTCAAA‑ACCGCCA‑CCCTTTCCCAACCCGGGTTTGCCTG‑CCCC‑T‑TT‑AAAAAGCCGGTGTTTACCCCCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCg                                                                                                                      <  2:834530/147‑1 (MQ=255)                                      agcggggAGCAAGCTTATTCTATTTAAAT‑TTATTACCGGTAAAATCAAGTCAAAA‑TACGGCT‑TCGTTTCGCACCCCTTTTTTGGCTC‑CCCC‑A‑TT‑AATAAGCCGGTTTATCCCCCCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTAc                                                                                                                   <  2:238417/145‑1 (MQ=255)                                         cGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAA‑TACGGCA‑TCCGTTCGCAATCCTGGTTTGGCTT‑CGCC‑A‑TT‑AATAAGCCAGGTTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAggg                                                                                                              >  2:930970/1‑151 (MQ=255)                                            gAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAA‑TACGGCA‑TCCGTTCGCAATCCTGGTTTGGCTT‑CGCC‑A‑TT‑AATAAGCCAGGTTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCg                                                                                                            >  2:185717/1‑150 (MQ=255)                                              gCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAA‑TACGGCA‑TCCGTTCGCAATCCTGGTTTGGCTT‑CGCC‑A‑TT‑AATAAGCCAGGTTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGAtt                                                                                                         >  2:508525/1‑151 (MQ=255)                                                aaGTTTTTTCCTTTTAAAA‑TTTTTCCGGGTAAAATCAAGTTAAAA‑TACGGCA‑TCCTTTCGCATTCCGGTTTTGGCTT‑CGCC‑A‑TT‑AAAAGGCCAGTTTTTACCCCTTTTATC‑CAAAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGAtt                                                                                                         <  1:195405/149‑1 (MQ=255)                                                aaGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAA‑TACGGCA‑TCCGTTCGCAATCCTGGTTTGGCTT‑CGCC‑A‑TT‑AATAAGCCAGGTTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCg                                                                                                       >  2:450907/1‑151 (MQ=255)                                                     tttttctttttAAATTTATTCCCGGT‑AAAACAAATTCAAACTCCGGCC‑CCCTTTCACATCCCTGGTTTGCTTC‑CCCC‑T‑TT‑AAAAAGCCAGTTTTTCTCCCCTTTACC‑CCAAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCa                                                                                                    <  2:1025931/143‑1 (MQ=255)                                                      aTTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAA‑TACGGCA‑TCCGTTCGCAATCCTGGTTTGGCTT‑CGCC‑A‑TT‑AATAAGCCAGGTTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCg                                                                                                 >  1:992834/1‑151 (MQ=255)                                                            ttttAAA‑TTATTCCCGGTAAAATCAATTCAAAA‑TCCGGCA‑TCCTTTCGCATCCCTGTTTTGGTTT‑CCCC‑A‑TT‑AAAAAGCCGGTTTTTACCCCCTTTATC‑CACAC‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAAc                                                                                             <  1:612015/149‑1 (MQ=255)                                                             tttAAA‑TTATTACCGGTGAAATCAAGTTCAAA‑TACGGCA‑TCCGTTCGCAATCCTGGTTTGGCTT‑CGCC‑A‑TT‑AATAAGCCAGGTTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTcc                                                                                          >  1:830993/1‑151 (MQ=255)                                                                aaa‑TTATTACCGGTGAAATCAAGTTCAAA‑TACGGCA‑TCCGTTCGCAATCCTGGTTTGGCTT‑CGCC‑A‑TT‑AATAAGCCAGGTTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCAtt                                                                                       >  1:828751/1‑151 (MQ=255)                                                                 aa‑TTATTACCGGTGAAATCAAGTTCAAA‑TACGGCA‑TCCGTTCGCAATCCTGGTTTGGCTT‑CGCC‑A‑TT‑AATAAGCCAGGTTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTg                                                                                      >  2:27004/1‑151 (MQ=255)                                                                     tattaCCGGTGAAATCAAGTTCAAA‑TACGGCA‑TCCGTTCGCAATCCTGGTTTGGCTT‑CGCC‑A‑TT‑AATAAGCCAGGTTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtg                                                                                   >  2:678918/1‑151 (MQ=255)                                                                      attaCCGGTGAAATCAAGTTCAAA‑TACGGCA‑TCCGTTCGCAATCCTGGTTTGGCTT‑CGCC‑A‑TT‑AATAAGCCAGGTTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtgt                                                                                  >  1:717489/1‑151 (MQ=255)                                                                      attaCCGGTGAAATCAAGTTCAAA‑TACGGCA‑TCCGTTCGCAATCCTGGTTTGGCTT‑CGCC‑A‑TT‑AATAAGCCAGGTTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtg                                                                                   >  2:7861/1‑150 (MQ=255)                                                                          ccGGTGAAATCAAGTTCAAA‑TACGGCA‑TCCGTTCGCAATCCTGGTTTGGCTT‑CGCC‑A‑TT‑AATAAGCCAGGTTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTg                                                                              >  2:679217/1‑151 (MQ=255)                                                                           cGGTTAAAACAAATTCCAA‑TCAGGCCTTCCTTTCGCACCCGTGGTTTGGCTT‑CCCC‑A‑TT‑AAAAACCCAGGTTTCACCGCATTTATT‑CCAAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGtt                                                                               <  1:682352/150‑1 (MQ=255)                                                                           cGGTGTAATCAAGCTGAAA‑TACGGCA‑TCCGATCGCAATCCTGGTTTGGCTT‑CGCC‑A‑TT‑AATTAGCCAGGCTTTATCACCTTTATT‑CACAA‑CCGCCGCATCAGCCCTTCCGTACCTGGCCGCTTCGCCAGCGTAACTCCATTGCTgtgt                                                                                >  2:635060/1‑148 (MQ=255)                                                                            ggggAAATAAAGTTAAAA‑TGGGGCA‑CCCGTCCGCAACCCGGTTTTGTCTT‑CCCC‑A‑TA‑AATACCCCGGGTTTTATCCCCTTTATA‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAg                                                                            <  1:715016/148‑1 (MQ=255)                                                                                aaaTCAAGTTCAAA‑TACGGCA‑TCCGTTCGCAATCCTGGTTTGGCTT‑CGCC‑A‑TT‑AATAAGCCAGGTTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGAt                                                                        >  2:741813/1‑151 (MQ=255)                                                                                     aattcaaaa‑TCCGGCA‑‑CCCTTCCCAAACCCTGTTTTGCTT‑CCCCAA‑TA‑AATAAGCCGGTTTTTATCCCCTTTTTC‑CAAAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCa                                                                      <  1:545413/142‑1 (MQ=255)                                                                                          cAAA‑TACGGCA‑TCCGTTCGCAATCCTGGTTTGGCTT‑CGCC‑A‑TT‑AATAAGCCAGGTTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGaca                                                               >  2:590855/1‑150 (MQ=255)                                                                                          cAAA‑TACGGCA‑TCCGTTCGCAATCCTGGTTTGGCTT‑CGCC‑A‑TT‑AATAAGCCAGGTTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGac                                                                >  2:442251/1‑149 (MQ=255)                                                                                             aaatccgcctccctttCGAAACCCGGTTTTGGTTT‑C‑CC‑A‑TTAAAAAAGCCAGTTTTTATCCCTTTTATC‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGaca                                                               <  1:547680/135‑1 (MQ=255)                                                                                               aaCGGCG‑TCTTTTCCCACTCGTGTTTTGGCTT‑CACC‑A‑TT‑AA‑AAGCCATTTTTTTAACCCTTTCTTCCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGa                                                          <  1:473210/150‑1 (MQ=255)                                                                                                       tccgtccgCAA‑CCTTTTTTGGGTTCCCCC‑A‑TT‑AATAACCCAGTTTTTACCCCTTTTATC‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTa                                                    <  2:791589/150‑1 (MQ=255)                                                                                                        ccttccacatcccTGTTTTGGCTT‑CCCC‑A‑AA‑AAAAACCCAGGTTTTTTCCCTTTTATC‑CACAC‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACtt                                                     <  1:92426/137‑1 (MQ=255)                                                                                                               cAATCCTGGTTTGGCTT‑CGCC‑A‑TT‑AATAAGCCAGGTTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaagaa                                            >  2:513576/1‑150 (MQ=255)                                                                                                                ccatcttggttGGCCT‑CGCC‑T‑TTAAAAAAGCAAGGTTTTA‑CCCCTTTTTT‑ACCAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaagaa                                            <  2:530747/141‑1 (MQ=255)                                                                                                                aaTCCTGGTTTGGCTT‑TGCC‑A‑TT‑AATAAGCCAGGTTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaaga                                             >  2:15631/1‑148 (MQ=255)                                                                                                                 aTCCTGGTTGGCTTT‑TGCC‑A‑TA‑AATAACCCAGGTTTTACCCCTTTTATC‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaaga                                             <  1:15632/147‑1 (MQ=255)                                                                                                                  tCCTGGTTTGGCTT‑CGCC‑A‑TT‑AATAAGCCAGGTTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaata                                         >  2:603483/1‑150 (MQ=255)                                                                                                                   ccTGTTTTGGCTC‑CCCA‑A‑TA‑AAAAAGCCAGTTTTTATCCCCTTTATC‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaata                                         <  1:123313/149‑1 (MQ=255)                                                                                                                    cctgTTTTGCTT‑CCCC‑CTTT‑AAAAGGCCTGTTTTTACCCCTTTTTTC‑CACAC‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACtg                                                     <  2:765606/134‑2 (MQ=255)                                                                                                                        tttGGCTT‑CGCT‑A‑TT‑AATAAGCCAGTTTTTATCCCTTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtct                                     <  1:992990/148‑1 (MQ=255)                                                                                                                        ttgggttt‑CGCC‑A‑TT‑AAAAAGCCAGTTTTTACCCCTTTTATC‑CACAC‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtcttc                                   <  1:10505/144‑1 (MQ=255)                                                                                                                            gcctcctcc‑A‑TT‑AATAAGCCGGTTTTAAACCCTTTTATC‑CAAAC‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGc                                <  2:487831/143‑1 (MQ=255)                                                                                                                            gCTT‑CGCC‑A‑TT‑AATAAGCCAGGTTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtt                              >  2:226676/1‑151 (MQ=255)                                                                                                                               ggctgcgaagtaaataaGCCTGGTTTTATCCCTTTTTTT‑CAAAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTg                                 <  1:748623/137‑1 (MQ=255)                                                                                                                                 cGCC‑A‑TT‑AATAAGCCAGGTTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACACCAGAAACTTATTGAAGAAATAATCTTCTGCtttttt                          >  2:543394/1‑151 (MQ=255)                                                                                                                                 cGCC‑A‑TT‑AAAAAGCCGGTTTTTACCCCTTTTATC‑CCCAC‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTg                                 <  2:177755/144‑1 (MQ=255)                                                                                                                                     ccttt‑AATAACCGGGGTTTTACCCCCTTTATT‑CTCAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCttttt                           <  1:660572/145‑1 (MQ=255)                                                                                                                                      a‑TT‑AATAAGCCAGGTTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTc                       >  1:41883/1‑150 (MQ=255)                                                                                                                                             taaGCCAGGTTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACAAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGCGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAAt                  >  2:671580/1‑150 (MQ=255)                                                                                                                                                cccgttttttACCGCTTTTATC‑CGCAC‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGAGTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAAt                  <  1:5749/141‑1 (MQ=255)                                                                                                                                                 ccAGTTTTTATCCCTTTTTTT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGt               <  1:353455/149‑1 (MQ=255)                                                                                                                                                  cAGGTTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTg              >  2:745639/1‑149 (MQ=255)                                                                                                                                                     tttttATCCCTTTTATC‑CACAC‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAg            <  1:414363/147‑1 (MQ=255)                                                                                                                                                     gTTTTATCACCTTTATT‑CACAA‑CCGCCGCACAAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGGTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCAt         >  2:706915/1‑151 (MQ=255)                                                                                                                                                      ttttATCTCCTTTATC‑CAAAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCAt         <  2:25289/150‑1 (MQ=255)                                                                                                                                                      ttttACCCCCTTTATT‑CACAC‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCAt         <  1:216011/150‑1 (MQ=255)                                                                                                                                                       tttATCCCCTTTATC‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCAtt        <  1:254440/150‑1 (MQ=255)                                                                                                                                                       tttATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTg       >  1:18386/1‑151 (MQ=255)                                                                                                                                                             aCCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGAtttt  >  2:1042839/1‑150 (MQ=255)                                                                                                                                                             |                                                                                                                                                        ATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAA‑TACGGCA‑TCCGTTCGCAATCCTGGTTTGGCTT‑CGCC‑A‑TT‑AATAAGCCAGGTTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTT  >  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 < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑