Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,508,657	0	A	G	42.5%	120.1 / 23.6	81	F386S (TTC→TCC)	yhfW	phosphopentomutase‑related metalloenzyme superfamily protein
Reads supporting (aligned to +/- strand):  ref base A (6/40);  new base G (34/0);  total (41/40)
Fisher's exact test for biased strand distribution p-value = 1.74e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGT‑AGCCCCCAC‑ATCAGAAAGCGTGGTGC‑GCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCC  >  NC_000913/3508512‑3508806                                                                                                                                                   |                                                                                                                                                          ccACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGaaaaa                                                                                                                                                       <  2:1020917/150‑1 (MQ=255)    cGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGG‑GGGAAAAATc                                                                                                                                                    >  2:198908/1‑150 (MQ=255)     gttgttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAACCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGCGGGGCGG‑GGgaaaa                                                                                                                                                       >  2:1183957/1‑146 (MQ=255)       tgttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGgaaaa                                                                                                                                                       >  1:559527/1‑144 (MQ=255)       tgttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGAGGCGCGA‑GGgaaaa                                                                                                                                                       >  1:18664/1‑144 (MQ=255)           tttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATCac‑acac                                                                                                                                             >  1:661122/1‑149 (MQ=255)              aaaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGCAGCCC‑GGGAGAATTCac‑acacacc                                                                                                                                          >  2:357445/1‑146 (MQ=255)                  ccGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGCGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGGGA‑GGGAAAACTCAA‑ACACCAT‑GGcccc                                                                                                                                   >  1:658551/1‑151 (MQ=255)                  ccGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATTAC‑AAACCGT‑TGcccc                                                                                                                                   >  1:1246695/1‑151 (MQ=255)                   cGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGT‑AGccccc                                                                                                                                  <  2:991259/151‑1 (MQ=255)                   cGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGAGACGCGC‑GGGAGAAAAcac‑acaagcaagacccc                                                                                                                                   >  2:885101/1‑140 (MQ=255)                     gagagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGT‑AGCCCCcac                                                                                                                                <  1:576652/151‑1 (MQ=255)                       gagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGT‑AGCCCCCAC‑At                                                                                                                             <  1:1040773/151‑1 (MQ=255)                       gagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGT‑AGCCCCCAC‑At                                                                                                                             <  2:347335/151‑1 (MQ=255)                        agaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGGGC‑CGGAAAAAAAAC‑ACACCGC‑AGcccccc                                                                                                                                 >  1:275053/1‑146 (MQ=255)                          aTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATTAA‑ACACGGT‑AGCCCCCAC‑ATCa                                                                                                                           >  1:202323/1‑150 (MQ=255)                          aTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTAGGGGGC‑GGGAAAAATTAA‑AAACGGG‑AGACCCCCC‑AACa                                                                                                                           >  1:698584/1‑150 (MQ=255)                           tCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATTAC‑ACAACGT‑AACCCCCCC‑ATc                                                                                                                            >  1:139307/1‑148 (MQ=255)                            cAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGG‑GGGAAAAATAAC‑ACACGGC‑AGGCCCccc‑acc                                                                                                                            >  2:825333/1‑145 (MQ=255)                             aaTATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATCAC‑ACACCGT‑AACCCCCCC‑CTCAGa                                                                                                                         >  1:263306/1‑149 (MQ=255)                                atGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGT‑AGCCCCCAC‑ATCAGAAAGc                                                                                                                     <  2:10811/150‑1 (MQ=255)                                  gCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGT‑AGCCCCCAC‑ATCAGAAAGCgt                                                                                                                   <  2:896504/150‑1 (MQ=255)                                   cAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATTAA‑ACACCGT‑AGGCCCCCC‑CCCCGAAAGCggg                                                                                                                  >  1:928025/1‑148 (MQ=255)                                   cAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGT‑AGCCCCCAC‑ATCAGAAAGCgtg                                                                                                                  <  2:405931/150‑1 (MQ=255)                                        gCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGT‑AGCCCCCAC‑ATCAGAAAGCGTGGTgc‑g                                                                                                            <  2:844414/150‑1 (MQ=255)                                        gCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGT‑AGCCCCCAC‑ATCAGAAAGCGTGGTgc‑g                                                                                                            <  1:191487/150‑1 (MQ=255)                                          tCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATTAC‑AAACAGT‑AACCCCCCC‑ATCCGaaaaacggggggg                                                                                                             >  2:298099/1‑138 (MQ=255)                                            aTGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATCAC‑ACAACGT‑AGGCCCCCC‑ACCAAAAAGAGGGGGgc‑g                                                                                                            >  1:1098533/1‑146 (MQ=255)                                                 cGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGG‑GGTGAAAATCAC‑ACACGGT‑AGCCCCCAC‑ATCAGAACGCGGGGGGA‑Gcaga                                                                                                        >  2:146555/1‑143 (MQ=255)                                                  gAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGCGGCGCGC‑GGGGAAAAATAC‑ACACGGA‑AGACCCCCC‑AACCGAAAGAGTGGGGT‑GGACAccaa                                                                                                    >  1:466519/1‑146 (MQ=255)                                                      gtgtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGT‑AGCCCCCAC‑ATCAGAAAGCGTGGTGC‑GCACACCGAGCTGCGt                                                                                             <  1:238500/151‑1 (MQ=255)                                                       tgtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGGGC‑GGGAAAAATTAC‑ACACGGG‑AAACCCCCC‑ATCCCAAAAAGGGGGGG‑GGACACCGAGCcgcg                                                                                              >  2:1033328/1‑149 (MQ=255)                                                       tgtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATCAC‑ACACGGT‑AGCCCCCAC‑ATCAGAAAACGGGGGGC‑GCACAcc                                                                                                      >  1:695563/1‑141 (MQ=255)                                                       tgtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATTAC‑ACACGGG‑AGCCCCCAC‑ACCAGAAAGAGGGGGGG‑GCCCACCCAGCCGc                                                                                               >  2:451846/1‑148 (MQ=255)                                                        gtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGT‑AGCCCCCAC‑ATCAGAAAGCGTGGTGC‑GCACACCGAGCTGCGTa                                                                                            <  1:1377372/150‑1 (MQ=255)                                                        gtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGT‑AGCCCCCAC‑ATCAGAAAGCGTGGTGC‑GCACACCGAGCTGCGTAg                                                                                           <  2:775336/151‑1 (MQ=255)                                                          aaCCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATTAC‑ACACGGG‑AGCCCCCCC‑ATCAGAAAGCGGGGGGG‑GCACACCCAGCTGCGGAg                                                                                           >  1:1127693/1‑149 (MQ=255)                                                             ccGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATTAC‑ACACGGT‑AGCCCCCAC‑CTCAGAAAGCGTGGGGG‑GCACACCCAGCCGCGTAGCg                                                                                         >  1:611427/1‑148 (MQ=255)                                                              cGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGT‑AGCCCCCAC‑ATCAGAAAGCGTGGTGC‑GCACACCGAGCTGCGTAGCGatca                                                                                     <  2:801816/151‑1 (MQ=255)                                                               gTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATTAC‑ACACGGG‑AGGCCCCCC‑ATCACAAAACGGGGGGC‑GCACACCCAGCTGCGTAGCGatca                                                                                     >  1:1422103/1‑150 (MQ=255)                                                               gTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATCAC‑ACACCGG‑AGCCCCCCC‑CCCAAAAAACGGGGGGG‑GCACACCCAGCTggggaggggtca                                                                                     >  2:249818/1‑139 (MQ=255)                                                                     tCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATTAC‑ACACGGG‑AACCCCCCC‑ATCAGAAAACGGGGGGC‑GCCCACCCAGCTGCGTAGCGATCATCCCtt                                                                               >  1:1210030/1‑150 (MQ=255)                                                                      cAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATCAC‑ACACGGT‑AGGCCCCAC‑ATCAGAAAACGTGGGGC‑GCACAACCGGCGGCGTAGCGATCATCCCCTg                                                                              >  1:1228103/1‑150 (MQ=255)                                                                      cAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGCGGCGCGG‑GGAGAAAATAACAACAC‑GA‑AGCCCCCCC‑ATAACAAAAAGGGGGGG‑GCGCACCAAGCCTCGGTGGGCTCATCCCTTg                                                                              >  1:470693/1‑150 (MQ=255)                                                                       aaaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGT‑AGCCCCCAC‑ATCAGAAAGCGTGGTGC‑GCACACCGAGCTGCGTAGCGATCATCCCTTgat                                                                            <  1:582688/151‑3 (MQ=255)                                                                        aaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGT‑AGCCCCCAC‑ATCAGAAAGCGTGGTGC‑GCACACCGAGCTGCGTAGCGATCATCCCTTGCTg                                                                           <  2:333030/151‑1 (MQ=255)                                                                             cATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATCAC‑ACACAGT‑AGCCCCCACAATCA‑AAAAAGGGGGGG‑GCACACCGAGCCGCGTAGCGATCATCCCCTTCTGATaa                                                                       >  2:1202893/1‑150 (MQ=255)                                                                             cATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGT‑AGCCCCCAC‑ATCAGAAAGCGTGGTGC‑GCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATaa                                                                       <  2:588034/150‑1 (MQ=255)                                                                               tAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAATTTAA‑ACACGGG‑AGCCCCCAC‑AACAGAAAGCGTGGTGG‑GGACACCCAGCTGCGTAGCGATCATCCCTTGCTGATaaa                                                                      >  1:846324/1‑149 (MQ=255)                                                                                aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGT‑AGCCCCCAC‑ATCAGAAAGCGTGGTGC‑GCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAAcc                                                                    <  1:9774/150‑1 (MQ=255)                                                                                aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGT‑AGCCCCCAC‑ATCAGAAAGCGTGGTGC‑GCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAAcc                                                                    <  1:110539/150‑1 (MQ=255)                                                                                aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGT‑AGCCCCCAC‑ATCAGAAAGCGTGGTGC‑GCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCa                                                                   <  1:502607/151‑1 (MQ=255)                                                                                 gggTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGT‑AGCCCCCAC‑ATCAGAAAGCGTGGTGC‑GCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAg                                                                  <  2:793536/151‑1 (MQ=255)                                                                                   gtgtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGT‑AGCCCCCAC‑ATCAGAAAGCGTGGTGC‑GCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCa                                                                <  2:1095063/151‑1 (MQ=255)                                                                                     gtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGT‑AGCCCCCAC‑ATCAGAAAGCGTGGTGC‑GCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCAc                                                               <  1:485346/150‑1 (MQ=255)                                                                                     gtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGT‑AGCCCCCAC‑ATCAGAAAGCGTGGTGC‑GCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACt                                                              <  2:33148/151‑1 (MQ=255)                                                                                       ctcCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATCTA‑ACACGGG‑AGCCCCCAC‑AACAGCAAACGTGGGGG‑GCCCACCCAGCTGCGGAGCGGTCATCCCCTTCTGATTAAACAACACCgg                                                            >  2:914969/1‑151 (MQ=255)                                                                                         ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTTGGGCGG‑GGGGAAAATTAC‑AAACGGTAAGCGCCCCC‑TTAAAAAAG‑GAGGGGG‑GGACACCCAGCCGCGGAGGAACCATCCCCTGTTGATAAAACacaacagg                                                            >  2:263041/1‑142 (MQ=255)                                                                                             cAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGT‑AGCCCCCAC‑ATCAGAAAGCGTGGTGC‑GCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGcaccac                                                      <  1:981242/151‑1 (MQ=255)                                                                                               aaCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGTGCGCGC‑GGAAAAAATAAC‑AAAAAGG‑AGCCCCCCC‑AACAAGAAGCGTGGGGCGGCTAACCG‑GGAGCGGCGGGATCCTCCCCTGCTGATAAAACAACACCGGcacc                                                        >  2:201239/1‑147 (MQ=255)                                                                                                aCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAATTTAC‑ACACCGT‑AGCCCCCAC‑ATCAGAAAACGCGGTGC‑GCACCCCCACCGGCGTAGCGATCATCCCCTGCTGGTTAACCAGCACCGGGACCCCtt                                                    >  1:241450/1‑150 (MQ=255)                                                                                                 ccGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGT‑AGCCCCCAC‑ATCAGAAAGCGTGGTGC‑GCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTcg                                                  <  1:1406485/151‑1 (MQ=255)                                                                                                  cGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGT‑AGCCCCCAC‑ATCAGAAAGCGTGGTGC‑GCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTcgc                                                 <  1:1113783/151‑1 (MQ=255)                                                                                                   ggagggagGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATCAC‑ACACCGG‑AGCCCCCCC‑ATCAGAAAACGGGGGGG‑GCACACCCCGCTGCGTAGCGATCATCCCCTGCTGATAAAACAACACCGGCACCACTTcgc                                                 >  1:717025/1‑150 (MQ=255)                                                                                                   ggagggagGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGCGC‑GGGAAAAATAAA‑ACACAGG‑ATGCCCCCC‑CTCAGAAAACGTGGGGC‑GCCCCCCGAGCCGCGTAGGGGTCATTCCTTGCTGATAAACCAGCACAGGggccaccttgc                                                 >  1:385231/1‑145 (MQ=255)                                                                                                    gagggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGT‑AGCCCCCAC‑ATCAGAAAGCGTGGTGC‑GCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCgg                                               <  2:602609/151‑1 (MQ=255)                                                                                                            aaagaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTAAC‑‑CACCGG‑AGCCCCCCC‑ACCAGAAAGCGTGGGGG‑GCACCCCCAGCCGCGTGGCGATCATCCCCTGCTGGTAAACCAGCACCCGGGCCACCTCGCGGGGAtgg                                         >  2:1179098/1‑149 (MQ=255)                                                                                                            aaagaaaagaGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGT‑AGCCCCCAC‑ATCAGAAAGCGTGGTGC‑GCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTAtggtg                                       <  2:1411313/151‑1 (MQ=255)                                                                                                            aaagaaaagaGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGT‑AGCCCCCAC‑ATCAGAAAGCGTGGTGC‑GCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTAtggtg                                       <  1:1189243/151‑1 (MQ=255)                                                                                                               gaaaagaGCGACCATTTTGTGGCGGTGGCGGGC‑GGGAAAAATTAC‑ACACCGG‑GACCCCCCC‑CTCATAAAGCGGGGGGG‑GGACACCGAGGCGCGGAGCCATCATCCCTCTCTGCAAAACCAGCACCGGGACCCCCTCGCggggcaggtgg                                      >  1:904078/1‑141 (MQ=255)                                                                                                               gaaaagaGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGT‑AGCCCCCAC‑ATCAGAAAGCGTGGTGC‑GCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCt                                    <  2:18664/151‑1 (MQ=255)                                                                                                                 aaagaGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATTAC‑ACACCGT‑AGCCCCCCC‑CTCAGAAAACGGGGGGC‑GCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAAACAACACTGGCACCACTTCGCGGGGATggtgggtg                                   >  1:516561/1‑150 (MQ=255)                                                                                                                    gagCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGT‑AGCCCCCAC‑ATCAGAAAGCGTGGTGC‑GCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGAc                               <  1:1073230/151‑1 (MQ=255)                                                                                                                      gCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGT‑AGCCCCCAC‑ATCAGAAAGCGTGGTGC‑GCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGAcc                              <  2:1119584/150‑1 (MQ=255)                                                                                                                        gACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGT‑AGCCCCCAC‑ATCAGAAAGCGTGGTGC‑GCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAAt                           <  1:1319150/151‑1 (MQ=255)                                                                                                                           cATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGT‑AGCCCCCAC‑ATCAGAAAGCGTGGTGC‑GCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGt                        <  2:187927/151‑1 (MQ=255)                                                                                                                               ttGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGT‑AGCCCCCAC‑ATCAGAAAGCGTGGTGC‑GCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCgg                     <  2:985482/150‑1 (MQ=255)                                                                                                                                  tggcggtggcgCGC‑GGAAAAATTCAC‑ACACGGT‑AGCCCCCAC‑ATCAGAAAGCGTGGTGC‑GCACACCGAGCTGCCTAGCGATCATCCCTTGCTGATAAACCAGCCCTGGCCCCACTTCGCGGGTATGGTGGCCGTGACCAATGGTCGGa                    >  2:947623/1‑148 (MQ=255)                                                                                                                                  tggcggtagcgGGC‑GCAAAACGTCAC‑ACACAGT‑CGTCAGCAC‑ATCAGAAAGCGTGGTGC‑GCACACCGAGCTGCATTGCGATCGTCCCTTCCTGGTAACCAAGCACTGGCACCATTTCGCGGGTATGGTGGCTGTGCCCTGTGTTCgg                     <  2:103002/147‑1 (MQ=255)                                                                                                                                    gcggtggcgCGC‑GGAAAAATTCAC‑ACACGGT‑AGCCCCCAC‑ATCAGAAAGCGTGGTGC‑GCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGtt               <  1:363637/151‑1 (MQ=255)                                                                                                                                     cggtggcgCGC‑GGAAAAATTCAC‑ACACGGT‑AGCCCCCAC‑ATCAGAAAGCGTGGTGC‑GCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTg              <  1:160072/151‑1 (MQ=255)                                                                                                                                      ggtggcgCGC‑GGAAAAATTCAC‑ACACGGT‑AGCCCCCAC‑ATCAGAAAGCGTGGTGC‑GCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGc             >  1:1102094/1‑151 (MQ=255)                                                                                                                                          gcgcgc‑GGAAAAATTCAC‑ACACGGT‑AGCCCCCAC‑ATCAGAAAGCGTGGTGC‑GCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTg         <  2:1425406/151‑1 (MQ=255)                                                                                                                                          gcgcgc‑GGAAAAATTCAC‑ACACGGT‑AGCCCCCAC‑ATCAGAAAGCGTGGTGC‑GCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTg         <  1:930433/151‑1 (MQ=255)                                                                                                                                                  gAAAAATTCAC‑ACACGGT‑AGCCCCCAC‑ATCAGAAAGCGTGGTGC‑GCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGcc  >  1:1116127/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                          CCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGT‑AGCCCCCAC‑ATCAGAAAGCGTGGTGC‑GCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCC  >  NC_000913/3508512‑3508806

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