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

TCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCCATTAATAAGCC‑A‑GGTTTT‑ATCACCTTT‑ATTCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTT  >  NC_000913/4537961‑4538258                                                                                                                                                         |                                                                                                                                                         tCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCCATTAATAAGCC‑A‑GGTTTT‑ATCAc                                                                                                                                                         <  1:732911/150‑1 (MQ=255)   cccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCCATTAATAAGCC‑A‑GGTTTT‑ATCACCtt                                                                                                                                                      >  1:525098/1‑151 (MQ=255)   cccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCCATTAATAAGCC‑A‑GGTTTT‑ATCACCtt                                                                                                                                                      >  2:649368/1‑151 (MQ=255)    ccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCCATTAATAAGCC‑A‑GGTTTT‑ATCACCttt                                                                                                                                                     <  2:322688/151‑1 (MQ=255)          aaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCCATTAATAAGCC‑A‑GGTTTT‑ATCACCTTT‑ATTcac                                                                                                                                              >  2:223238/1‑151 (MQ=255)           aCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCCATTAATAAGCC‑A‑GGTTTT‑ATCACCTTT‑ATTcac‑a                                                                                                                                            <  2:820786/151‑1 (MQ=255)            ccGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCCATTAATAAGCC‑A‑GGTTTT‑ATCACCTTT‑ATTCAC‑aa                                                                                                                                           <  2:744149/151‑1 (MQ=255)                cAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCCATTAATAAGCC‑A‑GGTTTT‑ATCACCTTT‑ATTCAC‑AA‑ccgc                                                                                                                                      <  2:562944/151‑1 (MQ=255)                      cGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCCATTAATAAGCC‑A‑GGTTTT‑ATCACCTTT‑ATTCAC‑AA‑CCGCCGCAcc                                                                                                                                >  1:5085/1‑151 (MQ=255)                       gCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCCATTAATAAGCC‑A‑GGTTTT‑ATCACCTTT‑ATTCAC‑AA‑CCGCCGCAcc                                                                                                                                >  2:1004383/1‑150 (MQ=255)                                     agaCGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCCATTAATAAGCC‑A‑GGTTTT‑ATCACCTTT‑ATTCAC‑AA‑CCGCCGCACCAGCCGTTCCGTAcc                                                                                                                  >  1:497613/1‑150 (MQ=255)                                     agaCGGAAGAAACCTTTTTCCGTTTAAAATTATGCCGGGTGAAATCAAGTCAAAATCCGCCACCCTTGCAC‑ATCCGGGTTTGGGCTCCGCCATAAATAGCCC‑G‑GTTTTT‑ATCCCCTTT‑ATTCAC‑AA‑CCGCCGCACCAGCCGTTcc                                                                                                                       <  2:285182/145‑1 (MQ=255)                                      gaggggAGCAGGCTTATCCCATTTTAAATTATTCCCGGTAAAATAAATTTAAAATACGGCATCCGTCCGAAATCCTGTTTT‑GGCTTCGCCATTAATAACCA‑A‑GGTTTT‑ATCCCCTTT‑ATTCAC‑AA‑CCGCCGCACCAGCCGTTCCGTAc                                                                                                                   <  1:488162/145‑1 (MQ=255)                                        cGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCCATTAATAAGCC‑A‑GGTTTT‑ATCACCTTT‑ATTCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAggg                                                                                                              >  1:978341/1‑151 (MQ=255)                                         gggAGCAAGCTTATCCCATTTAAAATTATTCCCGGTAAAATCAAGTTAAAATCCGGCATCCTTTCGCAATCCTGTTTT‑GGCTTCCCCATAAATAAGCC‑A‑GGTTTT‑ATCCCCTTT‑TTCCAA‑AC‑CCGCCGCACCAGCCGTTCCGTACCAgg                                                                                                               <  1:593358/149‑1 (MQ=255)                                           gAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCCATTAATAAGCC‑A‑GGTTTT‑ATCACCTTT‑ATTCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCg                                                                                                            >  1:371981/1‑150 (MQ=255)                                               aaGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCCATTAATAAGCC‑A‑GGTTTT‑ATCACCTTT‑ATTCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCg                                                                                                       >  2:791915/1‑151 (MQ=255)                                                   tttttCCATTTTAAATTATTCCCGGTAAAACAAATTTCAAATCCGGCACCCGTTCGCAATCCTGTTTT‑GGCTTCGCCATTAATAAGCC‑A‑GTTTTT‑ATCCCCTTT‑ATCCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc                                                                                                     <  1:838071/146‑1 (MQ=255)                                                   ttATTCCATTTAAAATTTTTCCGGGGGAAATCAAGTTAAAATCCGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCCATAAATAACCC‑A‑GGTTTT‑ACCCCCTTT‑ATCCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc                                                                                                     <  1:679447/149‑1 (MQ=255)                                                    tATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCCATTAATAAGCC‑A‑GGTTTT‑ATCACCTTT‑ATTCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGc                                                                                                  >  1:18652/1‑151 (MQ=255)                                                    tATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCCATTAATAAGCC‑A‑GGTTTT‑ATCACCTTT‑ATTCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGc                                                                                                  >  2:723132/1‑151 (MQ=255)                                                      ttCTATTTAAAATTATTCCCGGTAAAATAAATTTCAAAACCGCCATCCTTCCGCAACCCTGTTTG‑GGCTCCCCCATAAAAAAGCA‑G‑GGTTTT‑ACCCCTTTT‑ATCCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGc                                                                                                  <  2:183357/149‑1 (MQ=255)                                                      ttCCATTTTAAATTATTCCCGGTAAAATCAATTTCAAATACGGCACCCGTTCGCAATCCTGTTTT‑GGTTTCGCCATTAATAAGCC‑A‑GGTTTT‑ATCCCCTTT‑ATTCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGc                                                                                                  <  1:42994/149‑1 (MQ=255)                                                      ttCCATTTAAAATTATTCCCGGTAAAATAAAGTTAAATTCCGCCATCCGTTCGCAATCCGGGTTG‑GGCTTCGCCATTAAAAACCC‑G‑GGTTTA‑ATCCCTTTT‑ATTCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCg                                                                                                 <  1:1019206/150‑1 (MQ=255)                                                        ccATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCCATTAATAAGCC‑A‑GGTTTT‑ATCACCTTT‑ATTCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTa                                                                                               >  2:781805/1‑150 (MQ=255)                                                         cATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCCATTAATAAGCC‑A‑GGTTTT‑ATCACCTTT‑ATTCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAAc                                                                                             >  1:658266/1‑151 (MQ=255)                                                           ttttAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCCATTAATAAGCC‑A‑GGTTTT‑ATCACCTTT‑ATTCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTc                                                                                           >  1:523332/1‑151 (MQ=255)                                                            tttAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCCATTAATAAGCC‑A‑GGTTTT‑ATCACCTTT‑ATTCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTcc                                                                                          >  1:19752/1‑151 (MQ=255)                                                              tAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCCATTAATAAGCC‑A‑GGTTTT‑ATCACCTTT‑ATTCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCAt                                                                                        >  1:69386/1‑151 (MQ=255)                                                               aaaTTATTCCCGGTGAAATCAATTACAAATCCGCCTCCCGTCGGCAACCCTGTTTT‑GGTTCCCCTATTAA‑AAGCCAA‑GGTTTT‑ATCACCTTT‑TTGCCA‑AT‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCAt                                                                                        <  1:932028/150‑1 (MQ=255)                                                               aaaTTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCCATTAATAAGCC‑A‑GGTTTT‑ATCACCTTT‑ATTCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCAtt                                                                                       >  2:520980/1‑151 (MQ=255)                                                                     ttaCCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCCATTAATAAGCC‑A‑GGTTTT‑ATCACCTTT‑ATTCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtgtg                                                                                 >  1:888432/1‑151 (MQ=255)                                                                     ttaCCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCCATTAATAAGCC‑A‑GGTTTT‑ATCACCTTT‑ATTCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtgtg                                                                                 >  2:63698/1‑151 (MQ=255)                                                                        ccGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCATCCCTGTTTG‑GGTTTCGCCATTAATAAGCC‑A‑GTTTTT‑TTCACCTTT‑ATTCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTg                                                                              <  1:733049/151‑1 (MQ=255)                                                                        ccGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCCATTAATAAGCC‑A‑GGTTTT‑ATCACCTTT‑ATTCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTg                                                                              >  1:280870/1‑151 (MQ=255)                                                                         cGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCCATTAATAAGCC‑A‑GGTTTT‑ATCACCTTT‑ATTCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTg                                                                              >  1:350242/1‑150 (MQ=255)                                                                           gTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCCATTAATAAGCC‑A‑GGTTTT‑ATCACCTTT‑ATTCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAgg                                                                           >  2:416225/1‑151 (MQ=255)                                                                              aaaTCAATTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCCCCATTAATAAGCC‑G‑GGTTTT‑ATCCCCTTT‑ATTCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGa                                                                         <  1:314375/150‑1 (MQ=255)                                                                              aaaTCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCCATTAATAAGCC‑A‑GGTTTT‑ATCACCTTT‑ATTCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGAt                                                                        >  1:51010/1‑151 (MQ=255)                                                                              aaaCCAATTTCAAACACGCCCTCCGTTCCCAATCCTGGTTG‑GGCTCCCCCATAAATAAGCC‑‑‑TGTTTT‑AACCCCTTTAATCCCC‑AC‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACGCCATTGCTGTGTTGAGGGAt                                                                        <  2:192640/151‑1 (MQ=255)                                                                                  cAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCCATTAATAAGCC‑A‑GGTTTT‑ATCACCTTT‑ATTCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaa                                                                     >  1:940542/1‑150 (MQ=255)                                                                                   aattataaaaACGGCATCCCTCCGCAAACCTGGTGT‑GGCTTCCCCATTAAAAAGCC‑A‑TTTTTT‑ACCCCTTTT‑AATCAC‑AA‑CCGCCCCCCCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaaa                                                                    <  1:868251/144‑1 (MQ=255)                                                                                   aaGTTAAAATCCGCCACCCTTCCGCATTCCTGTTTT‑GGCTCCCCCA‑TAAAAAACC‑C‑GGTTTT‑TACCCCTTT‑TTTCCCAAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaaa                                                                    <  2:398546/150‑1 (MQ=255)                                                                                    attaaaaaTCCGGCCCCCGTCCACAGTCCGGGTTT‑GCCTCCGCAATAAAAAAGCC‑A‑GTTTTT‑ACCCCCTTT‑ATTAAA‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaaaa                                                                   <  2:789108/145‑1 (MQ=255)                                                                                      ttCAAATACGGCATCCGTTCGCAATCCTGTTTT‑GGCTTCCCCATTAATAAGCC‑A‑GGTTTT‑ATCCCCTTT‑ATTCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAg                                                                  <  1:898057/149‑1 (MQ=255)                                                                                         aaaTCCGGCACCCGTCCGCAATCCTGTTTT‑GGCTTCGCCATTAATAACCC‑A‑GGTTTT‑ATCCCTTTT‑ATTCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGac                                                                <  1:351616/148‑1 (MQ=255)                                                                                         aaaTCCGCCACCCGTTCGCAACCCGGTTTT‑GGCTTCGCAATTAATAACCC‑A‑GTTTTT‑ATCCCCTTT‑ATTCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGaca                                                               <  1:109389/149‑1 (MQ=255)                                                                                         aaaTACGGCTTCCTTTCGCATTCTGGTTTT‑GGCTTCCCCATAAAAAAGCC‑G‑GGTTTT‑ACCCCTTTT‑ATCCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaac                                                             <  1:387870/151‑1 (MQ=255)                                                                                         aaaTACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCCATTAATAAGCC‑A‑GGTTTT‑ATCACCTTT‑ATTCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaac                                                             >  1:685390/1‑151 (MQ=255)                                                                                         aaaTACGGCATCCGTTCGAAATCCGGTTTT‑GGCTTCGCCATTAATAAGCC‑A‑GTTTTT‑ATCCCCTTT‑ATTCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGaca                                                               <  1:399430/149‑1 (MQ=255)                                                                                         aaaTACGGCACCCGTCCGCAATCCTGTTTG‑GCCTTCGCCATAAATAAGCC‑A‑GTTTTT‑ATCCCCTTT‑ATTCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaa                                                                     <  2:525861/143‑1 (MQ=255)                                                                                            tACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCCATTAATAAGCC‑A‑GGTTTT‑ATCACCTTT‑ATTCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaac                                                             >  2:737100/1‑148 (MQ=255)                                                                                              cGGCATCCGTTCGCAATCCTGGTTT‑GGCTCCGCCATTAATAAGCC‑A‑GTTTTT‑ATCCCCTTT‑ATCCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGa                                                          <  1:1017784/149‑1 (MQ=255)                                                                                                 ccccctttCGCATCCCTGTTTT‑GGTTCCCCCATAAATAAGCC‑G‑GTTTTT‑ATCCCTTTT‑ATCCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGa                                                          <  2:935391/140‑1 (MQ=255)                                                                                                    ccGTTCGCAATCTGGTTTT‑GGCTTCGCCATTAAAAACCC‑A‑GGTTTT‑ATCCCCTTT‑ATTCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTa                                                    <  1:158824/149‑1 (MQ=255)                                                                                                    ccGTTCGCAATCCTGGTTT‑GGCTTCGCCATTAATAAGCC‑A‑GGTTTT‑ATCACCTTT‑ATTCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTAtt                                                  >  2:402847/1‑151 (MQ=255)                                                                                                    ccGTTCGCAATCCTGGTTT‑GGCTTCCCCATTAAAAACCC‑G‑GGTTTT‑ACCCCTTTT‑ATTCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTa                                                    <  1:27168/149‑1 (MQ=255)                                                                                                      tttCACAATCATGTTTT‑GGTTTCG‑TATTAAAAAACC‑AGGGTTTT‑ATC‑CCTTT‑AATCCC‑AATCCGCCGCGCCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTa                                                    <  1:878568/146‑1 (MQ=255)                                                                                                      gTTCGCAATCCTGGTTT‑GGCTTCGCCATTAATAAGCC‑A‑GGTTTT‑ATCACCTTT‑ATTCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTg                                                 >  1:62085/1‑150 (MQ=255)                                                                                                           cAACCCGGGTTG‑GGCTTCGCTATTAAAAACCG‑A‑GGTTTC‑ACCCCCTTT‑TTCGAA‑AC‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAGCAGAAACtt                                                     <  2:603083/141‑1 (MQ=255)                                                                                                            aaTCCTGGTTT‑GCTTTCCCTATTAATAACCC‑G‑GGTTTT‑ATCCCCTTT‑ATTCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaaga                                             <  1:1004414/148‑1 (MQ=255)                                                                                                               ccTGGTTT‑GGCTTCGCCATTAATAAGCC‑A‑GGTTTT‑ATCACCTTT‑ATTCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaataat                                       >  2:861035/1‑151 (MQ=255)                                                                                                               ccTGGTTG‑GGCTTCGCAATAAATAAGCC‑A‑GGTTTT‑ACCCCCTTT‑ATTCAA‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaataat                                       <  1:1063080/151‑1 (MQ=255)                                                                                                               ccTGGTCA‑GCCTTCACCATAAATAAGCC‑A‑GTTTTT‑ACCCCCTTA‑TTACAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTa                                                    <  2:290258/138‑1 (MQ=255)                                                                                                                  gtttt‑GGCTCCCCCATAAAAAAGCC‑G‑GGTTTT‑ATCCCCTTA‑ATCCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaataat                                       <  2:824672/146‑1 (MQ=255)                                                                                                                  ggTTT‑GGCTCCCCAATTAAAAGGCC‑A‑GTTTTT‑ATCCCTTTT‑ATCCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaataat                                       <  2:639759/148‑1 (MQ=255)                                                                                                                    ttt‑GGCTCCGCCATTAATAAGCC‑G‑GGTTTT‑ATCCCCTTT‑ATCCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtc                                      <  1:1048642/147‑1 (MQ=255)                                                                                                                     ttggtctccGCCATAAAAAACCC‑A‑GTTTTT‑TCCCCTTTT‑ATTCCA‑AT‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaataa                                        <  2:743677/137‑1 (MQ=255)                                                                                                                        ggctccccaaTTAAT‑AGCC‑C‑GGTTTTCACCCCTTTT‑ATTCCC‑AC‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtcttct                                  <  2:633823/139‑1 (MQ=255)                                                                                                                        ggCTTCGCCATTAATAGGCC‑A‑GGTTTA‑ATCCCTTTT‑ATTCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTg                                 <  1:655965/149‑1 (MQ=255)                                                                                                                        ggCTCCGCCATTAATAACCC‑A‑GTTTTT‑ATCCCCTTT‑ATTCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGc                                <  1:655521/150‑1 (MQ=255)                                                                                                                        gccTTCCCCATTAAAAAGCC‑A‑GTTTTT‑ACCCCTTTT‑ATCCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTg                                 <  2:801127/147‑1 (MQ=255)                                                                                                                          cTTCGCCATTAATAAGCC‑A‑GGTTTT‑ATCACCTTT‑ATTCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACCTATTGAAGAAATAATCTTCTGCttt                             >  2:232918/1‑151 (MQ=255)                                                                                                                          cTTCGCCATAAATAAGCC‑A‑GTTTTT‑ATCCCTTTT‑ATCCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGc                                <  2:940126/148‑1 (MQ=255)                                                                                                                           ttCGCCATTAATAAGCC‑A‑GGTTTT‑ATCACCTTT‑ATTCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttt                            >  2:15694/1‑151 (MQ=255)                                                                                                                            tCGCCATTAATAAGCC‑A‑GGTTTT‑ATCACCTTT‑ATTCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCttttt                           >  2:192670/1‑151 (MQ=255)                                                                                                                             ccccATAAAAAGGCC‑A‑TTTTTT‑ACCCCTTTT‑ATCCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttt                            <  2:946271/147‑1 (MQ=255)                                                                                                                             ccccATAAAAAAGCC‑A‑GTTTTT‑ATCCCTTTT‑ATCCAA‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCTATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttttt                          <  1:649389/149‑1 (MQ=255)                                                                                                                              gCCATTAATAAGCC‑A‑GGTTTT‑ATCACCTTT‑ATTCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTg                         >  1:814881/1‑151 (MQ=255)                                                                                                                               cattaaataaGCC‑A‑GGTTTT‑ATCTCCTTT‑ATTCCC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGt                        <  1:801063/146‑1 (MQ=255)                                                                                                                                  ttAATAAGCC‑A‑GTTTTT‑ATCCCTTTT‑ATCCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCa                      <  1:956556/150‑1 (MQ=255)                                                                                                                                   taataaGCC‑A‑GGTTTT‑ATCACCTTT‑ATTCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaa                     >  2:288347/1‑150 (MQ=255)                                                                                                                                     ataaGCC‑A‑GTTTTT‑ATCCCCTTT‑ATTCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAAt                  <  1:892634/151‑1 (MQ=255)                                                                                                                                       aaccc‑A‑GTTTTT‑ATCCCCTTT‑ATCCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATa                 <  1:70680/147‑1 (MQ=255)                                                                                                                                          cc‑A‑GGTTTT‑ATCCCCTTT‑ATTCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTg              <  1:595380/150‑1 (MQ=255)                                                                                                                                               ggTTTT‑ATCACCTTT‑ATTCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGc           >  1:310416/1‑150 (MQ=255)                                                                                                                                                 tttt‑ATCACCTTT‑ATTCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCAtt        >  2:850517/1‑151 (MQ=255)                                                                                                                                                      gtCATCTTT‑ATTCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGTCAGCGTAACTCCATTGCGGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCTAAATAGTGAGCATTGAtt    >  2:341528/2‑151 (MQ=255)                                                                                                                                                      aTCACCTTT‑ATTCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGAtt    >  1:110238/1‑151 (MQ=255)                                                                                                                                                        cACCTTT‑ATTCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGAtttt  >  2:862772/1‑151 (MQ=255)                                                                                                                                                         |                                                                                                                                                         TCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCCATTAATAAGCC‑A‑GGTTTT‑ATCACCTTT‑ATTCAC‑AA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTT  >  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 < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑