Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,376,505	0	A	C	41.2%	146.3 / 11.8	81	T158P (ACC→CCC)	ycjR	putative TIM alpha/beta barrel enzyme
Reads supporting (aligned to +/- strand):  ref base A (36/11);  new base C (0/33);  total (36/45)
Fisher's exact test for biased strand distribution p-value = 2.52e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAAT‑GGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTA‑AA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGCGG  >  NC_000913/1376355‑1376655                                                                                                                                                          |                                                                                                                                                       cGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAAT‑GGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTA‑AA‑CCGCTATCAGGATCATATGATCAaca                                                                                                                                                        <  2:1394821/151‑1 (MQ=255)  gCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAAT‑GGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTA‑AA‑CCGCTATCAGGATCATATGATCAacac                                                                                                                                                       >  1:1287789/1‑151 (MQ=255)    ttACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAAT‑GGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTA‑AA‑CCGCTATCAGGATCATATGATCAACAccc                                                                                                                                                     >  2:1096994/1‑151 (MQ=255)     tACCGCCGATGACCTGGCCGCGTAGCCTGGACGGCGACCGCAAAAT‑GGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTA‑AA‑CCGCTATCAGGATCATATGATCAACACCCt                                                                                                                                                    >  2:922649/1‑151 (MQ=255)      aCCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAAT‑GGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTA‑AA‑CCGCTATCAGGATCATATGATCAACACCCt                                                                                                                                                    <  1:1177207/150‑1 (MQ=255)      aCCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAAT‑GGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTA‑AA‑CCGCTATCAGGATCATATGATCAACACCCTc                                                                                                                                                   <  1:289225/151‑1 (MQ=255)       ccgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAAT‑GGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTA‑AA‑CCGCTATCAGGATCATATGATCAACACCCTCg                                                                                                                                                  <  2:961350/151‑1 (MQ=255)            gATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAAT‑GGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTA‑AA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGAt                                                                                                                                             <  2:1031112/151‑1 (MQ=255)             aTGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAAT‑GGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTA‑AA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATg                                                                                                                                            >  1:1302284/1‑151 (MQ=255)             aTGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAAT‑GGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTA‑AA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATg                                                                                                                                            <  1:1124119/151‑1 (MQ=255)                aCCTCGCCGCGTAGCCTGGACGGCGACCGCAAAAT‑GGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTA‑AA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGccc                                                                                                                                         <  2:418098/151‑1 (MQ=255)                aCCTCGCCGCGTAGCCTGGACGGCGACCGCAAAAT‑GGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTA‑AA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGccc                                                                                                                                         >  1:754817/1‑151 (MQ=255)                    cgccgcGTAGCCTGGACGGCGACCGCAAAAT‑GGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTA‑AA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCccgccg                                                                                                                                     >  2:511046/1‑151 (MQ=255)                     gccgcGTAGCCTGGACGGCGACCGCAAAAT‑GGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTA‑AA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGt                                                                                                                                    >  2:1073467/1‑151 (MQ=255)                          gtggccgggcgggCGCCCGAAAAAT‑GGTGAGTGTTTCCCTGCGCTTACGGAAACGGGCCGCCGCGCGTCCGGGACCGGGGGTGTATCTCGACCCGTTA‑AC‑CCGCTATCAGGATAATATGATAAACCCCCTCGCCGATGCCCGCCGTTACAt                                                                                                                               <  1:338717/140‑1 (MQ=255)                          gTAGCCTGGACGGCGACCGCAAAAT‑GGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTA‑AA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACAt                                                                                                                               >  1:67781/1‑151 (MQ=255)                              ccgggAGGGCGACCGAAAAAG‑GGTGGGTGATCCCCGGCGGTTACTGAAAAAGGTCGCCGCGCGTACGGGACCGGGGGTTTATCTCGACCCGTAA‑AC‑CCGCTATCAGGATCATATGACAACCCCCCCCGCCGATGCCCGCCGTTACAtcgtc                                                                                                                           <  2:241401/148‑1 (MQ=255)                              ccTGGACGGCGCCCGCAAAAT‑GGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGACCCTTTA‑AC‑CCGCTATCAGGATCATATGATAAACCCCCCCGCCGATGCCCGCCGTTACAtcgtc                                                                                                                           <  1:249931/151‑1 (MQ=255)                                    cGGCGCCCGAAAAAG‑GGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCCCCGCGCGTCCCGGACCCGGGGTGTATCTCGACCCTTTA‑AC‑CCGCTATCAGGATCATATGATCACCCCCCCCCCCGATGCCCGCCGTTACATCGTCGAAAAc                                                                                                                     <  1:458740/151‑1 (MQ=255)                                           ccgaaaaTGGGTGAGTGTTTCCCTGCGCTTCCGGGAACAGCCCCCCGCGCGTCCCGGACCCGGGGTGTACCTC‑AACCCTTA‑AACCCGCTACCAGGATAATATGACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCtt                                                                                                               <  1:457680/148‑1 (MQ=255)                                             aaaaaG‑GGTGAGTTATCCCCTGCGCGTACTGAAACAGGTCGCCGCGCGTCCCGGACCCGTGGTTTATCTCGACCCTTTA‑AC‑CCGCTATCAGGATAATATGATAAACCCCCGCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTaaa                                                                                                            <  1:339834/150‑1 (MQ=255)                                              aaaaT‑GGTGAGTTATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCTTACCGGACCGGGGGTGTTCCCCGAACCGTAA‑AC‑CCGCTACCAGGATAAAATGATCAACACCCTCGCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTaaan                                                                                                           <  1:157723/151‑2 (MQ=255)                                              aaaaG‑GGGGAGTGATTCCCGGCGCGTACGGGAAGGGGTCGCCGCGCTTCCGGGAGCTTTGGTGTATCACGAGCGGTTA‑AA‑CCGCTACAGGGATTATATAACAACCCCCCTCGCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAAc                                                                                                           <  2:985274/151‑1 (MQ=255)                                                 aT‑GGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTA‑AA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATg                                                                                                        >  1:1419426/1‑151 (MQ=255)                                                   gggTGTGTTATCCCCTGCGGTTTCGGAAACGGGTCGCCGCGCGTCCCGGACCTGGGGTTTATCTCAAACCTTTA‑AC‑CCGCTACCAGCATCATATGACAACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGt                                                                                                       <  2:500989/150‑1 (MQ=255)                                                        agtgaTTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTA‑AA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTAc                                                                                                     >  2:698915/1‑148 (MQ=255)                                                           gaTCCCCTGCGCTTACTGAAACAGGTCGCCGCGCGTCCCGGACCCGTGGTGTATCCCGACCCGTTA‑AC‑CCGCTATCAGGATAATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGAtt                                                                                                <  1:55531/150‑1 (MQ=255)                                                            atccccGGCGCGTACGGAAACAGGTCCCCGCGGGTACCGAACCCGTGGTGTATCTCGAACCGTTA‑AC‑CCGCAATCAGGATAATATGATCAACCCCCTCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTAt                                                                                              <  2:961916/148‑1 (MQ=255)                                                                ccggcgcGTGTTGGAACAGGTCGCCCCGCGTCCCGGACCCGTGGTTTATCTCGACCCTTTA‑AC‑CCGTTATCAGGATCATATGATCAACCCCCGCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGc                                                                                          <  1:109764/148‑1 (MQ=255)                                                                 cTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTA‑AA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCg                                                                                         >  1:713440/1‑151 (MQ=255)                                                                 cTGCGAGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTA‑AA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCg                                                                                         >  1:662596/1‑151 (MQ=255)                                                                          cTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTA‑AA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATc                                                                                >  1:1408996/1‑151 (MQ=255)                                                                            ggAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTA‑AA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATcac                                                                              >  1:873091/1‑151 (MQ=255)                                                                             gAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTA‑AA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATcac                                                                              >  1:1194224/1‑150 (MQ=255)                                                                               agggggccccggGCGGACGGGACCCGTGTTTTCTCTC‑AACCTTTACAC‑CCGCTATCAGGATCATATGAACAAACCCCTCGCCGAGGTCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATg                                                                           <  2:985816/143‑1 (MQ=255)                                                                                cAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTA‑AA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTCTCACATg                                                                           >  1:1100463/1‑150 (MQ=255)                                                                                cAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTA‑AA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATg                                                                           >  1:472730/1‑150 (MQ=255)                                                                                 gggTCGCCGCGCGTACCGGACCCGTGGTGTATCTCGACCCTTTA‑AC‑CCGCTATCAGGATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGaa                                                                         <  2:1151170/150‑1 (MQ=255)                                                                                  ggTCCCCGCGCGTCCCGAACCCGGGGTGTTTCTCGAACCTTAA‑AC‑CCGCTATCAGGATAATATAATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAat                                                                        <  2:1243999/151‑1 (MQ=255)                                                                                      gccgcGCGTACGGGAGCGGGGGTGTATCGCGAACCGTAA‑AC‑CCGCTATCAGAATAATATAATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCg                                                                    <  1:1092023/151‑1 (MQ=255)                                                                                      gccgcGCGTACGGGAACCGTGGTGTATCTCGAACCGTAA‑AC‑CCGTTATCAGGATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCg                                                                    <  2:933597/151‑1 (MQ=255)                                                                                       ccgcGCGTACCGGAACGGTGTTGTATCTC‑AACCCTTA‑AACCCGTTATAAGGATCATATGACAACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCga                                                                   <  2:187021/151‑1 (MQ=255)                                                                                          cgcgTACCGGAACCGTGTTGTCTCCCAACCCTTAA‑AA‑CCGTTATCAGATTAAAATAATCACCACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                                <  1:397265/151‑1 (MQ=255)                                                                                             gTACCGGAACCGTGGTGTATCTCGAACCGTTA‑AA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                                >  1:1262798/1‑148 (MQ=255)                                                                                               cccGGAACCGGGGTGTATCTCGAACCGTTA‑AA‑CCGCTATCAGGATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATa                                                           <  2:588858/150‑1 (MQ=255)                                                                                               aCCGGAACGGTGGTGTATCCCGACCCGTAA‑AC‑CCGCTATCAGGATAATATGAACAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGAt                                                            <  1:1141460/150‑1 (MQ=255)                                                                                               aCCGGAACCGTGGTGTATCTCGAACCGTTA‑AA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGAt                                                            >  2:556906/1‑150 (MQ=255)                                                                                                   gAACGGTGGTGTATCTCGAACCGTTA‑AC‑CCGCTATCAGGATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCt                                                       <  2:29239/151‑1 (MQ=255)                                                                                                     cccGTGGTTTACCTCGACCCGTAA‑AA‑CCGTTATAAGGAACATTTGATAAACCCCCCCGCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTgg                                                     <  1:300746/150‑1 (MQ=255)                                                                                                        gtggtgTATCTCGAACCGTTA‑AA‑CCGCTATCAGGATCATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcgc                                                  <  2:178143/151‑1 (MQ=255)                                                                                                          ggtgTATCTC‑AACCCTTA‑AACCCGCTATCAGGATCATATGATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAg                                                <  2:1287788/151‑1 (MQ=255)                                                                                                             ttacctcGACCCTTAA‑AC‑CCGCTATAAGAATCATATGACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgc                                              <  1:1231854/146‑1 (MQ=255)                                                                                                             gtATCTCGAACCGTTA‑AA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgcg                                             >  1:13997/1‑151 (MQ=255)                                                                                                             gtATCTCGAACCGTTA‑AA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgcg                                             >  2:1276472/1‑151 (MQ=255)                                                                                                             gtATCTCGAACCGTAA‑AC‑CCGCTATAAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgc                                              <  1:1227658/150‑1 (MQ=255)                                                                                                               atctcaaccttttA‑AC‑CCGCTATCGGGTTCATATGACAACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgcgc                                            <  2:169670/139‑1 (MQ=255)                                                                                                                    gAACCTTTA‑AC‑CCGCTATCGGGATAATATAATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATg                                      <  1:370605/151‑1 (MQ=255)                                                                                                                    gAACCGTTA‑AA‑CCGCTATCAGGATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATg                                      <  2:973087/151‑1 (MQ=255)                                                                                                                     aaCCGTTA‑AA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATg                                      >  1:1043460/1‑150 (MQ=255)                                                                                                                        cGTTA‑AA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACa                                   >  1:947017/1‑150 (MQ=255)                                                                                                                           tA‑AA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcg                               >  2:200831/1‑151 (MQ=255)                                                                                                                            a‑aa‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgc                              >  2:876475/1‑151 (MQ=255)                                                                                                                            a‑aa‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgc                              >  1:602590/1‑151 (MQ=255)                                                                                                                             aaacccGCTAACAGGATCATATGACAAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAAcc                                <  2:754810/146‑1 (MQ=255)                                                                                                                              aa‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgcg                             >  1:610003/1‑151 (MQ=255)                                                                                                                                 ccacTAACAGAACCATATAATCACCCCCCCCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGAc                           <  1:734655/148‑1 (MQ=255)                                                                                                                                   gCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCt                         >  1:1269395/1‑151 (MQ=255)                                                                                                                                   ccaagcAGGATCTTATGATCACCGCCCCCGCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCt                         <  1:1095909/146‑1 (MQ=255)                                                                                                                                    gtaacAGGATTATTTAATCAACACCCAGGCCGTTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGAcc                          <  2:166747/145‑1 (MQ=255)                                                                                                                                      aTCAGGATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGc                       <  2:1192817/150‑1 (MQ=255)                                                                                                                                      aTCAGGATAATATGACAACCCCCCCCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCt                      <  1:1081433/151‑1 (MQ=255)                                                                                                                                      aTCAGAATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCt                      >  1:1448759/1‑151 (MQ=255)                                                                                                                                        aaGGATCATATGATCACCCCCCTCGCCGTTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTc                     <  2:1069957/149‑1 (MQ=255)                                                                                                                                           gATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTc                 >  1:1300155/1‑151 (MQ=255)                                                                                                                                           gATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTc                 >  1:159159/1‑151 (MQ=255)                                                                                                                                           aaTCATAGGACAAACACCCTCGCAGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTc                 <  2:1065265/150‑1 (MQ=255)                                                                                                                                               atatGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtgt             >  2:376680/1‑151 (MQ=255)                                                                                                                                               atatGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtgt             >  2:1280707/1‑151 (MQ=255)                                                                                                                                                 atGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtgtg            <  2:67777/150‑1 (MQ=255)                                                                                                                                                    aTCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCat         >  1:204778/1‑150 (MQ=255)                                                                                                                                                     tCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCatat       >  1:434653/1‑151 (MQ=255)                                                                                                                                                      cAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATAtt      >  2:1375461/1‑151 (MQ=255)                                                                                                                                                       aaCACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTg     >  1:1098765/1‑151 (MQ=255)                                                                                                                                                        acacCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGc    >  1:274978/1‑151 (MQ=255)                                                                                                                                                          acCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGCgg  >  1:820397/1‑151 (MQ=255)                                                                                                                                                          acCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGCgg  >  2:120270/1‑151 (MQ=255)                                                                                                                                                          |                                                                                                                                                       CGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAAT‑GGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTA‑AA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGCGG  >  NC_000913/1376355‑1376655

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