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

TCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCA‑TTTTAAATTATTACCGGTGAAATC‑AAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAAT‑AAGCCA‑GGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTT  >  NC_000913/4537961‑4538258                                                                                                                                                          |                                                                                                                                                      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:391614/1‑151 (MQ=255)                                                                                                              aaTCCTGGTTTGGCTTCGCC‑ATTAAT‑AAGCCA‑GGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaat                                          >  2:145224/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:505033/1‑149 (MQ=255)                                                                                                                                          aaGCCA‑GGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAAt                  >  2:579641/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)                                                                                                                                                   ttttATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGc           <  2:309342/148‑1 (MQ=255)                                                                                                                                                   ttttATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCa          <  1:579544/149‑1 (MQ=255)                                                                                                                                                     ttATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTg       >  1:805971/1‑150 (MQ=255)                                                                                                                                                      tATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGAt     >  1:935336/1‑151 (MQ=255)                                                                                                                                                         cACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGAtttt  >  1:1100545/1‑151 (MQ=255)                                                                                                                                                          |                                                                                                                                                      TCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCA‑TTTTAAATTATTACCGGTGAAATC‑AAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCC‑ATTAAT‑AAGCCA‑GGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTT  >  NC_000913/4537961‑4538258

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: ‑