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

ATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTT‑CGCAATCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTA‑TTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTT  >  NC_000913/4537960‑4538258                                                                                                                                                          |                                                                                                                                                        aTCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTT‑CGCAATCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGG‑TTTTATCa                                                                                                                                                         <  1:501380/150‑1 (MQ=255)  tCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTT‑CGCAATCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGG‑TTTTATCAcc                                                                                                                                                       >  1:561595/1‑151 (MQ=255)   ccccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTT‑CGCAATCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGG‑TTTTATCa                                                                                                                                                         <  2:555218/148‑1 (MQ=255)     ccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTT‑CGCAATCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGG‑TTTTATCACC‑ttt                                                                                                                                                   <  1:671487/151‑1 (MQ=255)      cGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTT‑CGCAATCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTa                                                                                                                                                  <  1:206634/151‑1 (MQ=255)         aaaaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTT‑CGCAATCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTA‑tt                                                                                                                                               <  2:340244/150‑1 (MQ=255)            aCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTT‑CGCAATCCTGGTTTGG‑C‑TTAGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTA‑TTca                                                                                                                                             >  1:135105/1‑149 (MQ=255)             ccGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTT‑CGCAATCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTA‑TTcaca                                                                                                                                           >  2:335848/1‑150 (MQ=255)              cGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTT‑CGCAATCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTA‑TTCACAAc                                                                                                                                         <  2:949920/151‑1 (MQ=255)                ccAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTT‑CGCAATCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTA‑TTCACAAccg                                                                                                                                       >  1:790900/1‑151 (MQ=255)                ccAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTT‑CGCAATCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTA‑TTCACAAccg                                                                                                                                       <  1:679486/151‑1 (MQ=255)                ccAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTT‑CGCAATCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTA‑TTCACAAccg                                                                                                                                       >  1:676595/1‑151 (MQ=255)                ccAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTT‑CGCAATCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTA‑TTCACAAccg                                                                                                                                       <  1:41981/151‑1 (MQ=255)                      aCGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTT‑CGCAATCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTA‑TTCACAACCGCCGCAc                                                                                                                                 <  1:95250/151‑1 (MQ=255)                       cGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTT‑CGCAATCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTA‑TTCACAACCGCCGCAc                                                                                                                                 <  2:42406/150‑1 (MQ=255)                          cATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTGCCGGTGAAATCAAGTTCAAATACGGCATCCGTT‑CGCAATCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTA‑TTCACAACCGCCGCACCAGc                                                                                                                             >  1:886438/1‑151 (MQ=255)                                     gagaGGGGAGCAAG‑TTTTTTCATTTTAAATTTTTTCCGGTTAAAATAAGGTACAATTCCGCCTCCCTTCCGCATTCTTGTTTTGC‑C‑TTCCCCATTAAAAAGCCAGT‑TTTTATCCCC‑TTTT‑TTCACAACCGCCGCACCAGCCGTTCCGTAc                                                                                                                   <  2:885875/150‑1 (MQ=255)                                     gagaCGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTT‑CGCAATCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTA‑TTCACAACCGCCGCACCAGCCGTTCCGTAcc                                                                                                                  >  2:650507/1‑151 (MQ=255)                                     gagaCGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTT‑CGCAATCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTA‑TTCACAACCGCCGCACCAGCCGTTCCGTAcc                                                                                                                  >  1:404050/1‑151 (MQ=255)                                          ggaaccaactttATCCAATTTAAAATTTTCCCCGGTAAATCCAATTTCAAAAACGGAACCCTTC‑CGCAACCCTGTTTTGT‑C‑CCCGCAATAAATAACCCAGT‑TTTTATC‑CC‑TTTATTTCACAACCGCCGCACCAGCCGTTCCGTACCAggg                                                                                                              <  2:38552/140‑1 (MQ=255)                                           ggAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTT‑CGCAATCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTA‑TTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCg                                                                                                            >  2:91190/1‑151 (MQ=255)                                               cAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTT‑CGCAATCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTA‑TTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTc                                                                                                        >  1:36579/1‑151 (MQ=255)                                                aaGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTT‑CGCAATCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTA‑TTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCg                                                                                                       >  1:636439/1‑151 (MQ=255)                                                    tttttCCATTTTAAATTATTACCGGTAAAATAAAGTTAAAATCCGGCATCCTTC‑CACAACCCTGGTTGGG‑C‑TCCGCCATTAAAAAGCCGGT‑TTTTACCCCT‑TTTA‑TCCCCACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCa                                                                                                    <  1:410424/147‑1 (MQ=255)                                                       tcccATTTTAAATTATTCCCGGTAAAACCAAGTTAAAATCCGCCACCCGTC‑CCCAACCCGGTTTTGG‑T‑TCCGCAATTAAAAGGCGAGT‑TTTTACCCCT‑TTTA‑TCCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGc                                                                                                  <  2:441987/147‑1 (MQ=255)                                                         ccATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTT‑CGCAATCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTA‑TTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTa                                                                                               >  1:355266/1‑150 (MQ=255)                                                           aTTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTT‑CGCAATCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTA‑TTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACt                                                                                            >  1:319180/1‑151 (MQ=255)                                                            tttaAAATTTTTCCGGGTAAAATAAATTTAAAATCGGCCACCGGTC‑CGCATTCCTGTTTTGG‑C‑TTCCCCATAAATAACCAAGT‑TTTTACCCCT‑TTTA‑TCCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAAc                                                                                             <  2:822157/149‑1 (MQ=255)                                                              ttAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTT‑CGCAATCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTA‑TTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTcc                                                                                          >  1:922154/1‑150 (MQ=255)                                                                aaaTTTTTCCCGGTAAAATCAAGTTAAAATCCGCCACCCTTT‑CGCAATCCGGTTTTGG‑C‑TTCGCCATAAATAAGCCAGG‑TTTTATCCCC‑TTTA‑TTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCAt                                                                                        <  1:933147/150‑1 (MQ=255)                                                                   ttattaCGGGTGAAACCAAGTTAAAATCCGGCCTCCTTT‑CGCAATCCTGTTTTGC‑C‑TTCGCCATTAATAAGCCGGT‑TTTTACCCCT‑TTTT‑TTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGGTTCGCCAGCGTAACTCCATTGc                                                                                     <  1:762353/150‑1 (MQ=255)                                                                        ccgggTAAAATCAATTTAAAATCCGCCACCCGTC‑CGCAACCTTGTTTTGG‑C‑TCCGCAATTAAAAAGCCGGT‑TTTTACCCCC‑TTTA‑TCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtgtg                                                                                 <  2:1067822/146‑1 (MQ=255)                                                                        aCCGGTGAAATCAAGTTCAAATACGGCATCCGTT‑CGCAATCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTA‑TTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGtt                                                                               >  2:332419/1‑151 (MQ=255)                                                                           gggaaaaTCAAGTACAAATGCGACACCCTTC‑CGCAATCCTGGTTTGG‑C‑TTCCCCATTAAAAAGCCAGT‑TTTTACCCCT‑TTTA‑TCCCAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtgtg                                                                                 <  2:276451/142‑1 (MQ=255)                                                                           ggTAAAATCAAGTTAAAATCCGCCACCCGTT‑CGCACTCCTGTTTTGG‑C‑TCCGCCATTAATAAGCCGGT‑TTTTATCCCC‑TTTA‑TTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGa                                                                             <  1:585224/150‑1 (MQ=255)                                                                           ggTAAAATAAAGTTAAAATCCGCCTTCCGTC‑CCCAATCCT‑GTTTGG‑CTTTCCCCATTAAAAAGCCAGT‑TTTTATCCCT‑TTTA‑TCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGa                                                                             <  2:53201/150‑1 (MQ=255)                                                                              gAAATCAAGTTCAAATACGGCATCCGTT‑CGCAATCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTA‑TTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAggg                                                                          >  1:1060750/1‑150 (MQ=255)                                                                               aattcatttacAAATACGCCATCCGTT‑CGAACCCGTGTTTTGG‑C‑TCCGCTATTAATAAGCCAGG‑TTTTTTCCCC‑TTTA‑TTCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGAt                                                                        <  1:490744/141‑1 (MQ=255)                                                                               aaaTCAAGTTAAAATACGGCACCCGTT‑CGCAATCCGGGTTTGG‑C‑TTCGCAATAAAAAAGCCAGT‑TTTTACCCCC‑TTTA‑TTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGAt                                                                        <  1:767748/151‑1 (MQ=255)                                                                               aaaTCAAGTTAAAATACGGCACCCGTT‑CGAATTCCTGGTTTGC‑C‑TTCGCCATTAATAAGCCAGT‑TTTTACCCCC‑TTTA‑TTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGAt                                                                        <  1:601834/151‑1 (MQ=255)                                                                                        ttaaaaACCGCCTCCGTT‑CCCAAACCCGGTTTGGCC‑TCCGCCA‑TAAAAAGCCCGGTTTTTACCCCT‑TTTA‑TCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGa                                                                 <  2:335876/148‑1 (MQ=255)                                                                                         cAAATACGGCATCCGTT‑CGCAATCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTA‑TTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGaca                                                               >  1:996211/1‑150 (MQ=255)                                                                                          aaaTCCGCCATCCGTT‑CGCAATCCTTTTTTGG‑T‑TTCGCCATTAATAAGCCAGT‑TTTTATCCCC‑TTTA‑TCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAg                                                                  <  1:272737/146‑1 (MQ=255)                                                                                           aatccgccttCCTTT‑CGCAA‑CCTGGTTTGT‑CTTTCCCCATTAA‑AAGCCCGG‑TTTTAACCCCTTTTA‑TCCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGaca                                                               <  2:1002256/139‑1 (MQ=255)                                                                                           aaTACGGCATCCGTT‑CGCAATCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTA‑TTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaac                                                             >  2:106091/1‑150 (MQ=255)                                                                                            aTACGGCATCCGTT‑CGCAATCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTA‑TTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaac                                                             >  2:221748/1‑149 (MQ=255)                                                                                                ggCACCCGTT‑CGCAATCCGGTTTTGG‑C‑TTCGCAATTAAAAAGCCAGG‑TTTTATCCCC‑TTTA‑TCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGaaa                                                        <  1:21173/150‑1 (MQ=255)                                                                                                gcctcccGTT‑CGAAATCCGGTTTTGG‑T‑TTCGCCATAAATAAGCCAGG‑TTTTACCCCC‑TTTA‑TCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAg                                                           <  1:949975/142‑1 (MQ=255)                                                                                                  cATCCGTC‑CGCATTCCGTGTTTGG‑C‑TTCGCCATTAATAACCCAGT‑TTTTATCCCC‑TTTA‑TTCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAAc                                                       <  2:736913/149‑1 (MQ=255)                                                                                                   acccttt‑CGCACCCCTGGTTGGG‑C‑TTCCCCATTAATAACCCAGT‑TTTTACCCCC‑TTTA‑TCCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGaa                                                         <  2:887508/141‑1 (MQ=255)                                                                                                    tCCGTT‑CGCAATCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTA‑TTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTAt                                                   >  1:258291/1‑151 (MQ=255)                                                                                                      ccttccgaaccccTGGTTTGG‑C‑TCCGCCTTAAATAACCCGGT‑TTTTACCCCC‑TTTT‑TCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGa                                                          <  2:19969/132‑1 (MQ=255)                                                                                                      cGTT‑CGCAATCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTA‑TTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTg                                                 >  1:133255/1‑151 (MQ=255)                                                                                                        tt‑CGCAATCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTA‑TTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTga                                                >  2:762826/1‑150 (MQ=255)                                                                                                        tt‑CGCAATCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTA‑TTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTga                                                >  1:915719/1‑150 (MQ=255)                                                                                                          gcGCAATCCTGTTTTGT‑C‑TCCGCAATAAAAACGCCAGG‑TTTTATCACT‑TTTA‑TCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaa                                               <  2:943230/149‑1 (MQ=255)                                                                                                           cgaaaccgggttttGG‑C‑TTCCCCATAAAAAGGCCGGT‑TTTTATCCCT‑TTTA‑TTCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTg                                                 <  1:597512/136‑1 (MQ=255)                                                                                                           cGCAATCCTGTTTGGG‑T‑TCCGCCATAAATAGGCCGGG‑TTTTACCCCT‑TTTT‑TCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaag                                              <  2:119641/150‑1 (MQ=255)                                                                                                              aaTCCTGTTTTGT‑T‑TTCCCCATTAAAAAGCCAGT‑TTTTATCACT‑TTTA‑TTCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaag                                              <  1:552335/147‑1 (MQ=255)                                                                                                              aaTCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTA‑TTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaat                                          >  1:94302/1‑151 (MQ=255)                                                                                                               aTCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTA‑TTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaat                                          >  1:640006/1‑150 (MQ=255)                                                                                                                    ggTTGGG‑C‑TCCGCCA‑TAAAAAGCCCGGTTTTTACCCCC‑TTTA‑CCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtctt                                    <  2:202272/151‑1 (MQ=255)                                                                                                                        tGG‑C‑TTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTA‑TTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGc                                >  1:482448/1‑151 (MQ=255)                                                                                                                              ttCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTA‑TTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttt                            >  2:889413/1‑151 (MQ=255)                                                                                                                              tcccccATTAAAAAGCCGGT‑TTTTATCCCT‑TTTA‑TCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtt                              <  2:263792/145‑1 (MQ=255)                                                                                                                                ccccATTAAAAAGCCAGT‑TTTTACCCCT‑TTTA‑TTCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttttt                          <  1:660635/149‑1 (MQ=255)                                                                                                                                  ccATTAATAAGCCAGG‑TTTTATCACC‑TTTA‑TTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGt                        >  2:726545/1‑151 (MQ=255)                                                                                                                                  ccATAAATAGCCCAGT‑TTTTACCCCC‑TTTA‑CCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGt                        <  2:274298/151‑1 (MQ=255)                                                                                                                                  caaTTAAAAAGCCTGT‑TTTTACCTCT‑TTTA‑CCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTg                         <  2:567332/148‑1 (MQ=255)                                                                                                                                    aTTAATAAGCCAGG‑TTTTATCACC‑TTTA‑TTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTc                       >  2:653822/1‑150 (MQ=255)                                                                                                                                     ttAATAAGCCAGG‑TTTTATCACC‑TTTA‑TTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCa                      >  2:88735/1‑150 (MQ=255)                                                                                                                                      taataaGCCAGG‑TTTTATCACC‑TTTA‑TTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaa                     >  2:2941/1‑150 (MQ=255)                                                                                                                                       aatacGCCAGG‑TTTTACCCCC‑TTTA‑TCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGt                        <  2:123214/146‑1 (MQ=255)                                                                                                                                       aataaCCCAGT‑TTTTATCCCT‑TTTA‑TCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCttt                             <  1:791802/141‑1 (MQ=255)                                                                                                                                          aaGCCAGG‑TTTTATCACC‑TTTA‑TTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAAt                  >  2:704385/1‑149 (MQ=255)                                                                                                                                            gCCAGG‑TTTTATCACC‑TTTA‑TTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTg              >  1:116060/1‑151 (MQ=255)                                                                                                                                              cgtgtttttCTCCCC‑TTTA‑TCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATa                 <  1:920578/142‑1 (MQ=255)                                                                                                                                                gg‑TTTTATCACC‑TTTA‑TTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGc           >  1:554263/1‑150 (MQ=255)                                                                                                                                                gg‑TTTTACCCCT‑TTTA‑TCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGc           <  2:453476/150‑1 (MQ=255)                                                                                                                                                     ttATCCCT‑TTTA‑TCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTg       <  2:403943/150‑1 (MQ=255)                                                                                                                                                         cACT‑TTTA‑TCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGAttt   <  1:635880/150‑1 (MQ=255)                                                                                                                                                         cACC‑TTTA‑TTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGAtttt  >  1:252220/1‑151 (MQ=255)                                                                                                                                                          aCC‑TTTA‑TTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGAtttt  >  2:107870/1‑150 (MQ=255)                                                                                                                                                          |                                                                                                                                                        ATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTT‑CGCAATCCTGGTTTGG‑C‑TTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTA‑TTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTT  >  NC_000913/4537960‑4538258

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