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

TCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCG‑GTGAAATCAAGTTCAAATACGGCA‑TCCGTTCGCAA‑T‑CCTGGTTTGGCTTCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTT  >  NC_000913/4537961‑4538258                                                                                                                                                             |                                                                                                                                                        tCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCG‑GTGAAATCAAGTTCAAATACGGCA‑TCCGTTCGCAA‑T‑CCTGGTTTGGCTTCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ATCAcc                                                                                                                                                       >  1:681838/1‑151 (MQ=255)   cccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCG‑GTGAAATCAAGTTCAAATACGGCA‑TCCGTTCGCAA‑T‑CCTGGTTTGGCTTCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ATCACC‑tt                                                                                                                                                    <  2:353677/151‑1 (MQ=255)      gCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCG‑GTGAAATCAAGTTCAAATACGGCA‑TCCGTTCGCAA‑T‑CCTGGTTTGGCTTCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ATCACC‑TTTAt                                                                                                                                                 >  1:724170/1‑151 (MQ=255)       cAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCG‑GTGAAATCAAGTTCAAATACGGCA‑TCCGTTCGCAA‑T‑CCTGGTTTGGCTTCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ATCACC‑TTTAtt                                                                                                                                                <  2:90716/151‑1 (MQ=255)        aaaaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCG‑GTGAAATCAAGTTCAAATACGGCA‑TCCGTTCGCAA‑T‑CCTGGTTTGGCTTCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ATCACC‑TTTATT‑c                                                                                                                                              >  2:255976/1‑151 (MQ=255)         agccccccAGCTACGCCATCTGGTGATGAGACGAGAGCAAGCTTAATCCATTTTAAATTAATACCG‑GTGAAATCAAGTTCAAATGCGGGA‑TCCGTTGGCAA‑T‑CCTGGTTTGGCTCCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ACCGCC‑TTTAtt                                                                                                                                                <  2:218362/143‑1 (MQ=255)         aaaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCG‑GTGAAATCAAGTTCAAATACGGCA‑TCCGTTCGCAA‑T‑CCTGGTTTGGCTTCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ATCACC‑TTTAtt                                                                                                                                                <  2:821711/149‑1 (MQ=255)          aaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCG‑GTGAAATCAAGTTCAAATACGGCA‑TCCGTTCGCAA‑T‑CCTGGTTTGGCTTCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ATCACC‑TTTATT‑cac                                                                                                                                            <  2:807269/151‑1 (MQ=255)             cGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCG‑GTGAAATCAAGTTCAAATACGGCA‑TCCGTTCGCAA‑T‑CCTGGTTTGGCTTCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ATCACC‑TTTATT‑CACAAc                                                                                                                                         >  2:452420/1‑151 (MQ=255)                 aGCTACGCCATCTGGTGATGAGACGGGAGCGAGCTTATTCCATTTTAAATTATTACCG‑GTGAAATCAAGTTCAAATACGGCA‑TCCGTTCGCAA‑T‑CCTGGTTTGGCTTCGCC‑ATT‑AATAAGCCAGG‑TTTT‑GTCACC‑TTTATT‑CACAAccgcc                                                                                                                                     <  1:79632/151‑1 (MQ=255)                     aCGCCATCTGGGGATGAGACGGGAGCAAGTTTATTCCATTTTAAATTATTCCCG‑GTGAAATAAATTTAAAATACGGCA‑TCCTTTCGCAA‑T‑CCTGTTTTGGCTTCGCC‑ATT‑AAAAAGCCAGT‑TTTT‑ATCCCC‑TTTATT‑CACAACcgccgc                                                                                                                                   <  2:497405/149‑1 (MQ=255)                        ccATCGGGTGAGGAGACGGGAGAAAGTTTTTTCCTTTTAAAATTATTCCCG‑GTAAAATCAAGTTAAAATCCGGCA‑TCCGTTCGCAA‑C‑CCTGGTTTGGCTTCGCC‑ATT‑AATAAGCCGGG‑TTTT‑ATCCCT‑TTTATC‑CACAACCGCCGCACCAg                                                                                                                              <  2:774871/151‑1 (MQ=255)                             tGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCG‑GTGAAATCAAGTTCAAATACGGCA‑TCCGTTCGCAA‑T‑CCTGGTTTGGCTTCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGtt                                                                                                                         >  2:4208/1‑151 (MQ=255)                              ggTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCG‑GTGAAATCAATTTCAAATACGGCA‑TCCGTTCGCAA‑T‑CCTGGTTTGGCTTCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ATCCCC‑TTTATT‑CACAACCGCCGCACCAGcc                                                                                                                            <  1:160786/147‑1 (MQ=255)                                  aggggaCGGCAGCAATCTTTTTCACTTTAAAATTATTACCG‑GTGAAAACAA‑TTACAATATCGCACTCCTTTC‑CAA‑TCCCTGGTTTGGCTCCCCC‑ATA‑AATAA‑CCAGG‑TTTTCATCCCC‑TTTATC‑AAAACCCGCCGCACCAGCCGTTCCGTa                                                                                                                    <  2:777998/147‑1 (MQ=255)                                               aaGTTTTTTCCATTTAAAATTATTCCGG‑GTAAAATCAATTTAAAATCCGGCT‑TCCGTCCGCAA‑T‑CCTGTTTTGGCTTCCCA‑ATT‑AATAAGCCAGT‑TTTT‑ATCCCC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTc                                                                                                        <  1:420807/150‑1 (MQ=255)                                               aaGCTTTTTCAATTTAAATTTTTTCCCG‑GTAAAATAAATTTAAAATCCGGCC‑CCCGTGCGCAA‑C‑CCGGTTTTGGCTTCCCC‑ATT‑AATAAGCCGGT‑TTTT‑ATCCCT‑TTTTTT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTc                                                                                                        <  2:940522/150‑1 (MQ=255)                                                 gCTTATTCCATTTTAAATTATTACCG‑GTGAAATCAAGTTCAAATACGGCA‑TCCGTTCGCAA‑T‑CCTGGTTTGGCTTCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc                                                                                                     >  1:207964/1‑151 (MQ=255)                                                   ttttttttttttAAATTTTTCCCGTGTAAAATCAAGTCAAAATACGCCC‑TCCTTCCCCAT‑T‑CCT‑GTTTGGTTTCCCCAATA‑AATCAGCGTGT‑TTTT‑ACCCCT‑TTTATC‑CAAAGCCGCCGCACCAGCCGTTCCGTACCAGGGCGAt                                                                                                          <  2:12589/137‑1 (MQ=255)                                                    tttcctttttaAAATTATTCCCG‑GTAAAATCAATTTAAATTCCGCCA‑CCCGTTCCCAC‑T‑CCTGTTTTGGC‑TCCCC‑ATTAAAAAAGCCGGG‑TTTT‑ATCCCC‑TTTTTC‑CAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGAtt                                                                                                         <  2:453414/137‑1 (MQ=255)                                                       tCCATTTTAAATTATTACCG‑GTGAAATCAAGTTCAAATACGGCA‑TCCGTTCGCAA‑T‑CCTGGTTTGGCTTCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGt                                                                                                >  2:457325/1‑150 (MQ=255)                                                        ccATTTAAAATTATTCCGG‑TTAAAATAAATTAAAAATAGGCAA‑TCCTTCCGCAC‑C‑CCTGGTTTGGCTCCCCA‑TTT‑AATAAGCCGGT‑TTTT‑ACCCCT‑TTTTTC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGt                                                                                                <  2:361444/149‑1 (MQ=255)                                                          aTTTTAAATTATTACCG‑GTGAAATCAAGTTCAAATACGGCA‑TCCGTTCGCAA‑T‑CCTGGTTTGGCTTCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACt                                                                                            >  2:141823/1‑151 (MQ=255)                                                           tttaAATTTATTCCGG‑GTAAAATAAATTTAAATTACGCCC‑TCCTTCCCCAT‑T‑CTTGTTTTGGCTTCGCA‑ATA‑AAAAAGCCGGG‑TTTT‑ACCCCT‑TTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAAc                                                                                             <  2:870176/149‑1 (MQ=255)                                                            tttAAATTTTTTCCG‑GGTAAATCAAATTAAAAATCGGCC‑TTCCTTCCCAA‑TCCCTTTTTTGGCTCCCCC‑TTA‑AATAACCCAGT‑TTTT‑ACCCCT‑TTTTTC‑CAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTaa                                                                                              <  2:681745/148‑1 (MQ=255)                                                            tttAAATTATTACCG‑GTGAAATCAAGTTCAAATACGGCA‑TCCGTTCGCAA‑T‑CCTGGTTTGGCTTCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTcc                                                                                          >  1:271385/1‑151 (MQ=255)                                                             ttAAATTATTACCG‑GTGAAATCAAGTTCAAATACGGCA‑TCCGTTCGCAA‑T‑CCTGGTTTGGCTTCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTcc                                                                                          >  2:166289/1‑150 (MQ=255)                                                              tAAATTATTACCG‑GTGAAATCAAGTTCAAATACGGCA‑TCCGTTCGCAA‑T‑CCTGGTTTGGCTTCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCAt                                                                                        >  2:366829/1‑151 (MQ=255)                                                               aaaTTATTCCGG‑GTAAAATCAAGTTAAAATCCGGCA‑‑CCCTTCCCAACC‑CCTGGTTTGCCTCCGCA‑ATT‑AAT‑AGCCAGGTTTTT‑ATCCCT‑TTTATC‑AACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCAtt                                                                                       <  2:736369/151‑1 (MQ=255)                                                                      ttcccgggaaaaTCAATTTAAAAT‑CCGCC‑CCCTTTCCCAA‑TCCCTGGTTGGCCTCCGCT‑A‑T‑AAAAAGCCAAG‑TTTT‑AACCCCTTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCt                                                                                    <  2:779739/139‑1 (MQ=255)                                                                                      gTTCAAATACGGCA‑TCCGTTCGCAA‑T‑CCTGGTTTGGCTTCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaaa                                                                    >  1:394166/1‑148 (MQ=255)                                                                                       tcaaaaTACGCCA‑CCCTTCCGAAA‑C‑CCTGGTTGGGCTCCGCC‑ATT‑AAAAAGCCAGT‑TTTT‑ACCCCT‑TTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGa                                                                 <  1:14504/147‑1 (MQ=255)                                                                                          aaattggacA‑TCCGTTCGCAT‑T‑CCGGGTTTGGCTTC‑CC‑CTT‑TATAAGCCCGGTTTTT‑ATCCCT‑TTTTTT‑CAAAACCGCCGCTCCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaa                                                              <  1:105183/142‑1 (MQ=255)                                                                                          aaaTCCGGCA‑TCTGTCCGCAA‑T‑CCTGTTTGGGCTTCGCC‑TTT‑AATAGCCCAGG‑TTTT‑A‑CCCC‑TTTATTCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaa                                                              <  2:402950/150‑1 (MQ=255)                                                                                          aaaTCCGGCA‑TCCTTCCGCAT‑T‑CCGGGTTTGCCTTCCCC‑ATA‑AAAAACCCAGT‑TTTT‑ATCCCC‑TTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaa                                                              <  1:74476/150‑1 (MQ=255)                                                                                                gccA‑TCCGTTCGCAA‑T‑CCTGTTTTGGCTTCGCC‑ATT‑AATAAGCCAGT‑TTTT‑ACCCCT‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGaaa                                                        <  1:739563/148‑1 (MQ=255)                                                                                                 gccttCCGTTCGCAT‑C‑CCTGTTTGGGCTTCGCC‑TTA‑AATAAGCCAGT‑TTTT‑ACCCCC‑TTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGa                                                          <  1:29744/144‑1 (MQ=255)                                                                                                 gccttCCGTTCACAA‑C‑CCTGTTTTGGTTTCCCC‑ATT‑AAAAACCCAGT‑TTTT‑ACCCCC‑TTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGaa                                                         <  1:227102/145‑1 (MQ=255)                                                                                                      ccGTTCGCAA‑T‑CCTGGTTTGGCTTCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTAtt                                                  >  2:166264/1‑151 (MQ=255)                                                                                                         ttCGCAA‑T‑CCTGTTTTGGCTTCGCC‑ATA‑AATAAGCCAGT‑TTTT‑ATCCCC‑TTTTTC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTga                                                <  1:349193/150‑1 (MQ=255)                                                                                                         ttCGCAA‑T‑CCTGGTTTGGCTTCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTga                                                >  2:22950/1‑150 (MQ=255)                                                                                                         ttCGCAA‑T‑CCTGGTTTGGCTTCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTga                                                >  1:657443/1‑150 (MQ=255)                                                                                                         ttCACAA‑T‑CTGGTTTTGGCTCCGCC‑TTT‑AATAAGCCAGG‑TTTT‑ATCACC‑TTTATC‑CAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTg                                                 <  1:768164/149‑1 (MQ=255)                                                                                                         tccGCAA‑T‑CCTGTTTGGTCTTCGCC‑AAT‑AAAAACCCGGT‑TTTT‑ATCCCT‑TTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTg                                                 <  2:559409/147‑1 (MQ=255)                                                                                                                   ctgggTTGGGCTTCGCC‑ATA‑AATAGCCAAGG‑TTTT‑ATCACC‑TTTCTC‑CAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaataa                                        <  2:758783/147‑1 (MQ=255)                                                                                                                     ttgTTAGGCTACCCC‑TTT‑AAAAAGAAACG‑TTTA‑ATCCCT‑TTCCTG‑CACAGCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACACCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaataat                                       <  1:21126/147‑1 (MQ=255)                                                                                                                     tgttttGCCTTCGCC‑ATT‑AATAACCCAGT‑TTTT‑ATCCCC‑TTTATT‑CACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaataat                                       <  1:925747/146‑1 (MQ=255)                                                                                                                      ggTTTGGCTTCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtctt                                    >  1:316631/1‑151 (MQ=255)                                                                                                                        tttGGCTTCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtcttct                                  >  2:350853/1‑151 (MQ=255)                                                                                                                           ggCTTCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCt                               >  1:14090/1‑151 (MQ=255)                                                                                                                              ttCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttt                            >  1:874396/1‑151 (MQ=255)                                                                                                                                tccccaataaataaGCCAGT‑TTTA‑ACCCCT‑TTTACC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtcttct                                  <  2:248455/136‑1 (MQ=255)                                                                                                                                  ccaaTT‑AATAACCCGGT‑TTAA‑ACCACC‑TTTATC‑ACCACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCttttt                           <  2:921941/146‑1 (MQ=255)                                                                                                                                      tt‑AATAAGCCAGT‑TTTT‑ATCCCC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCttttt                           <  1:953879/145‑1 (MQ=255)                                                                                                                                         aataaGCCAGT‑TTTT‑ATCCCC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaa                     <  1:736582/149‑1 (MQ=255)                                                                                                                                             acccAGT‑TTTT‑ACCCCC‑TTTATT‑CCAATCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATa                 <  2:798890/147‑1 (MQ=255)                                                                                                                                                       tt‑ATCCCC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTg       <  1:503342/150‑1 (MQ=255)                                                                                                                                                            cACT‑TTTATC‑CAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGAt     <  2:275119/148‑1 (MQ=255)                                                                                                                                                            cACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGAtttt  >  2:312588/1‑151 (MQ=255)                                                                                                                                                             |                                                                                                                                                        TCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCG‑GTGAAATCAAGTTCAAATACGGCA‑TCCGTTCGCAA‑T‑CCTGGTTTGGCTTCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTT  >  NC_000913/4537961‑4538258

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