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

CCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAAT‑A‑CGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTG  >  NC_000913/4537963‑4538253                                                                                                                                                          |                                                                                                                                                  cccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAAT‑A‑CGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAGCCAGG‑TTTT‑ATCACCtt                                                                                                                                               <  1:124226/151‑1 (MQ=255)  ccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAAT‑A‑CGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAGCCAGG‑TTTT‑ATCACCttt                                                                                                                                              <  1:916271/151‑1 (MQ=255)   cGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAAT‑A‑CGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTa                                                                                                                                             <  1:968164/151‑1 (MQ=255)    gCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAAT‑A‑CGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTAt                                                                                                                                            <  2:361753/151‑1 (MQ=255)     cAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAAT‑A‑CGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTAtt                                                                                                                                           >  2:522078/1‑151 (MQ=255)     cAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAAT‑A‑CGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTAt                                                                                                                                            <  1:384870/150‑1 (MQ=255)      aaaaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAAT‑A‑CGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑c                                                                                                                                         <  2:22453/151‑1 (MQ=255)      aaaaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAAT‑A‑CGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑c                                                                                                                                         <  1:181842/151‑1 (MQ=255)        aaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAAT‑A‑CGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑cac                                                                                                                                       >  1:635547/1‑151 (MQ=255)          ccGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAAT‑A‑CGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑CACaa                                                                                                                                     <  1:366824/151‑1 (MQ=255)                    cGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTCCCGGTGAAATCAATTTCAAAT‑A‑CGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAGCCAGT‑TTTT‑ATCCCCTTTATC‑CACAACcgccgc                                                                                                                              <  1:665890/148‑1 (MQ=255)                      ccATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAAT‑A‑CGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCAcc                                                                                                                           >  2:919099/1‑149 (MQ=255)                            ggTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAAT‑A‑CGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTc                                                                                                                   >  1:883289/1‑151 (MQ=255)                             gTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAAT‑A‑CGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTcc                                                                                                                  >  1:543025/1‑151 (MQ=255)                               gatgaGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAAT‑A‑CGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGtt                                                                                                                    >  1:11275/1‑147 (MQ=255)                                    gagggggGCAAGTTTATTCTTTTTAAAA‑TTATTCCCGTTAAAATCAATTTAAAAT‑C‑CGGCCCCCGTTCCCAACCCTGGTTTGG‑CTT‑CCCCATTAATAAGCCAGT‑TTTT‑CTCCCTTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAg                                                                                                           <  1:121557/148‑1 (MQ=255)                                        ggAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAAT‑A‑CGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCg                                                                                                       >  2:675857/1‑151 (MQ=255)                                               gCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAAT‑A‑CGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc                                                                                                >  1:848877/1‑151 (MQ=255)                                                   aTTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAAT‑A‑CGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCg                                                                                            >  2:870627/1‑151 (MQ=255)                                                    ttCCATTTAAAA‑TTATTCCCGGTAAAATAAAG‑TAAAAT‑ACCGGCATCCTTTCGCATCCCGGTTTTGC‑CTC‑CGCCATAAAAAAGCCGGG‑TTTT‑ACCCCCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGc                                                                                             <  1:861988/149‑1 (MQ=255)                                                       ccttttAAATTTATTCCCGTTAAAATCAATTTAAAAT‑C‑CGCCCCCCTTTCCCATCCTTGTTTTGG‑TTC‑CCCCATAAAAAACCCAGT‑TTTT‑ACCCCTTTTATC‑CACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGt                                                                                           <  1:722011/147‑1 (MQ=255)                                                        aTTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAAT‑A‑CGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACt                                                                                       >  2:489570/1‑151 (MQ=255)                                                         ttttAAA‑TTATTACCGGTGAAATCAAGTTCAAAT‑A‑CGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTc                                                                                      >  1:30527/1‑151 (MQ=255)                                                          tttAAA‑TTATTACCGGTGAAATCAAGTTCAAAT‑A‑CGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTcc                                                                                     >  2:806858/1‑151 (MQ=255)                                                             aaa‑TTATTCCCGGTAAAATCAAGTTAAAAT‑C‑CGCCACCCGTTCCAAATCCTGGTTTGG‑CTC‑CGCAATTAATAAGCCGGG‑TTTT‑ATCCCCTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCAtt                                                                                  <  2:262678/151‑1 (MQ=255)                                                              aa‑TTATTCCCGGTAAAATCAATTTAAAAT‑C‑CGGCATCCGTTCGCAACCCTGTTTTGG‑CTT‑CGCCATAAATAACCCAGT‑TTTT‑ATCCCCTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTg                                                                                 <  2:628157/151‑1 (MQ=255)                                                                    ttcccGGTAAAATCAGGACAAAAT‑A‑CGGCCCCCGTTCGCACCCGTGTTTGGG‑CTC‑CTCTATAAAAACGCCGGT‑TTTT‑ACCCCCTTAGCC‑CAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCt                                                                               <  1:434176/145‑1 (MQ=255)                                                                    ttcccGGTAAAATAAAGTTAAAAT‑A‑CGGCACCCGTTCGCAACCCGGTTTTGG‑CTT‑CCCCATTAATAACCCGGG‑TTTT‑ATCCCTTTTATT‑CAAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtgt                                                                             <  2:1015725/147‑1 (MQ=255)                                                                      aCCGGTGAAATCAAGTTCAAAT‑A‑CGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGtt                                                                          >  1:223186/1‑151 (MQ=255)                                                                            gAAATCAAGTTCAAAT‑A‑CGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAggg                                                                     >  2:710270/1‑150 (MQ=255)                                                                            gAAATCAAGTTCAAAT‑A‑CGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAggg                                                                     >  1:855293/1‑150 (MQ=255)                                                                            aaaaTC‑AGTTCAAATCA‑CGCCCCCCGTCCCAAATCTTGTTTTGG‑‑TT‑CCCCCATAAAAAGCCCGGTTTTT‑ACCCCCTTTATC‑CCAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAgg                                                                      <  1:516247/148‑1 (MQ=255)                                                                             aaataaattcaaaaT‑C‑CGGCATCCGTCCGCAATCCTGTTTTGG‑CTT‑CGCCATTAATAAGCCAGT‑TTTT‑ATCCCTTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGAt                                                                   <  2:899212/140‑1 (MQ=255)                                                                             aaaTCAAGTTAAAAT‑A‑CGGCACCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATAAAAAACCCAGT‑TTTT‑ACCCCCTTTATC‑CACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAggg                                                                     <  2:842152/149‑1 (MQ=255)                                                                                     ttaaaaT‑C‑CGCCATCCGTTCGCAACCCTGTTTTGG‑CTC‑CGCAATTAAAAAGCCAGT‑TTTT‑ATCCCCTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAg                                                             <  2:124052/146‑1 (MQ=255)                                                                                         aaT‑A‑CGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaac                                                        >  1:882877/1‑150 (MQ=255)                                                                                               cGGCACCCGTTTCTAATCCCTGTTTGG‑CTTACGCCACAAATACGTGAGG‑TTTT‑ATCCCTTTTATC‑AACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGaa                                                    <  1:298776/151‑1 (MQ=255)                                                                                                ggccccttttCGCCATCTTGTTTTGG‑CTT‑CGCCTAAAATAACCGTGT‑TTTT‑ATCCCCTTTACC‑CACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGaa                                                    <  1:100509/141‑1 (MQ=255)                                                                                                ggCATCCTTCCGCATTCCTGGTTTGG‑CTT‑CGCCATAAATAACCCAGT‑TTTT‑ATCCCCTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGaaa                                                   <  2:674909/150‑1 (MQ=255)                                                                                                ggCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAAc                                                  >  1:618303/1‑151 (MQ=255)                                                                                                     cctttccaaacccTGTTTTGG‑CTT‑CGCAATAAATAAGCAAGG‑TTTT‑ACCCCTTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGaaa                                                   <  2:552728/134‑1 (MQ=255)                                                                                                     ccGTTCGAAACCCGGTTTGGG‑TTT‑CCCCTTAAAAAACCCAGT‑TTTT‑ACCCCTTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTAtt                                             <  2:667437/151‑1 (MQ=255)                                                                                                        tccGCAATCC‑GGTTTGGCCTT‑CCCCATAAAAAAGCCATT‑TTTT‑ACCACTTTTACA‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTg                                            <  1:338938/147‑1 (MQ=255)                                                                                                          cgaaacccTGTTTTGG‑TTT‑CCCCATTAAAAACCCAGG‑TTTT‑ACCCCTTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaaga                                        <  2:320548/145‑1 (MQ=255)                                                                                                             aaTCCGGTTTTGG‑TTC‑CCCCATTAAAAAGCCGGT‑TTTT‑ATCCCTTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaat                                     <  2:204189/151‑1 (MQ=255)                                                                                                                ccTGGTTTGT‑CTT‑CCCCATTAAAAAGCCAGT‑TTTT‑ACCACCTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaat                                     <  2:554206/148‑1 (MQ=255)                                                                                                                  tGGTTTGG‑CTT‑CGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtct                                >  2:55005/1‑151 (MQ=255)                                                                                                                    tttggg‑CTC‑CGCCATAAAAAAGCCAGT‑TTTT‑ATCCCCTTTATC‑CAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaataat                                  <  2:347728/143‑1 (MQ=255)                                                                                                                     tttCG‑CTT‑CCCCATAAAT‑AGCCTGG‑TTTTCCCCCCCTTTATA‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCAATGCAGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaataat                                  <  1:60807/146‑1 (MQ=255)                                                                                                                       tGG‑CTT‑CGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGc                           >  1:287130/1‑151 (MQ=255)                                                                                                                            tt‑CGCCATAAAAAAGCCAGT‑TTTT‑ACCCCTTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCt                          <  2:98053/148‑1 (MQ=255)                                                                                                                             ttccccttaaaaaaGCAAGT‑TTTT‑ATCCCCTTTTTC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCttt                        <  1:602970/138‑1 (MQ=255)                                                                                                                             t‑CGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCttttt                      >  2:824710/1‑151 (MQ=255)                                                                                                                               cGCCTTTAAAAAGCCAGT‑TTTT‑ACCCCCTTTATC‑CAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCttt                        <  2:192343/148‑1 (MQ=255)                                                                                                                                  cATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTc                  >  2:368448/1‑151 (MQ=255)                                                                                                                                    ttAATAAGCCGGT‑TTTT‑ATCCCTTTTATT‑CACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCa                 <  1:675745/150‑1 (MQ=255)                                                                                                                                    ttAATAAGCAGGG‑TTTT‑ATCACTTTAATC‑CAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCa                 <  2:882507/150‑1 (MQ=255)                                                                                                                                    ttAATAACCCGGT‑TTTT‑ACCCCCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTg                    <  2:608766/147‑1 (MQ=255)                                                                                                                                      aataaGCCAGG‑TTTT‑ATCCCCTTTATC‑CACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaaaa              <  2:161890/151‑1 (MQ=255)                                                                                                                                      aataaGCCAGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaaa               >  1:700883/1‑150 (MQ=255)                                                                                                                                      aataaGCCAGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaaa               >  2:351695/1‑150 (MQ=255)                                                                                                                                         acccaggg‑TTTT‑ATCCCTTTTATC‑CAAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATa            <  1:467927/144‑1 (MQ=255)                                                                                                                                             cgagttttt‑atcccttttatc‑ccAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTg         <  1:937465/145‑1 (MQ=255)                                                                                                                                             cAGG‑TTTT‑TGCCCCTTTATTCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGt          <  1:539421/149‑1 (MQ=255)                                                                                                                                                  tttt‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACAAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAACCTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCAtt   >  1:901188/1‑151 (MQ=255)                                                                                                                                                  tttt‑ACCCCCTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCAt    <  1:1022585/150‑1 (MQ=255)                                                                                                                                                   ttt‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTg  >  1:902036/1‑151 (MQ=255)                                                                                                                                                   ttt‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTg  >  1:543904/1‑151 (MQ=255)                                                                                                                                                    tt‑ATCCCCTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTg  <  2:520926/150‑1 (MQ=255)                                                                                                                                                          |                                                                                                                                                  CCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAA‑TTATTACCGGTGAAATCAAGTTCAAAT‑A‑CGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTG  >  NC_000913/4537963‑4538253

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