Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,538,101	0	G	T	40.6%	122.1 / 36.1	64	T222N (ACC→AAC)	nanM	N‑acetylneuraminic acid mutarotase
Reads supporting (aligned to +/- strand):  ref base G (31/7);  new base T (0/26);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 7.38e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGCTTATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCA‑TTTTAAATTATTACCGGTGAAATC‑AAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAAT‑AAGCCA‑GGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGC  >  NC_000913/4537954‑4538249                                                                                                                                                         |                                                                                                                                                     ttGCTTATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCA‑TTTTAAATTATTACCGGTGAAATC‑AAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAAT‑AAGCCA‑GGttt                                                                                                                                                   >  2:787782/1‑151 (MQ=255)  tGCTTATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCA‑TTTTAAATTATTACCGGTGAAATC‑AAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAAT‑AAGCCA‑GGttt                                                                                                                                                   <  2:703245/150‑1 (MQ=255)      tATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCA‑TTTTAAATTATTACCGGTGAAATC‑AAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAAT‑AAGCCA‑GGTTTTATc                                                                                                                                               >  2:877848/1‑150 (MQ=255)        tCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCA‑TTTTAAATTATTACCGGTGAAATC‑AAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAAT‑AAGCCA‑GGTTTTATCAcc                                                                                                                                            >  2:857051/1‑151 (MQ=255)         ccccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCA‑TTTTAAATTATTACCGGTGAAATC‑AAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAAT‑AAGCCA‑GGTTTTATCACCt                                                                                                                                           <  2:972662/151‑1 (MQ=255)          cccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCA‑TTTTAAATTATTACCGGTGAAATC‑AAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAAT‑AAGCCA‑GGTTTTATCACCtt                                                                                                                                          >  1:53839/1‑151 (MQ=255)                aaaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCA‑TTTTAAATTATTACCGGTGAAATC‑AAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAAT‑AAGCCA‑GGTTTTATCACCTTTAtt                                                                                                                                      >  1:643354/1‑149 (MQ=255)                    cGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCA‑TTTTAAATTATTACCGGTGAAATC‑AAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAAT‑AAGCCA‑GGTTTTATCACCTTTATTCACAAc                                                                                                                                >  1:288625/1‑151 (MQ=255)                     gCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCA‑TTTTAAATTATTACCGGTGAAATC‑AAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAAT‑AAGCCA‑GGTTTTATCACCTTTATTCACAAcc                                                                                                                               >  1:698383/1‑151 (MQ=255)                      ccAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCA‑TTTTAAATTATTACCGGTGAAATC‑AAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAAT‑AAGCCA‑GGTTTTATCACCTTTATTCACAAccg                                                                                                                              <  1:335046/151‑1 (MQ=255)                             ccccATCGGGTGATGAGACGGGAGCAAGTTTATTCCA‑TTTAAAATTATTCCCGGTAAAATC‑AAGTTAAAATCCGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAAT‑AAGCCA‑GTTTTTACCCCCTTTTTCCACAACCGCCGCAcc                                                                                                                       <  1:834937/149‑1 (MQ=255)                             cGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCA‑TTTTAAATTATTACCGGTGAAATC‑AAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAAT‑AAGCCA‑GGTTTTATCACCTTTATTCACAACCGCCGCAc                                                                                                                        <  1:1050417/150‑1 (MQ=255)                                    tGGTGATGAGACGGGAGCAAGCTTATTCCA‑TTTTAAATTATTACCGGTGAAATC‑AAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAAT‑AAGCCA‑GGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGtt                                                                                                                >  2:626546/1‑151 (MQ=255)                                                 ggAGCAAGTTTTTTCAT‑TTTAAAATTATTCCCGGTGAATTA‑AAGTTCAAATACGGCATCCGTCCGCAATCCTGGTTTGGTTTCGCC‑ATTAAT‑AAGCCA‑GTTTTTACCCCCTTTATTCCCAACCGCCGCACCAGCCGTTCCGTACCAGGGCg                                                                                                   <  2:340981/151‑1 (MQ=255)                                                 ggAGAAAGCTTATCCCT‑TTTTAATTGTTTCCCGGTAAAATC‑AGGTTAAAATACGGCTTCCGTCCGAACTCCTGGTTTGGCTTCGCC‑ATTAAT‑AAGCGG‑GTTTTTACCCCCTTTTTCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCg                                                                                                   <  2:1113846/151‑1 (MQ=255)                                                 ggAACAAGTTTTTTCCTTTTTTAAATTATTCCCGGGAAAATA‑AATTTAAATTACGCCTCCCTTCCGCACTCCTGTTTTGTCTCCGCA‑ATTAAT‑AGCCCG‑GTTTTTATCCCTTTTATCCAAACCCGCCGCACCAGCCGTTCCGTACCAgg                                                                                                      <  1:50567/149‑1 (MQ=255)                                                   aGCAAGCTTATTCCA‑TTTTAAATTATTACCGGTGAAATC‑AAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAAT‑AAGCCA‑GGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGAt                                                                                                 >  2:130164/1‑151 (MQ=255)                                                     cAAGCTTATTCCA‑TTTTAAATTATTACCGGTGAAATC‑AAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAAT‑AAGCCA‑GGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTc                                                                                               >  2:949658/1‑151 (MQ=255)                                                       aagttttttcattttaAAATTATTACGGGTAAAATA‑AAGTAAAAAACCGCCATCCGTTCGCATCCCGGGTTGGGCTCCGCC‑ATAAAT‑AAGCCA‑GTTTTTACCACCTTTATCCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTc                                                                                               <  2:609871/138‑1 (MQ=255)                                                       aGCTTATTCCA‑TTTTAAATTATTACCGGTGAAATC‑AAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAAT‑AAGCCA‑GGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGc                                                                                             >  1:94370/1‑151 (MQ=255)                                                           tttcccA‑TTTAAAATTATT‑CCCGTGAAAACAAAGTCCAAAGACGGCTTCCGTCCACAATCCTGTTTTGGTTGCTCA‑ATAAAA‑AACCAA‑GGTTTTATCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCa                                                                                           <  2:1081779/145‑1 (MQ=255)                                                               ccA‑TTTTAAATTATTACCGGTGAAATC‑AAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAAT‑AAGCCA‑GGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTa                                                                                      >  1:443626/1‑150 (MQ=255)                                                                     ttAAATTATTACCGGTGAAATC‑AAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAAT‑AAGCCA‑GGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTcc                                                                                 >  2:701976/1‑150 (MQ=255)                                                                         attattACCGGTGAAATC‑AAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAAT‑AAGCCA‑GGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGc                                                                            >  1:19976/1‑151 (MQ=255)                                                                           tattaCCGGTGAAATC‑AAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAAT‑AAGCCA‑GGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtg                                                                          >  1:602341/1‑151 (MQ=255)                                                                             ttaCCGGTGAAATC‑AAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAAT‑AAGCCA‑GGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtgtg                                                                        >  2:144702/1‑151 (MQ=255)                                                                                 cGGTGAAATC‑AAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAAT‑AAGCCA‑GGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTg                                                                     >  2:737104/1‑150 (MQ=255)                                                                                      aaaTC‑AAGTTAAAATACGCCATCCGTTCGCAATCCTGTTTTGGCTTCGCC‑ATTAAT‑AAGCCA‑GTTTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGa                                                                <  2:366152/150‑1 (MQ=255)                                                                                        aTC‑AAGTTCAAATCCGGCCCCCGTCCCCAACCCGGGTTGGGCTGCCCC‑ATTAAT‑AACCCA‑GTTTTTACCCCTTTTATCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATc                                                              <  1:298558/150‑1 (MQ=255)                                                                                          c‑AAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAAT‑AAGCCA‑GGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaa                                                            >  2:1086268/1‑150 (MQ=255)                                                                                            aattaaaaaTACGCGTTCCGTTCGCAACCCTGTTTGGCCTTCCCC‑ATAAAT‑AAGCAA‑GTTTTTACCTCTTTTATCCAAATCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaaaa                                                          <  1:1064066/145‑1 (MQ=255)                                                                                            aaGTTAAAATCCGCAATCCTTTCGAATCCCTGGTTTGCCTTCGCC‑ATTATT‑ACGCCA‑GTTTTTACCACCTTTATTCCCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaaa                                                           <  2:232685/150‑1 (MQ=255)                                                                                               ttCAAATCCGGCACCCTTTCGC‑ATCCCTGTTTTGCTTCCCCAATTAAT‑AACCCA‑GTTTTTACCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAg                                                         <  2:150706/149‑1 (MQ=255)                                                                                               ttCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAAT‑AAGCCA‑GGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGac                                                       >  2:155705/1‑151 (MQ=255)                                                                                                   aaTACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAAT‑AAGCCA‑GTTTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaaca                                                   <  2:991293/151‑1 (MQ=255)                                                                                                   aaTACGGCACCCGTTCGAAATCCTGTTTGGGCTTCGCC‑ATTAAA‑AAGCCA‑GTTTTTACCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaaca                                                   <  2:1091436/151‑1 (MQ=255)                                                                                                             ccGTCCCCATTCCTGTTTTGGCTTCCCC‑ATAAAA‑AAGCCA‑GTTTTTACCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTAt                                          <  1:68990/150‑1 (MQ=255)                                                                                                                    cactcttgtttgggCTTCGCA‑ATTAAA‑AAGCCA‑GTTTTTATCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTa                                           <  1:186533/130‑1 (MQ=255)                                                                                                                    cAATCCTGGTTTGGCTTCGCC‑ATTAAT‑AAGCCA‑GGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaagaa                                   >  1:866048/1‑150 (MQ=255)                                                                                                                     aaTCTTGTTTTGGTTGCGCT‑ATAAAA‑AAGCCA‑GGTTTTATCACCTTTATCCCAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGa                                                 <  1:78769/135‑1 (MQ=255)                                                                                                                     aaTCCTGGTTTGGCTTCGCC‑ATTAAT‑AAGCCA‑GGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaat                                 >  2:145224/1‑151 (MQ=255)                                                                                                                     aaTCCTGGTTTGGCTTCGCC‑ATTAAT‑AAGCCA‑GGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaat                                 >  2:391614/1‑151 (MQ=255)                                                                                                                      aTCCTGGTTTGGCTTCGCC‑ATTAAT‑AAGCCA‑GGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaat                                 >  1:848048/1‑150 (MQ=255)                                                                                                                                ggttccccc‑ATAAAT‑AGGCCT‑GTTTTCATCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTg                        <  1:204043/142‑1 (MQ=255)                                                                                                                                ggCTTCGCC‑ATTAAT‑AACCCA‑GTTTTTATCCCCTTTATCCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGc                       <  2:435621/150‑1 (MQ=255)                                                                                                                                ggCTCCCCC‑ATCAAT‑AGGCCT‑GTTTTTAACCCCTTTACTACCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAAGAATCTTCTg                        <  1:76491/149‑1 (MQ=255)                                                                                                                                   ttCGCC‑ATTAAT‑AAGCCA‑GGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttt                   >  1:557522/1‑151 (MQ=255)                                                                                                                                    tcccc‑ATTAAAAAAGCCG‑GTTTTTACCCCCTTTTTCCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtt                     <  1:1073878/146‑1 (MQ=255)                                                                                                                                     cGCC‑ATTAAT‑AAGCCA‑GTTTTTACCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCt                      <  2:94380/146‑1 (MQ=255)                                                                                                                                     cGCC‑ATAAAA‑AACCCA‑GTTTTTACCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttttt                 <  2:602438/151‑1 (MQ=255)                                                                                                                                     cGCA‑ATTAAA‑AAGCCG‑GTTTTTACCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttttt                 <  1:689954/151‑1 (MQ=255)                                                                                                                                      gCC‑ATTAAT‑AAGCCA‑GGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTg                >  1:884817/1‑151 (MQ=255)                                                                                                                                      gCC‑ATTAAT‑AAGCCA‑GGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTg                >  2:519670/1‑151 (MQ=255)                                                                                                                                       ccccttaaaaaacccA‑GTTTTTACCCCTTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttttt                 <  2:576120/138‑1 (MQ=255)                                                                                                                                        cccataaa‑aaGCCATGTTTTTATCCCTTTAATTCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttttt                 <  1:378752/146‑1 (MQ=255)                                                                                                                                         caataaaaaaGCCA‑GTTTTTATCCCTTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGt               <  1:13516/143‑1 (MQ=255)                                                                                                                                             aat‑aaGCCA‑GGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaaa           >  1:887438/1‑150 (MQ=255)                                                                                                                                             aat‑aaCCCG‑GTTTTTATCCCTTTTATCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaa            <  1:626443/149‑1 (MQ=255)                                                                                                                                               t‑aaGCCA‑GGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAAt         >  1:842469/1‑150 (MQ=255)                                                                                                                                                 aaGCCA‑GGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAAt         >  2:579641/1‑149 (MQ=255)                                                                                                                                                 aaGCCA‑GGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAAt         >  2:505033/1‑149 (MQ=255)                                                                                                                                                 aaGCCA‑GGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAAt         >  2:964811/1‑149 (MQ=255)                                                                                                                                                    ccA‑GTTTTTACCCCCTTAATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGt      <  1:180931/149‑1 (MQ=255)                                                                                                                                                    ccA‑GGTTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGt      <  2:683904/149‑1 (MQ=255)                                                                                                                                                        ggTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGc  >  1:482049/1‑150 (MQ=255)                                                                                                                                                        ggTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGc  >  2:723841/1‑150 (MQ=255)                                                                                                                                                         |                                                                                                                                                     TTGCTTATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCA‑TTTTAAATTATTACCGGTGAAATC‑AAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAAT‑AAGCCA‑GGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGC  >  NC_000913/4537954‑4538249

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