Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,508,662	0	T	A	40.7%	131.8 / 38.6	83	E384D (GAA→GAT)	yhfW	phosphopentomutase‑related metalloenzyme superfamily protein
Reads supporting (aligned to +/- strand):  ref base T (8/40);  new base A (33/0);  total (42/41)
Fisher's exact test for biased strand distribution p-value = 2.12e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGT‑AGCCCCCAC‑ATCAGAAAGCGTGGTGC‑GCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCAC  >  NC_000913/3508515‑3508811                                                                                                                                                     |                                                                                                                                                          cGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGG‑GGGAAAAATc                                                                                                                                                         >  2:198908/1‑150 (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                                                                                                                                  <  2:347335/151‑1 (MQ=255)                    gagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGT‑AGCCCCCAC‑At                                                                                                                                  <  1:1040773/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:110539/150‑1 (MQ=255)                                                                             aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGT‑AGCCCCCAC‑ATCAGAAAGCGTGGTGC‑GCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAAcc                                                                         <  1:9774/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                                            <  1:1189243/151‑1 (MQ=255)                                                                                                         aaagaaaagaGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGT‑AGCCCCCAC‑ATCAGAAAGCGTGGTGC‑GCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTAtggtg                                            <  2:1411313/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              <  1:930433/151‑1 (MQ=255)                                                                                                                                       gcgcgc‑GGAAAAATTCAC‑ACACGGT‑AGCCCCCAC‑ATCAGAAAGCGTGGTGC‑GCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTg              <  2:1425406/151‑1 (MQ=255)                                                                                                                                               gAAAAATTCAC‑ACACGGT‑AGCCCCCAC‑ATCAGAAAGCGTGGTGC‑GCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGcc       >  1:1116127/1‑151 (MQ=255)                                                                                                                                                    aTTCAC‑ACACGGT‑AGCCCCCAC‑ATCAGAAAGCGTGGTGC‑GCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCAc  <  1:167707/151‑1 (MQ=255)                                                                                                                                                    aTTCAC‑ACACGGT‑AGCCCCCAC‑ATCAGAAAGCGTGGTGC‑GCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCAc  >  2:883423/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                          CGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAC‑ACACGGT‑AGCCCCCAC‑ATCAGAAAGCGTGGTGC‑GCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCAC  >  NC_000913/3508515‑3508811

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: ‑