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

ATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTT  >  NC_000913/4537960‑4538257                                                                                                                                                            |                                                                                                                                                       aTCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ATCAc                                                                                                                                                       <  2:583916/151‑1 (MQ=255)  tCCCCGCAAAACCACCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ATCAcc                                                                                                                                                      >  2:826758/1‑151 (MQ=255)   ccccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTA‑TCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ATCACC‑t                                                                                                                                                    <  1:796287/150‑1 (MQ=255)    cccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ATCACC‑tt                                                                                                                                                   <  1:258829/151‑1 (MQ=255)      cGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ATCACC‑TTTa                                                                                                                                                 <  1:554609/151‑1 (MQ=255)        cAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ATCACC‑TTTAtt                                                                                                                                               <  2:381467/151‑1 (MQ=255)             ccGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ATCACC‑TTTATT‑CACaa                                                                                                                                         <  2:970076/151‑1 (MQ=255)                ccAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ATCACC‑TTTATT‑CACAAccg                                                                                                                                      <  2:435368/151‑1 (MQ=255)                 cAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ATCACC‑TTTATT‑CACAAccgc                                                                                                                                     >  1:251236/1‑151 (MQ=255)                     tACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ATCACC‑TTTATT‑CACAACcgccgc                                                                                                                                  >  2:55084/1‑150 (MQ=255)                            tCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ATCACC‑TTTATT‑CACAACCGCCGCACCAGCCg                                                                                                                          >  1:928744/1‑151 (MQ=255)                            tCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ATCACC‑TTTATT‑CACAACCGCCGCACCAGCCg                                                                                                                          >  2:516908/1‑151 (MQ=255)                                          gggAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAATACGACATCCGTTCGCTATCTTGGCTTGG‑CTTCAAC‑ATT‑AATAAGCCAGG‑TTTT‑ATCACC‑TTTATT‑CACAACCGCTGCACAATACGTTCTGCACCAGGGc                                                                                                            >  1:587521/1‑151 (MQ=255)                                           ggaagaaGGTTTTTCCCTTTTAAATTTATTACCGGGAAAATGAGGTTCAATCCCGGCTCCTGTCCGCAACCCTGGTTTGC‑CTCCCCC‑TTA‑AAAACGCCGGG‑TTTT‑ATCTCC‑TTTATT‑AACACCCGCCGCACCAGCCGTTCCGTACCAGGGc                                                                                                            <  2:370591/146‑1 (MQ=255)                                           ggAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCg                                                                                                           >  1:209480/1‑151 (MQ=255)                                           ggAGAAAGTTTTTCCCTTTTTAAA‑TTTTTCCCGGTAAAATCAATTACAAATCCGGCACCCTTTCGCATCCCTGGTTTGC‑CTCCCCC‑ATA‑AATAACCCAGG‑TTTT‑ACCCCT‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCg                                                                                                           <  1:375789/151‑1 (MQ=255)                                                aagtttttccctttttAAT‑TTTTTCCCGGGAAAATCAGTTTAAAATGCGGCCTCCTTTCGCAACCCTGTTTTGG‑TTCCCCC‑TTT‑AATAAGCC‑GG‑TTTT‑TCCCCC‑TTTTTTCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGAt                                                                                                         <  2:74882/136‑1 (MQ=255)                                                       ttCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCg                                                                                                >  1:649178/1‑150 (MQ=255)                                                            ttttAAA‑TTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTc                                                                                          >  1:368515/1‑151 (MQ=255)                                                            tttaAAT‑TTATTCCCGGTAAAATCAAGTTAAAATACGGCACCCGTCCGCAACCCTGGTTTGG‑TTCCGCA‑ATT‑AATAAGCCGGG‑TTTT‑ATCCCC‑TTTATT‑CAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTa                                                                                              <  2:994850/147‑1 (MQ=255)                                                              ttAAA‑TTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTcc                                                                                         >  1:993396/1‑150 (MQ=255)                                                                aaa‑TTATTCCCGGTAAAATCAAGTAAAAATCGGGAATCGGTCCGAAATCCTGTTTGGG‑CTTCCCA‑ATA‑AATAGGCCAGT‑TTTT‑ATCCCC‑TTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCAt                                                                                       <  2:122624/150‑1 (MQ=255)                                                                 aatttattcCCGGTAAAATCAAGTTAAAATCCGGCACCCGTCCGCAATCCTGTTTGGG‑TTTCGCC‑ATA‑AATAAGCCGGT‑TTTT‑ACCCCT‑TTTATC‑CAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCAt                                                                                       <  2:530603/147‑1 (MQ=255)                                                                  a‑ttattACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGc                                                                                    >  2:40255/1‑151 (MQ=255)                                                                     tattaCCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtg                                                                                  >  2:701109/1‑151 (MQ=255)                                                                          ccGGTAAAATAAATTTAAAATCCGGCTTCC‑TTCCCAAATCCGGTTTGGTCTCCCCC‑ATT‑AATAAGCCGGT‑TTTT‑ATCCCC‑TTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCt                                                                                   <  1:73824/145‑1 (MQ=255)                                                                                aaaTCAATTTCAATTACGGCATCCGTTCGCAATCCTGTTTGGG‑CTTCGCC‑ATT‑AATAACCCGGT‑TTTT‑ACCCCT‑TTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAgg                                                                          <  2:11122/148‑1 (MQ=255)                                                                                 aaTAAAGTAAAAATCCGGCCCCCGTCCCCAACCCTGGTTGGG‑CTCCCCC‑TTT‑AATAGCCCGGT‑TTTT‑ACCCCT‑TTTATC‑CACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATc                                                                      <  1:404619/151‑1 (MQ=255)                                                                                           aaaTCCGGCTTCCTTTCGAAATCCTGTTTGGG‑CTCCGCC‑ATA‑AAAAAGCCAGT‑TTTT‑ATCCCT‑TTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaac                                                            <  2:29611/151‑1 (MQ=255)                                                                                           aaaTCCGCCACCCGTTCGCAATCCGGGTTTGG‑CTTCCCC‑TTT‑AAAAACCCAGT‑TTTT‑ACCCCT‑TTTTTT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGaca                                                              <  2:955365/149‑1 (MQ=255)                                                                                           aaaTACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaac                                                            >  1:118383/1‑151 (MQ=255)                                                                                            attcggccctccgtccgCATCCCGGTTTTGT‑TTTCGCT‑TTA‑AATACGCCTGT‑TTTT‑ACCCCC‑TTTATC‑CCCACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGaca                                                              <  2:309764/139‑1 (MQ=255)                                                                                            aaTACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaac                                                            >  2:382733/1‑150 (MQ=255)                                                                                                  gcccccgttccgcaTCCTGGTTTGGTCTCCCCC‑ATA‑AAT‑ACCCCGG‑TTTT‑AACCCCTTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGa                                                         <  1:1001361/135‑1 (MQ=255)                                                                                                   cacccGTCCGCAACCCGGGTTGGG‑C‑TGGCC‑CTT‑TAAGAGCCCGG‑TTTTAATCACC‑TTTATC‑CACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAAc                                                      <  2:409005/146‑1 (MQ=255)                                                                                                    aTCCGTTCGCAATCCTGGTTTGG‑CTTCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACtt                                                    >  1:470706/1‑150 (MQ=255)                                                                                                       cGTTCGCAATCATGGTTTGG‑CTTCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTg                                                >  1:775950/1‑151 (MQ=255)                                                                                                         ttCGCAATCCTGGTTTGG‑CTTCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTga                                               >  2:520557/1‑150 (MQ=255)                                                                                                         ttCGCAATCCTGGTTTGG‑CTTCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTga                                               >  2:974545/1‑150 (MQ=255)                                                                                                         tcccaaccccTGTTTGGG‑CTTC‑CC‑ATTAAAAAACCCAGT‑TTTT‑ACCCCT‑TTTATC‑CACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACt                                                     <  1:701003/136‑1 (MQ=255)                                                                                                              attcttgttttGG‑CTTCCCA‑ATT‑AATACGCCGGT‑TTTT‑ATCCCT‑TTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaaga                                            <  2:320614/140‑1 (MQ=255)                                                                                                              aaTCCTGGTTTGG‑CTTCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaat                                         >  2:333567/1‑151 (MQ=255)                                                                                                                atcgggTTTGG‑GCTCCCC‑AAT‑AAAAAGCCGGG‑TTTT‑AACCCC‑TTTATTCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaat                                         <  1:55080/146‑1 (MQ=255)                                                                                                                 ccTGTTTTGG‑CTTCCCC‑ATT‑AATAAGCCAGT‑TTTT‑ATCCCC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaataat                                      <  2:18006/151‑1 (MQ=255)                                                                                                                 ccTGTTTTGG‑CTTCCCA‑ATA‑AAAAACCCAGT‑TTTT‑ACCCCC‑TTTATC‑CACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATCGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTg                                                <  1:1009220/141‑1 (MQ=255)                                                                                                                 ccTGGTTTGG‑CTCCGCC‑TTT‑AATAGCCCAGT‑TTTT‑ACCCCT‑TTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTga                                               <  1:174641/142‑1 (MQ=255)                                                                                                                      tttGG‑CTTCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtcttct                                 >  2:327663/1‑151 (MQ=255)                                                                                                                      tttGG‑CTCCGCC‑ATT‑AATAAGCCAGT‑TTTT‑ATCACC‑TTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtcttc                                  <  2:486891/150‑1 (MQ=255)                                                                                                                       ttGG‑CTTCCCCAATA‑AATAACCGAGG‑TTTA‑ATCTCT‑TTTATC‑CAAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaataa                                       <  1:561241/145‑1 (MQ=255)                                                                                                                         ttggttccgc‑ATT‑AAAAAACCCGGTTTTT‑ATCCCC‑TTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtcttct                                 <  2:451746/141‑1 (MQ=255)                                                                                                                         gg‑CTTCCCC‑ATT‑AATAAGCCAGT‑TTTT‑ACCCCC‑TTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTg                                <  2:239391/149‑1 (MQ=255)                                                                                                                         gg‑CTGCGCA‑CTT‑AAAAAGCTGGG‑TTTT‑TTCTCA‑TTTATT‑CACAACCGCCGCCCCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTg                                <  2:251296/149‑1 (MQ=255)                                                                                                                               cGCC‑ATT‑AATAAGCCAGT‑TTTT‑ATCCCT‑TTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttttt                         <  2:984389/151‑1 (MQ=255)                                                                                                                                cccattt‑AATAAGCCGGT‑TTTT‑ATCCCT‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttttt                         <  2:72370/146‑1 (MQ=255)                                                                                                                                  c‑ATT‑AATAAGCCAGG‑TTTT‑ATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTc                      >  1:599824/1‑151 (MQ=255)                                                                                                                                   ccttaaaaaaGCCGGT‑TTTT‑ATCCCC‑TTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGt                       <  2:802181/143‑1 (MQ=255)                                                                                                                                     ataaatagGCCGGT‑TTTT‑ATCCCC‑TTAATC‑AACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTc                      <  1:122539/147‑1 (MQ=255)                                                                                                                                        aataaGCCAGG‑TTTT‑ATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGCGCGATTCGCCAGCGTAACTCCATTGCTGCGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaaa                   >  1:106356/1‑150 (MQ=255)                                                                                                                                         ataaGCCAGG‑TTTT‑ATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAAt                 >  2:12824/1‑151 (MQ=255)                                                                                                                                         aaaaagcaggttttt‑ttCCCT‑TTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaaaa                  <  1:561626/140‑1 (MQ=255)                                                                                                                                           aaGCCAGT‑TTTT‑ATCCCC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATa                <  2:748929/150‑1 (MQ=255)                                                                                                                                           aaGCCAGG‑TTTT‑ATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAAt                 >  1:698282/1‑149 (MQ=255)                                                                                                                                              ccAGT‑TTTT‑ATCCCC‑TTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTg             <  1:194311/150‑1 (MQ=255)                                                                                                                                                          tCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGAttt  >  2:242995/1‑151 (MQ=255)                                                                                                                                                            |                                                                                                                                                       ATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCC‑ATT‑AATAAGCCAGG‑TTTT‑ATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTT  >  NC_000913/4537960‑4538257

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