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

ATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCAT‑CCGTTCGCAATCCTGGTTT‑GG‑C‑TTCGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGAT  >  NC_000913/4537960‑4538255                                                                                                                                                            |                                                                                                                                                    aTCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCAT‑CCGTTCGCAATCCTGGTTT‑GG‑C‑TTCGCCATTAATAAGCCAGG‑TTTT‑ATCa                                                                                                                                                     <  2:43018/150‑1 (MQ=255)  tCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCAT‑CCGTTCGCAATCCTGGTTT‑GG‑C‑TTCGCCATTAATAAGCCAGG‑TTTT‑ATCAcc                                                                                                                                                   >  2:226286/1‑151 (MQ=255)        cAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCAT‑CCGTTCGCAATCCTGGTTT‑GG‑C‑TTCGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTAtt                                                                                                                                             <  2:741880/151‑1 (MQ=255)         aaaaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCAT‑CCGTTCGCAATCCTGGTTT‑GG‑C‑TTCGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑c                                                                                                                                           >  2:774042/1‑151 (MQ=255)          aaaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCAT‑CCGTTCGCAATCCTGGTTT‑GG‑C‑TTCGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑c                                                                                                                                           >  2:869913/1‑150 (MQ=255)            aCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCAT‑CCGTTCGCAATCCTGGTTT‑GG‑C‑TTCGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑caca                                                                                                                                        <  1:890281/151‑1 (MQ=255)                ccAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCAT‑CCGTTCGCAATCCTGGTTT‑GG‑C‑TTCGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑CACAAccg                                                                                                                                    >  2:964902/1‑151 (MQ=255)                     tACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCAT‑CCGTTCGCAATCCTGGTTT‑GG‑C‑TTCGCCATTAATAAGCCAGG‑TTTT‑ATCCCCTTTATT‑CACAAccgccg                                                                                                                                 <  1:468592/149‑1 (MQ=255)                     tACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCAT‑CCGTTCGCAATCCTGGTTT‑GG‑C‑TTCGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑CACAACcgccgc                                                                                                                                >  1:389216/1‑150 (MQ=255)                         ccATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCAT‑CCGTTCGCAATCCTGGTTT‑GG‑C‑TTCGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCAcc                                                                                                                             >  2:874003/1‑149 (MQ=255)                           aTCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCAT‑CCGTTCGCAATCCTGGTTT‑GG‑C‑TTCGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCACCAGcc                                                                                                                         >  2:248911/1‑151 (MQ=255)                             cTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCAT‑CCGTTCGCAATCCTGGTTT‑GG‑C‑TTCGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGt                                                                                                                       >  2:845708/1‑151 (MQ=255)                               ggggatgaGACGGGAGCAAGCTTTTTCCATTTTAATTTTTTACCGGTAAAATCAATTTCAAATACGGCTC‑CCGTCCGCAATCCGGGTTG‑GG‑C‑TTCGCCATTAATAAGCCAGT‑TTTT‑ACCCCCTTTATC‑CACAACCGCCGCACCAGcc                                                                                                                         <  1:1062275/144‑1 (MQ=255)                                       gaCGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCAT‑CCGTTCGCAATCCTGGTTT‑GG‑C‑TTCGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTAcc                                                                                                               >  2:59515/1‑149 (MQ=255)                                          gggAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCAT‑CCGTTCGCAATCCTGGTTT‑GG‑C‑TTCGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGc                                                                                                          >  1:427881/1‑151 (MQ=255)                                          gggAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCAT‑CCGTTCGCAATCCTGGTTT‑GG‑C‑TTCGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTACT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGc                                                                                                          >  2:2290/1‑151 (MQ=255)                                          gggAGCAAGCTTATTCCATTTAAAATTATTCCCGGTAAAATCAATTCAAAATCCGGCAC‑CCTTTCCCATCCCTGTTTG‑GC‑C‑TTCCCCATAAAAAACCCAGG‑TTTT‑ACCCCCTTTATC‑CACACCCGCCGCACCAGCCGTTCCGTACCAGGGc                                                                                                          <  1:338798/151‑1 (MQ=255)                                          gggAGAAAGTTTATCCCTTTTAAAATTATTCCCGGTAAAACCAATTTAAAATCCGGCAC‑CCTTCCGCAACCCTGGTTG‑GG‑T‑TCCCCCATTAAAAAGCCAGT‑TTTT‑TCCCCCTTTTTA‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGc                                                                                                          <  1:591948/151‑1 (MQ=255)                                          gggAGAAAGCTTTTTCATTTTAAAATTATTCCCGGTAAAACCAATTTAAAATACGGCAT‑CCGTTCGCAACCCTGTTTT‑GG‑C‑TTCGCCATTAATAAGCCAGT‑TTTT‑ACCCCCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAgg                                                                                                            <  2:976789/149‑1 (MQ=255)                                               cAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCAT‑CCGTTCGCAATCCTGGTTT‑GG‑C‑TTCGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTc                                                                                                     >  2:703991/1‑151 (MQ=255)                                                aggCTTCTTCAATTTAAAATTATTCCCGGTAAAATAAATTGCAAATACGGTTG‑CTGTTCGAAACCCTGTTGT‑GG‑C‑TTCGCCATGAA‑AACGCAGG‑TTTTCATCCCTTTTAGA‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTc                                                                                                     <  2:1066216/148‑1 (MQ=255)                                                aaGTTTTTTCCTTTTAAAATTATTCCGGGTAAAACAAATTAAAAATACGCCAC‑CCGTCCGAAACCCTGTTTT‑GG‑C‑TTCGCCTTTAATAAGCCAGG‑TTTT‑ATCCCTTTTATT‑CAAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGAtt                                                                                                      <  1:795901/149‑1 (MQ=255)                                                aaGTTTATTCCATTTTAAATTATTCCCGGTAAAATCAAGTTCAAATACGGCAC‑CCTTTCGCAATCCTGTTTT‑GG‑C‑TTCGCAATTAATAAGCCAGG‑TTTT‑ATCCCTTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTc                                                                                                     <  2:411511/150‑1 (MQ=255)                                                      aTTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCAT‑CCGTTCGCAATCCTGGTTT‑GG‑C‑TTCGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCg                                                                                              >  1:826990/1‑151 (MQ=255)                                                       ttCCATTTTAAATTATTCCCGGGAAAATCAAGTTAAAATACGGCAT‑CCGTTCGCAATCCGGGTTT‑GG‑T‑TTCGCCATTAAAAAGCCAGT‑TTTT‑ATCCCTTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGc                                                                                               <  2:951806/149‑1 (MQ=255)                                                       ttCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCAT‑CCGTTCGCAATCCTGGTTT‑GG‑C‑TTCGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGt                                                                                             >  2:256963/1‑151 (MQ=255)                                                       ttCCATTTAAAATTATTCCCGGTGAAATCAAGTTCAAATACGGCAT‑CCGTTCGCAATCCTGTTTT‑GG‑C‑TTCGCCATTAATAAGCCAGG‑TTTT‑ATCCCCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCg                                                                                              <  2:780022/150‑1 (MQ=255)                                                          cATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCAT‑CCGTTCGCAATCCTGGTTT‑GG‑C‑TTCGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAAc                                                                                          >  2:194695/1‑151 (MQ=255)                                                            tttaAAATTATTCCCGGTAAAATAAATTTAAAATCCGGCAC‑CCTTTCGAAATCCTGTTTT‑GG‑T‑TTCCCCATTAATAAGCCAGG‑TTTT‑ATCCCCTTTATC‑CAAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTaa                                                                                           <  2:1114506/148‑1 (MQ=255)                                                             tttAAATTATTACCGGTGAAATCAAGTTCAAATACGGCAT‑CCGTTCGCAATCCTGGTTT‑GG‑C‑TTCGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTcc                                                                                       >  1:914365/1‑151 (MQ=255)                                                                aaaTTATCCCCGGAAAAACAATGTTCAAAT‑CCGCATCCCGTTC‑CAATCCTGGTTT‑GGCC‑TTCCCCTTAAATACCCGGGT‑TTTT‑CTCCCTTTTTTC‑CAAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCAtt                                                                                    <  1:378279/151‑1 (MQ=255)                                                                     attcCCGGTAAAATCAATTACAAATACGGCAT‑CCGTTCGCAATCCTGTTTT‑GG‑C‑TTCCCCATAAATAACCCGGT‑TTTT‑ACCCCCTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtgt                                                                               <  2:673035/151‑1 (MQ=255)                                                                      ttcccGGTAAAATCAAGTTCAAATCCGGCAT‑CCGTCCGCAATCCTGTTTT‑GG‑C‑TTCGCCATTAAAAAGCCAGT‑TTTT‑ATCCCCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtgt                                                                               <  2:167996/147‑1 (MQ=255)                                                                       taCCGGTGAAATCAAGTTCAAATACGGCAT‑CCGTTCGCAATCCTGGTTT‑GG‑C‑TTCGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTgtgt                                                                             >  1:620044/1‑151 (MQ=255)                                                                       taCCGGTAAAATCAAGTTAAAATACGGCAT‑CCGTTCGAAATCCTGTTTT‑GG‑C‑TTCGCCATTAATAAGCCAGG‑TTTT‑ACCCCCTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtgtg                                                                              <  2:680275/150‑1 (MQ=255)                                                                         ccGGTAAAATCAAGTTCAAATACGGCAT‑CCGTTCGCAACCCTGGTTT‑GG‑T‑TCCGCCATTAATAAGCCAGT‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTCCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtgt                                                                               <  1:934646/147‑1 (MQ=255)                                                                         ccGGTAAAATAAATTTAAAATACGGCAT‑CCGTTCGCAATCCTGTTTG‑GG‑C‑TTCGCCATAAATAAGCCTGT‑TTTT‑ATCCCCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGtt                                                                            <  2:511776/150‑1 (MQ=255)                                                                          cGGTGAAATCAAGTTCAAATACGGCAT‑CCGTTCGCAATCCTGGTTT‑GG‑C‑TTCGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTg                                                                           >  1:491505/1‑150 (MQ=255)                                                                           ggTGAAATCAAGTTCAAATCCGGCAT‑CCTTTCGCAACCCTGGTTG‑GG‑C‑TTCCCCATTAATAACCCGGT‑TTTT‑ATCCCCTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAg                                                                         <  2:347966/151‑1 (MQ=255)                                                                           ggTAAAATCATGTTAAAATGCGGC‑C‑CCGGTCCCAATTCCGGTTT‑TG‑CTTTCGCCATTAAAAACCCAGT‑TTTT‑CCCCCTTTTTTC‑CCCAGCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTg                                                                           <  1:867373/149‑1 (MQ=255)                                                                              gAAATCAAGTTCAAATACGGCAT‑CCGTTCGCAATCCTGGTTT‑GG‑C‑TTCGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAggg                                                                       >  2:332057/1‑150 (MQ=255)                                                                               aaaTCAAGTTCAAATACGGCAT‑CCGTTCGCAATCCTGGTTT‑GG‑C‑TTCGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGAt                                                                     >  1:1012261/1‑151 (MQ=255)                                                                                aaTCAAGTTAAATTCCGGCAC‑CCTTCCGAAATCCGGTTTT‑GG‑C‑TCCGCCATTAATAAGCCGGG‑TTTT‑ACCCCCTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATc                                                                    <  2:946681/151‑1 (MQ=255)                                                                                  tCAAGTTCAAATACGGCAT‑CCGTTCGCAATCCTGGTTT‑GG‑C‑TTCGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCa                                                                   >  2:162647/1‑150 (MQ=255)                                                                                   cAAGTTCAAATACGGCAT‑CCGTTCGCAATCCTGGTTT‑GG‑C‑TTCGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaa                                                                  >  2:713338/1‑150 (MQ=255)                                                                                    aaGTTCAAATACGGCAT‑CCGTTCGCAATCCTGGTTT‑GG‑C‑TTCGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaaa                                                                 >  1:403917/1‑150 (MQ=255)                                                                                    aaGGTCAAATACGACAT‑CCGTTCGCAATCCTGGTTT‑GG‑C‑TTCGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaaa                                                                 >  1:608372/1‑150 (MQ=255)                                                                                      gTTCAAATACGGCAT‑CCGTTCGCAATCCTGGTTT‑GG‑C‑TTCGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaaa                                                                 >  2:856078/1‑148 (MQ=255)                                                                                       ttCAAATACGGCAT‑CCGTTCGCAATCCTGGTTT‑GG‑C‑TTCGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGac                                                             >  2:623927/1‑151 (MQ=255)                                                                                          aaaTCCGGCAC‑CCGTCCGCAATCCTGGTTT‑GG‑T‑TTCGCAATTAAAAAGCCAGT‑TTTT‑ACCCCCTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGaca                                                            <  2:494131/149‑1 (MQ=255)                                                                                          aaaTACGGCCT‑CCGTTCGCAATCCTGTTTT‑GG‑T‑TTCGCCATTAATAAGCCAGT‑TTTT‑ATCCCCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaa                                                           <  2:538059/150‑1 (MQ=255)                                                                                          aaaTACGGCAT‑CCGTTCGCAATCCTGGTTT‑GG‑C‑TTCGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaac                                                          >  1:172452/1‑151 (MQ=255)                                                                                           aaTCCGGCAC‑CCGTTCGCAATCCTGTTTT‑GG‑C‑TTCCCCATAAAAAGGCCAGT‑TTTT‑ACCCCCTTTATC‑CCCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaaca                                                         <  1:586870/151‑1 (MQ=255)                                                                                              aCGGCAT‑CCGTTCGCAATCCTGGTTT‑GG‑C‑TTCGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGa                                                       >  2:1089168/1‑150 (MQ=255)                                                                                                  cAT‑CCGTTCGCAATCCTGGTTT‑GG‑C‑TTCGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACtt                                                  >  2:555941/1‑151 (MQ=255)                                                                                                      cctttCGCACCCCTGTTTT‑GG‑C‑TTCCCAATAAATAAGCCAGT‑TTTT‑ATCCCTTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACt                                                   <  2:33238/144‑1 (MQ=255)                                                                                                      cctttCGCAATCCTGGTTT‑GG‑T‑TTCCCCATTAAAAAGCCGGT‑TTTT‑ACCCCTTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTAt                                                <  1:147162/147‑1 (MQ=255)                                                                                                        gTTCGCAATCCTGGTTT‑GG‑C‑TTCGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTg                                              >  1:506596/1‑150 (MQ=255)                                                                                                         ttCGCAATCCTGTTTT‑GG‑T‑TTCGCCATTAAAAAGCCGGT‑TTTT‑ATCCCCTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTga                                             <  2:714713/150‑1 (MQ=255)                                                                                                         ttCGC‑ATCCTGGTTTGGG‑C‑TTCCCCATAAAAAAGCCAGT‑TTTT‑ATCCCTTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTg                                              <  1:912144/149‑1 (MQ=255)                                                                                                                 ccTGTTTT‑GG‑C‑TTCCCCATTAATAAGCCAGT‑TTTT‑ATCCCCTTTATT‑CCCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaata                                      <  2:24206/149‑1 (MQ=255)                                                                                                                 ccTGTTTT‑GC‑C‑TCCCCAATAAATAGGCCGGT‑TTTT‑ACCCCCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGaaa                                        <  1:1072838/147‑1 (MQ=255)                                                                                                                      ttt‑GG‑T‑TTCGCCTTTATTAAGCCAGT‑TTTT‑ATCCCCTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtc                                   <  2:768477/147‑1 (MQ=255)                                                                                                                          gg‑C‑TTCGCCATTAAAAACCCAGT‑TTTT‑ATCCCTTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTg                              <  2:730407/149‑1 (MQ=255)                                                                                                                                 cGCCATTAAAAAGCCAGG‑TTTT‑ACCCCTTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttt                         <  2:294279/149‑1 (MQ=255)                                                                                                                                  gCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTg                      >  2:535428/1‑151 (MQ=255)                                                                                                                                  cccATTAAAAAGCCCGGTTTTT‑ATCCCTTTTTTC‑CAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttttt                       <  1:975899/150‑1 (MQ=255)                                                                                                                                  cccAATAAAAAGGCAGG‑TTTT‑TACCCCTTTTTTCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttttt                       <  1:39512/150‑1 (MQ=255)                                                                                                                                   ccATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGt                     >  2:737794/1‑151 (MQ=255)                                                                                                                                      ttAAAAAGCCAGT‑TTTT‑ATCCCTTTTATC‑CACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGt                     <  2:263163/148‑1 (MQ=255)                                                                                                                                         ataaGCCAGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaaaa                >  1:649126/1‑150 (MQ=255)                                                                                                                                         ataaGCCAGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAAt               >  1:398711/1‑151 (MQ=255)                                                                                                                                         aaaacccAGT‑TTTT‑ACCACTTTTTGC‑CACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAAt               <  1:541432/146‑1 (MQ=255)                                                                                                                                           aaGCAAGT‑TTTT‑ATCCCTTTAATT‑CACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATa              <  1:994566/150‑1 (MQ=255)                                                                                                                                             gCCAGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTg           >  2:782469/1‑151 (MQ=255)                                                                                                                                                aGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGc        >  1:183231/1‑151 (MQ=255)                                                                                                                                                    tttt‑ATCCCCTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCAt      <  2:566900/150‑1 (MQ=255)                                                                                                                                                       t‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGAt  >  1:1070038/1‑151 (MQ=255)                                                                                                                                                            |                                                                                                                                                    ATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCAT‑CCGTTCGCAATCCTGGTTT‑GG‑C‑TTCGCCATTAATAAGCCAGG‑TTTT‑ATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGAT  >  NC_000913/4537960‑4538255

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