Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,376,505	0	A	C	46.3%	115.6 / 34.2	82	T158P (ACC→CCC)	ycjR	putative TIM alpha/beta barrel enzyme
Reads supporting (aligned to +/- strand):  ref base A (35/9);  new base C (0/38);  total (35/47)
Fisher's exact test for biased strand distribution p-value = 6.67e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGCGG  >  NC_000913/1376358‑1376655                                                                                                                                                         |                                                                                                                                                       ttACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAACAccc                                                                                                                                                     >  1:1033360/1‑151 (MQ=255)   aCCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAACACCCt                                                                                                                                                    <  1:839094/150‑1 (MQ=255)     cgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAACACCCTCGc                                                                                                                                                 <  1:39469/151‑1 (MQ=255)     cgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAACACCCTCGc                                                                                                                                                 <  2:135773/151‑1 (MQ=255)          aTGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAACACCCTCGCCGATg                                                                                                                                            <  1:490167/151‑1 (MQ=255)          aTGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAACACCCTCGCCGATg                                                                                                                                            >  2:460504/1‑151 (MQ=255)                      cgTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTAc                                                                                                                                 >  2:772273/1‑150 (MQ=255)                         aGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACAtcg                                                                                                                             >  2:407080/1‑151 (MQ=255)                            ccggACGGCGACCGCAAAATGGGGAGGTTTTCCCTGCTCGTACGGAAACAGGTGGCCGCGGTGACCGGAC‑CC‑GTGT‑TGTATCTCGACCCTTAACC‑CCGCTCCCAGATTT‑ATATGATCAACCCCCTCGCCGATGCCCGCCGTTACAtcgtc                                                                                                                           <  1:269524/148‑1 (MQ=255)                              ggcgggCGCCCGCAAATGGGGGGGTTATTCCCTGCGGGGACGGGAACGGGTCGCCGCGCGTCCGGGAC‑CG‑GGGG‑TTTATCCCGACCCGTTAAA‑CCGTTATCAGGATA‑TTATGACAACCCCCCCCGCGGATGCCCGCCGTTACATCGTCGaa                                                                                                                        <  1:589514/147‑1 (MQ=255)                                 cGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAAc                                                                                                                     >  2:791792/1‑151 (MQ=255)                                  ggCGCCCGAAAAAGGGTGGGTGTTTCCGGGCGCTTACTGAAACGGTTCCCCGCGCGTCCCGGAC‑CC‑GTGG‑TGTCTCTCGACCTTTTAAC‑CCGCTATAAGAAAA‑ATAGGAAAAACCCCCCCGCGGAGGCCCGCCGTTACATCGTCGAAAACg                                                                                                                    <  1:980220/151‑1 (MQ=255)                                  ggCGACCGAAAAATGGTGAGTGATCCCCTGCGCGTACTGAAACGGGTCGCCGCGCGTCCCGAAC‑CG‑GTGG‑TGTATCCCGACCCGTAAAC‑CCGTTATCAGGATC‑ATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAAc                                                                                                                     <  2:610480/150‑1 (MQ=255)                                      aCCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCt                                                                                                                >  2:1020566/1‑151 (MQ=255)                                      aCCGAAAAGTGGTGAGTGATTCCCTGCGCGTACTGGAACAGGGCCCCGCGCGTCCCGGAC‑CG‑GTGG‑TGTATCTCGACCCGTTAAA‑CCGCTATCAGAATC‑ATATGATCACCCCCCGCGCCGATCCCCGCCGTTACATCGTCGAAAACGa                                                                                                                   <  2:173416/148‑1 (MQ=255)                                       ccGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGTTCtt                                                                                                               >  2:719416/1‑151 (MQ=255)                                       ccGAAAAGTGGTATGTGTTCCCCTGCGCGTACTGAAAGAGTCCGCCGCGCG‑ACCCGAA‑CCGGTGT‑TGTTCCTCTAA‑CGTTAAACCCGTTATCGGGATA‑ATTTGAACAACCCCCGCGCCGATGCCCGCCGTTACATCGTCGAAAACGa                                                                                                                   <  2:853284/147‑1 (MQ=255)                                            aatgggTGAGCGTTCCCCTGCGGGTACGGGAACAGGTCGCCGCGCGTCCCGAAA‑CG‑GGGG‑TTTACCTCGACCCTTTAAA‑CCGTTATAAGAATC‑ATATGATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACa                                                                                                          <  2:1059154/147‑1 (MQ=255)                                            aaaGGGTGAGTGTTTCCCGGCGCGTACTGGAACAGGTCGCCGCGCGTCCCGGAA‑CG‑GTGG‑TGTATCTCAAACCGTAAAC‑CCGCTATAAGGATC‑ATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACa                                                                                                          <  2:51622/151‑1 (MQ=255)                                             aaTGGTGGGTGATTCCCTGCGCGTACGGAAACGGGTCGCGGCGCGTCCCGGAA‑CC‑GGGG‑TTTACCCCGACCCTTAAAC‑CCGCTATCGGAATT‑ATATGATAAACCCCCTCGCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACAt                                                                                                         <  2:561754/151‑1 (MQ=255)                                                ggtgggtgTTCCCCGGCGCGTACGGGAACGGGTCGCCGCGGTTCCCGGCC‑CC‑GTGG‑TTTTTCCCGACCCGTTAAC‑CCGTTATCAGGATA‑ATATAAACAACCCCCTCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACAt                                                                                                         <  1:381785/148‑1 (MQ=255)                                                ggTGAGTGTTCCCCTGCGGTTACGGAAACGGGTCGCCGCGCGT‑CCGGAACCC‑GGGG‑TGTATCTCGACCCTTAAAC‑CCGCTATAAGGATC‑ATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGt                                                                                                       <  2:126778/150‑1 (MQ=255)                                                  tgagtgaTTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTAc                                                                                                     >  1:110023/1‑150 (MQ=255)                                                  tgagtgaTTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTAc                                                                                                     >  1:15678/1‑150 (MQ=255)                                                       gtttCCCTGCGGGTGCTGAAAGGGGTCGCCGCGCTTACCGGAA‑CC‑GTGG‑TGTATCTCGAACCGTTAAC‑CCGTTATCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTa                                                                                               <  2:94212/149‑1 (MQ=255)                                                       aatcccgtgggcttACTGAAACGGTTCGCCGCGGGTCCCGGAA‑CG‑GTGG‑TTTCTCTCAAGCCTTTAAC‑CCGCTATGGGGTTC‑ATTTAAACACCCCCCCCCCGGTGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTa                                                                                               <  2:792073/138‑1 (MQ=255)                                                         ttCCCGGCGCGTACTGGAACGGGTCCCCGCGGGTCCCGAAC‑CC‑GTGT‑TGTATCCCGACCCGTTAAC‑CCGCTATCGGGATA‑ATATAACCACCCCCCCCGCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTAt                                                                                              <  1:166003/150‑1 (MQ=255)                                                           cccTGCGCGTGCGGAAGCGGGCCGCCGCGCGCACCGGAC‑CG‑GTGG‑TGTACCTCGACCTGTTAAA‑CCGCTACCAGGATA‑ATATGAACAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATc                                                                                             <  1:103815/149‑1 (MQ=255)                                                           cccTGCGCGTACTGGAACAGGTCGCCGCGCGTCCCGGAC‑CG‑GTGG‑TGTATCTCGACCCTTAAAC‑CCTCTATAAGGATC‑ATATGATCAACCCCCTCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCgg                                                                                           <  2:431168/151‑1 (MQ=255)                                                           cccTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCgg                                                                                           >  2:417413/1‑151 (MQ=255)                                                           cccGGCGCTTACGGGACCAGGCCCCCGCGC‑TTCCCGAA‑CC‑GGGGTTTTCTCCCGACCCTTAAAC‑CCGCTACCAGAACA‑ATAAAAACACCCCCCCCGCCAATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCg                                                                                            <  1:282972/150‑1 (MQ=255)                                                            ccTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGc                                                                                          >  1:617797/1‑151 (MQ=255)                                                             cTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCg                                                                                         >  2:744161/1‑151 (MQ=255)                                                              tGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCg                                                                                         >  1:841361/1‑150 (MQ=255)                                                              tGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCg                                                                                         >  1:982038/1‑150 (MQ=255)                                                                    tACGGGAACAGGTCGCCGCGCGTACCGGAC‑CC‑GTGG‑TGTATCTCGAACCGTTAAC‑CCGCTATCAGGATC‑ATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCt                                                                                   <  2:68002/150‑1 (MQ=255)                                                                        ggAACGGGTCGCCGCGCGTCCCGGAC‑CG‑GTGG‑TGTATCTCGACCCGTTAAC‑CCGCTATCAGGATC‑ATATGACCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATcac                                                                              <  2:1040336/151‑1 (MQ=255)                                                                        gaaacgggcccccgcGCGTCCCGGAC‑CG‑GTGT‑TTTATC‑CCAACCCTTAAACCCGTTAACAGGATA‑ATATGACCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATcac                                                                              <  1:723673/140‑1 (MQ=255)                                                                               gTCGCCGCGCGTCCGGGAA‑CC‑GTGG‑TTTATCTCGAACCGTTAAC‑CCGCTATCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAata                                                                       <  2:231923/151‑1 (MQ=255)                                                                                 cgccgcGCGTCCCGAAC‑CG‑GTGT‑TTTATCCCAACCCGTTAAC‑CCGCTATCAGGATA‑ATATAACAACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAat                                                                        <  2:774549/148‑1 (MQ=255)                                                                                 cgccgcGCGTACCGGAA‑CC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATc                                                                     >  1:384026/1‑151 (MQ=255)                                                                                 cgccgcGCGTACCGGAA‑CC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATc                                                                     >  1:546085/1‑151 (MQ=255)                                                                                     gcgcgTACCGGAA‑CC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCga                                                                   >  2:163023/1‑149 (MQ=255)                                                                                      cgcgTACCGGAA‑CC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCga                                                                   >  2:918559/1‑148 (MQ=255)                                                                                         gtcccGGAC‑CC‑GTGT‑TGTTTCTCAACCCTTTAAC‑CCGCTATCAGGATC‑ATTTGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaagaa                                                               <  2:847457/146‑1 (MQ=255)                                                                                           cccGGAA‑CT‑GTGG‑TTTCTCGCGACCCGTAAAC‑CCGCTATCA‑GATCAATATAAACAACACCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGAt                                                            <  1:796388/149‑1 (MQ=255)                                                                                             cGGAACCC‑GTGG‑TGTTTCTCAAACCGTAAAC‑CCCCTTTAAGGTTC‑ATATAATCACCCCCCTCGCCAATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                                <  2:936506/145‑1 (MQ=255)                                                                                             cGGAA‑CC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGCTTTCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAAc                                                         >  1:447924/1‑151 (MQ=255)                                                                                               gaaccc‑GGGG‑TTTATCTC‑AACCCTTAAACCCGTTATCAGAATA‑ATATGATAACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGAt                                                            <  1:203563/143‑1 (MQ=255)                                                                                                aa‑CG‑GTGG‑TGTATCTCAACCCGTTAAC‑CCGCTATCAGGAAT‑ATATGATAACCACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAAc                                                         <  2:241249/148‑1 (MQ=255)                                                                                                aa‑CC‑GTGT‑TGTATCTCAAACCGTAAAC‑CCGCTATAAGAATC‑ATATAATAAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTg                                                      <  2:222445/151‑1 (MQ=255)                                                                                                aa‑CC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTg                                                      >  2:1071159/1‑151 (MQ=255)                                                                                                 a‑CC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTgg                                                     >  2:388668/1‑151 (MQ=255)                                                                                                    c‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcg                                                   >  1:749709/1‑151 (MQ=255)                                                                                                     cggtgttttATCGCGACCTGTAAAC‑CCGATAACAGAATA‑ATATGATAAACCCCCCCGCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATa                                                           <  1:127563/136‑1 (MQ=255)                                                                                                        gg‑ttTACCTCGACCCGTTAAC‑CCCCTATCAGAATC‑ATATGACCAACCCCCTCCCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCa                                                 <  2:305220/150‑1 (MQ=255)                                                                                                                tcccgaCCCTTAAACCCGTTATAAGGATC‑ATATGATCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgc                                              <  2:685662/142‑1 (MQ=255)                                                                                                                ctcGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGc                                         >  2:741051/1‑151 (MQ=255)                                                                                                                ctcGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGc                                         >  2:88878/1‑151 (MQ=255)                                                                                                                 tcGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGc                                         >  2:216673/1‑150 (MQ=255)                                                                                                                 ctcaaCCGTTAAACCCCCTATCAGTATA‑ATATGAAAACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgcgc                                            <  2:542760/145‑1 (MQ=255)                                                                                                                 ccggaCCCTTAAACCCACTATCAGCATC‑ATATAACCACCCCCCGCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTg                                          <  2:541106/146‑1 (MQ=255)                                                                                                                  cGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCAt                                       >  2:479276/1‑151 (MQ=255)                                                                                                                     aCCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGAc                                    >  1:242832/1‑151 (MQ=255)                                                                                                                         ttAAA‑CCGCTATCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAAcc                                >  1:281846/1‑151 (MQ=255)                                                                                                                           aaa‑CCCCTATCAGGTTC‑ATATTACCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcg                               <  2:423626/150‑1 (MQ=255)                                                                                                                             acccacTACCAGGATA‑TTATAACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgcg                             <  1:334775/146‑1 (MQ=255)                                                                                                                               ccGCTATCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGAc                           >  1:40410/1‑151 (MQ=255)                                                                                                                               ccGCTATCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGAc                           >  2:117406/1‑151 (MQ=255)                                                                                                                                    aTCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCt                      >  2:501286/1‑151 (MQ=255)                                                                                                                                    aTCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCt                      >  2:222465/1‑151 (MQ=255)                                                                                                                                    aTCAGAATA‑ATATAACCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGc                       <  1:1030331/150‑1 (MQ=255)                                                                                                                                       aGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCgg                   >  1:323063/1‑151 (MQ=255)                                                                                                                                        ggATC‑ATATGATCAACCCCCTCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTc                     <  2:619628/148‑1 (MQ=255)                                                                                                                                        ggATC‑ATATGACAACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTc                     <  2:487661/148‑1 (MQ=255)                                                                                                                                        gaaTC‑ATATGACAAACCCCCGCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCt                         <  1:126550/142‑1 (MQ=255)                                                                                                                                         ggatattatGACCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGc                       <  1:791775/140‑1 (MQ=255)                                                                                                                                         gaataatatGACAAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCg                    <  1:417408/144‑1 (MQ=255)                                                                                                                                                atGAACACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtgtg            <  1:407075/150‑1 (MQ=255)                                                                                                                                                   aTCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCat         >  1:95954/1‑150 (MQ=255)                                                                                                                                                   aTCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCat         >  1:299384/1‑150 (MQ=255)                                                                                                                                                      aaCCCCCCCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTg     <  2:883169/151‑1 (MQ=255)                                                                                                                                                         acCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGCgg  >  1:166459/1‑151 (MQ=255)                                                                                                                                                         |                                                                                                                                                       TTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGG‑TGTATCTCGAACCGTTAAA‑CCGCTATCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGCGG  >  NC_000913/1376358‑1376655

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