Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,228,222	0	T	G	43.5%	110.7 / 13.0	69	pseudogene (388/1287 nt)	mdtQ	pseudogene; putative channel/filament proteins
Reads supporting (aligned to +/- strand):  ref base T (10/29);  new base G (30/0);  total (40/29)
Fisher's exact test for biased strand distribution p-value = 3.75e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGC‑GTTCCGCCGCCCGTGC‑TTTAACCGTA‑CCCAGGCGGGCAGT‑AACCT‑CCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCC  >  NC_000913/2228072‑2228367                                                                                                                                                        |                                                                                                                                                       aagaagTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGt                                                                                                                                                        <  2:917780/151‑1 (MQ=255)  agaagTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGtt                                                                                                                                                       >  1:933068/1‑151 (MQ=255)  agaagTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGtt                                                                                                                                                       <  1:914972/151‑1 (MQ=255)    aagTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTg                                                                                                                                                     <  2:43477/151‑1 (MQ=255)    aagTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACGCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGgtt                                                                                                                                                      <  2:916240/150‑1 (MQ=255)      gTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCt                                                                                                                                                   <  2:840720/151‑1 (MQ=255)        aaTCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGttt                                                                                                                                                      >  1:445201/1‑146 (MQ=255)         aTCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGgtt                                                                                                                                                      >  2:915419/1‑145 (MQ=255)         aTCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTcgc                                                                                                                                                <  2:121688/151‑1 (MQ=255)          tCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGGTTTCTcgc                                                                                                                                                >  2:919013/1‑150 (MQ=255)             cGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAACAATTGGC‑GGGGTTGTTcgc                                                                                                                                                >  1:122364/1‑147 (MQ=255)              gTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGC‑GTTcc                                                                                                                                          <  2:806627/151‑1 (MQ=255)              gTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGC‑GTTcc                                                                                                                                          <  2:771437/151‑1 (MQ=255)              gTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGGTTTCttcc                                                                                                                                                >  2:643910/1‑143 (MQ=255)              gTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGGTTTCTccc‑ggt                                                                                                                                            >  2:71738/1‑147 (MQ=255)               ttCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGGTTGCTCGC‑GGTcc                                                                                                                                          >  1:435777/1‑150 (MQ=255)                   gATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGC‑GTTccgcc                                                                                                                                       <  1:701405/149‑1 (MQ=255)                        aGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGC‑GTTCCGCCGCCCGt                                                                                                                                 <  2:669092/150‑1 (MQ=255)                         gCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGGTTTCTCGC‑GTTCCGccccccccg                                                                                                                                >  2:758291/1‑147 (MQ=255)                            ggCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGC‑GTTCCGCCGCCCGTGC‑tt                                                                                                                            <  1:170194/150‑1 (MQ=255)                                   ggTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGC‑GTTCCGCCGCCCGTGC‑TTTAACCGTa                                                                                                                    <  2:812999/151‑1 (MQ=255)                                          aTAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGGTTTCTCCC‑GTTCCCCCGCCCCTGC‑CTTtacccta‑ccc                                                                                                                >  1:813061/1‑147 (MQ=255)                                           tAATGGTGTTCTGCTCTTTTTCTATTTGGTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGGTTTCGCCC‑GGTCCGCCCCCCCGGC‑TTTTAcccgaacccggc                                                                                                             >  2:634045/1‑140 (MQ=255)                                            aaTGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGC‑GTTCCGCCGCCCGTGC‑TTTAACCGTA‑CCCAGGCg                                                                                                           <  2:945092/150‑1 (MQ=255)                                            aaTGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGGTTGGTCGC‑GTTCCGCCGCCCCGGC‑TTTTACCGGA‑CCCcggcgg                                                                                                          >  2:126107/1‑151 (MQ=255)                                            aaTGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGGTTGCTTGC‑GGTTCGCCCCCCCTGC‑CTTTACCGTA‑CCCcggcgg                                                                                                          >  2:65011/1‑151 (MQ=255)                                             aTGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGGTTGCTCGC‑GTTCCCCCCCCCGTGC‑TTTAAACGGA‑CCCAGGCggg                                                                                                         >  1:990586/1‑151 (MQ=255)                                               ggTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGC‑GTTCCGCCGCCCGTGC‑TTTAACCGTA‑CCCAGGCGGGc                                                                                                        <  1:499331/150‑1 (MQ=255)                                               ggTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGC‑GTTCCGCCGCCCGTGC‑TTTAACCGTA‑CCCAGGCGGGc                                                                                                        <  1:123796/150‑1 (MQ=255)                                                  gtTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGGTTTCTCGCGGTTCGCCCCCCCGGCC‑TTTAACC‑CA‑CCCCGGCggg                                                                                                         >  2:577975/1‑146 (MQ=255)                                                     cTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGC‑GGTCCGCCCCCCCCGC‑TTTTACCCTA‑CCCCGGCGGGCAGT‑aa                                                                                                  >  2:50413/1‑149 (MQ=255)                                                       gCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTTCTCGC‑GGTCCCCCCCCCCCGG‑TTTTAACCTA‑CCCACGgggggcggaacc                                                                                                 >  2:312403/1‑141 (MQ=255)                                                       gCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGC‑GTTCCGCCGCCCGTGC‑TTTAACCGTA‑CCCAGGCGGGCAGT‑AACCT‑c                                                                                             <  2:239613/151‑1 (MQ=255)                                                           tttttCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGGTTGCTCGC‑GGTCCCCCGCCCCTGC‑CTTTACCGCA‑CCCAGGCGGGCAGT‑AACCT‑Ccccc                                                                                         >  1:546516/1‑151 (MQ=255)                                                             tttCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGGTTGCTCGC‑GTTCCGCCGCCCCCGC‑TTTTACCCCA‑CCCCGGCGGGCCGG‑AAACT‑ccccccc                                                                                       >  2:849963/1‑150 (MQ=255)                                                                 tATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGGTTGCTCCC‑GTTTCGCCGCCCCGGC‑TTTTAACCTA‑CCCACGCGGGGAGT‑AACCT‑cccccccaat                                                                                    >  2:583202/1‑146 (MQ=255)                                                                  aTTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGC‑GTTCCGCCGCCCGTGC‑TTTAACCGTA‑CCCAGGCGGGCAGT‑AACCT‑CCGCCCGAtt                                                                                    <  2:1059984/149‑1 (MQ=255)                                                                   tttGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGGTTTCTCGC‑GGTCCCCCGCCCCTGCTTTTAACCGAA‑CCC‑GGCGGGCCGG‑AAACC‑CCCCCCGcttctt                                                                                 >  2:684490/1‑151 (MQ=255)                                                                   tttGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGGTTGCTCCC‑GTTCCGCCGCCCGTGC‑CTTTACCCCA‑CCCCGGGGGGCAGT‑AACCT‑CCGccccgctctt                                                                                 >  2:404152/1‑144 (MQ=255)                                                                    ttGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTTCTCGC‑GGTCCGCCCCCCGTGC‑TTTAACCGTA‑CCCCGGCGGGCAGT‑AACCT‑CCGCCCGATTCttt                                                                                >  1:1003992/1‑151 (MQ=255)                                                                    ttGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGC‑GTTCCGCCGCCCGTGC‑TTTAACCGTA‑CCCAGGCGGGCAGT‑AACCT‑CCGCCCGATTCtt                                                                                 <  1:229463/150‑1 (MQ=255)                                                                     tgctgcAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGGTTGGTCCC‑GGTCCGCCGCCCCTGCGTTTAACC‑TA‑CCCACGCGGGCAGT‑AACCC‑CCCCCCCGTTCttt                                                                                >  2:850147/1‑150 (MQ=255)                                                                     tgctgcAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGATACGCCTCCAGCCAGCAATTGGC‑GGGGTTGCTCGC‑CTTCCGCCCCCCCTGC‑GTTTAACGTA‑CCCCAGAGGGCGAA‑AACCA‑CCGCCCCATTCtttt                                                                               >  2:679010/1‑150 (MQ=255)                                                                      gctgcAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGGTTCGC‑‑TTTCGCCCCCCCTGC‑TTTTACCCCA‑CCCCGGCGGGGAGG‑CAACT‑CCCCCCCGTTTTTTcc                                                                              >  2:206566/1‑151 (MQ=255)                                                                      gctgcAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGC‑GTTCCGCCGCCCGTGC‑TTTAACCGTA‑CCCAGGCGGGCAGT‑AACCT‑CCGCCCGATTCTTTcc                                                                              <  2:505710/151‑1 (MQ=255)                                                                        tgcAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGGTTTCTCGC‑GTTCCCCCGCCCCTGC‑TTTTACCGTA‑CCCCCGCGGGCAGT‑AACCC‑CCCCCCCATTCTTTcccc                                                                            >  2:342690/1‑151 (MQ=255)                                                                        tgcAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGGTTTCTCGC‑GGTTCGCCCCCCCTGC‑TTTTACCGTA‑CCCCGGCGGGGAGG‑AAACT‑CCCCCCCCGTCTTTtccc                                                                            >  2:1022721/1‑151 (MQ=255)                                                                          cAACACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGC‑GTTCCGCCGCCCGTGC‑TTTAACCGTA‑CCCAGGCGGGCAGT‑AACCT‑CCGCCCGATTCTTTCCCCAg                                                                          <  1:81238/151‑1 (MQ=255)                                                                           aaGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGC‑GTTCCGCCGCCCGTGC‑TTTAACCGTA‑CCCAGGCGGGCAGT‑AACCT‑CCGCCCGATTCTTTCCCCAGa                                                                         <  1:745548/151‑1 (MQ=255)                                                                             gACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTTTACCC‑GTTCCCCCGCCCCCGC‑TTTTAACGTACCCCCGGGGGGCGGA‑AACCC‑CCCCCCCATTTTTTtcc                                                                             >  2:748752/1‑146 (MQ=255)                                                                               ccGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑CGGGTTGCCCGC‑GGTCCGCCGCCCCTGC‑TTTAACCGTA‑CCCAGGTGGGGAGG‑AAACT‑CCCCCCCATTCTTTtcccc                                                                           >  1:12457/1‑144 (MQ=255)                                                                                cGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGGTTGCTCCC‑GGTCCGCCGCCCCGGC‑TTTAACCGCA‑CCCCGGGGGGCCGG‑AACCT‑CCGCCCCATTTTTTTCCCCGATATc                                                                     >  1:14948/1‑150 (MQ=255)                                                                                   gtTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGC‑GTTCCGCCGCCCGTGC‑TTTAACCGTA‑CCCAGGCGGGCAGT‑AACCT‑CCGCCCGATTCTTTCCCCAGATATCgag                                                                  <  1:467547/150‑1 (MQ=255)                                                                                     ttAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTTCTCGC‑GTTCCGCCGCCCCTGC‑TTTAAACCTA‑CCCAGGCGGGCCGT‑AAACT‑CCGCCCGATTCTTTTCCCCGATATCGAGATg                                                               >  1:683556/1‑151 (MQ=255)                                                                                        aCGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGGTTGCTCGC‑GGTCCCCCCCCCCCGG‑TTTTAACGAA‑CCCAAGCGGGGATT‑AAACTCCCGCCCGATTTTTTCCCCaaaattcggattgc                                                              >  1:729592/1‑137 (MQ=255)                                                                                          gccgccTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGC‑GTTCCGCCGCCCGTGC‑TTTAACCGTA‑CCCAGGCGGGCAGT‑AACCT‑CCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGc                                                          <  1:507447/151‑1 (MQ=255)                                                                                            cgccTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGGTTTCCCGC‑GTTTCGCCCCCCGCGC‑TTTTACCCTA‑CCCCGGCGGGGCGT‑AAACT‑CCGCCCCGTTTTTTTCCCCGATATCAAGGTGCCAGccc                                                        >  2:495706/1‑151 (MQ=255)                                                                                             gccTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGC‑GTTTCGCCCCCCGTGC‑TTTTACCCTA‑CCCACGCGGGCAGA‑AACCT‑CCCCCCGCTTTTTTCCCCCGatattgggaaggcacacc                                                        >  1:690145/1‑136 (MQ=255)                                                                                              ccTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGC‑GTTCCGCCGCCCGTGC‑TTTAACCGTA‑CCCAGGCGGGCAGT‑AACCT‑CCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCccgc                                                      <  1:35108/151‑1 (MQ=255)                                                                                              ccTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGGTTTCTCGC‑GTTTCGCCGCCCGGGC‑TTTAACCGTA‑CCCCGGCGGGGAGG‑AACCC‑CCCCCCGCTTCTTTTCCCCGATATCGAGATgcccgccccc                                                      >  2:852230/1‑149 (MQ=255)                                                                                               cTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGGTTGCTTCC‑GTTCCGCCGCCCGTGC‑TTTTAACGGC‑ACCAGGCGGGGAGT‑AAACC‑CCGCCCCAGTTTTTTCCCAGCTATTGAGATGCCAGCccgcc                                                     >  1:710738/1‑151 (MQ=255)                                                                                                tGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGGTTGCTCGC‑GTTCCCCCGCCCCGGCTTTTAACCG‑A‑CCCCGGCGGGGAATAAACCC‑CCCCCCGTTTTTTTCCCC‑GATATCGAGATGCCAGCccgcc                                                     >  2:973511/1‑150 (MQ=255)                                                                                                 gggTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCGGCCAGCAATTGGC‑GGGGTTTCTCGC‑GGTCCGCCGCCCCCGC‑TTTTACCGTA‑ACCCAGCGGGCCGG‑AAACT‑CCCCCCCCTTTTTTCCCCCGATATCGCGATGCCCGcccccc                                                     >  2:880834/1‑149 (MQ=255)                                                                                                 gggTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGC‑GTTCCGCCGCCCGTGC‑TTTAACCGTA‑CCCAGGCGGGCAGT‑AACCT‑CCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCccgcc                                                     <  2:278963/149‑1 (MQ=255)                                                                                                 gggTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGC‑GTTCCGCCGCCCGTGC‑TTTAACCGTA‑CCCAGGCGGGCAGT‑AACCT‑CCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGt                                                   <  1:709551/151‑1 (MQ=255)                                                                                                     ttGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGGTTGCTCTC‑GTTCCGCCGCCCCTGC‑TTTTACCGTA‑CCCAGGCGGGCAGT‑AACCC‑CCCCCCCATTCTTTTCCCAGATATCGAGATGCCAGCCCGCCGtt                                                  >  1:975453/1‑148 (MQ=255)                                                                                                             cccAGTACGGGCGGGCTCCGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGC‑GTTCCGCCGCCCGTGC‑TTTACCCGTA‑CCCAGGCGGGCAGT‑ACCCT‑CCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAg                                        <  1:821622/150‑1 (MQ=255)                                                                                                                     aGGGGGGCTACGCTGCCAGCCAGCAATGGGC‑GGGTTTGCTCGC‑GTTCCGCCGCCCGTGC‑TTTAACCGTA‑CCCAGGCGGGCATT‑ACCCT‑CCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTg                                <  2:805521/150‑1 (MQ=255)                                                                                                                          ggCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGC‑GTTCCGCCGCCCGTGC‑TTTAACCGTA‑CCCAGGCGGGCAGT‑AACCT‑CCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTa                           >  2:833397/1‑150 (MQ=255)                                                                                                                               cGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGC‑GTTCCGCCGCCCGTGC‑TTTAACCGTA‑CCCAGGCGGGCCGT‑AACCT‑CCGCCCGATTCTTTCCCCAGATATCGAGATGCCCGCCCGCCGTTAAGCCAAAAGTACCCTTTGTGTAcc                         >  2:161288/1‑147 (MQ=255)                                                                                                                                    ccagccagCAATGGGC‑GGTTTTGCTCGC‑GTTCGCCCGCCGGTCC‑TTTACCCGGC‑CGCAGGGGGGAAGA‑ACCCC‑CCGCCCGATGCTTCCCCCAGATATCGAGATGCCGGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCg                 <  2:552906/150‑1 (MQ=255)                                                                                                                                           gcAATTGGC‑GGGTTTGCTCGC‑GTTCCGCCGCCCGTGC‑TTTAACCGTA‑CCCAGGCGGGCAGT‑AACCT‑CCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTa           <  1:758412/149‑1 (MQ=255)                                                                                                                                             aaTTGGC‑GGGTTTGCTCGC‑GTTCCGCCGCCCGTGC‑TTTAACCGTA‑CCCAGGCGGGCAGT‑AACCT‑CCGCCCGATTCTTTTCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCCAAAGTACCGTTGGTGTACCCCGGGCCGGGCGTACCt        >  2:484510/1‑150 (MQ=255)                                                                                                                                                 ggC‑GGTTTGGCTCGC‑GTCCCGCCGCCGGTGC‑TTTACCCGTC‑CCCAGGCGGGCAGT‑ACCCC‑CCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGc      <  2:656111/148‑1 (MQ=255)                                                                                                                                                     gggTTTGCTCGC‑GTTCCGCCGCCCGTGC‑TTTAACCGTC‑CCCAGGGGGGAAGT‑ACCCT‑CCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGcc  <  1:770358/149‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                       AAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGC‑GTTCCGCCGCCCGTGC‑TTTAACCGTA‑CCCAGGCGGGCAGT‑AACCT‑CCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCC  >  NC_000913/2228072‑2228367

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