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.9%	115.0 / 14.6	62	T158P (ACC→CCC)	ycjR	putative TIM alpha/beta barrel enzyme
Reads supporting (aligned to +/- strand):  ref base A (29/7);  new base C (0/26);  total (29/33)
Fisher's exact test for biased strand distribution p-value = 2.41e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGA‑TCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCAT  >  NC_000913/1376355‑1376648                                                                                                                                                          |                                                                                                                                                cGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGA‑TCATATGATCAaca                                                                                                                                                 <  2:530046/151‑1 (MQ=255)    ttACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGA‑TCATATGATCAACAcc                                                                                                                                               <  1:419583/150‑1 (MQ=255)     tACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGA‑TCATATGATCAACAcc                                                                                                                                               <  1:156875/149‑1 (MQ=255)       ccgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGA‑TCATATGATCAACACCCTCg                                                                                                                                           <  2:736057/151‑1 (MQ=255)            gATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGA‑TCATATGATCAACACCCTCGCCGAt                                                                                                                                      >  2:35311/1‑151 (MQ=255)              tGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGA‑TCATATGATCAACACCCTCGCCGATGc                                                                                                                                    >  1:110329/1‑151 (MQ=255)                 ccTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGA‑TCATATGATCAACACCCTCGCCGATGCccg                                                                                                                                 >  1:274611/1‑151 (MQ=255)                  cTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGA‑TCATATGATCAACACCCTCGCCGATGCccgc                                                                                                                                >  2:1050705/1‑151 (MQ=255)                    cgccgcGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGA‑TCATATGATCAACACCCTCGCCGATGCccgcc                                                                                                                               <  1:69104/150‑1 (MQ=255)                    cgccgcGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGA‑TCATATGATCAACACCCTCGCCGATGCccgcc                                                                                                                               <  2:1000449/150‑1 (MQ=255)                     gccgcGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGA‑TCATATGATCAACACCCTCGCCGATGCccgcc                                                                                                                               <  1:800215/149‑1 (MQ=255)                     gccgcGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGA‑TCATATGATCAACACCCTCGCCGATGCCCGCCGt                                                                                                                             >  2:174547/1‑151 (MQ=255)                          gTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGA‑TCATATGATCAACACCCTCGCCGATGCCCGCCGTTACAt                                                                                                                        >  1:557189/1‑151 (MQ=255)                           tAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGA‑TCATATGATCAACACCCTCGCCGATGCCCGCCGTTACAtc                                                                                                                       >  1:537558/1‑151 (MQ=255)                                 ggACGGCGCCCGAAAAAGGGTGGGTTTTTCCCGGCGCGTACGGAAACGGGCCGCCGCGCGTCCCGGAC‑CCGTGTTTTTTCTCGACCCTTAAAC‑CCGCTATCAGGT‑TCATAGGATCACCCCCCCCGCCGATGCCCGCCGTTACAtcgt                                                                                                                     <  1:1104714/147‑1 (MQ=255)                                 ggACGGCGACCGCAAAAGGGTGAGTGATCCCGTGGGGGTCCGGGAACAGGTCGCCGCGCGAACGGGAC‑CGGGGGTTTTTCTCAACCCGTTAAC‑CCTTTAACGGGA‑TAATATGACAAACCCCCCCGCCGATGCCCGCCGTTACAtcgtc                                                                                                                    <  1:477566/148‑1 (MQ=255)                                        gACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGA‑TCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATc                                                                                                          >  1:538492/1‑151 (MQ=255)                                          ccGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGA‑TCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCtt                                                                                                        >  2:331685/1‑151 (MQ=255)                                              aaagTGGTGGGTGGTCCCTTGCGCGTGCGGGAACAGGTGGCCGCGCGTCCCGGAC‑GTGGGGTTTATCGCGACCCTTTATC‑CCGCAACAGGGA‑TAATATGTGCACCACCCGCGTCGATGCCCGCCGTTACATCGTCGAAAACGATCtt                                                                                                        <  1:1113275/147‑1 (MQ=255)                                               aaaGGGTGGGTTTTTCCCGGCGCTTCCGGAAAAAGGTCGCCGCGCGTCCGGGAC‑CGGGGGTGTATCCCGACCCTTTAAC‑CCGCTATCGGGA‑TCATATTAACAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACa                                                                                                   <  2:208781/151‑1 (MQ=255)                                               aaaGGGTGGGTGTTCCCCGGCGCGTCCGGAAACGGGCCGCCGCGCGTCCGGGAC‑CCGTGGTGTCCCCCGACCCTTAAAC‑CCGCTACCAGAA‑TAATATGAACACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACa                                                                                                   <  1:538228/151‑1 (MQ=255)                                                     tgagtgaTTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGA‑TCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTAc                                                                                              >  2:384763/1‑150 (MQ=255)                                                       agtgaTTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGA‑TCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTAc                                                                                              >  1:541132/1‑148 (MQ=255)                                                          gtttCCCTGCGCGTAGTGGACCGGGTCGCCGCGCGTACGGGAA‑CCGTGGTCTATCTCAAACGTTTAAC‑CCGTTATCAGAA‑CAATATCACCACCCCCCCCGCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGAtt                                                                                         <  1:140253/148‑1 (MQ=255)                                                          gaTTCCCTGCGCTTACTGGAACAGGTCGCCGCGCGTCCCGGAA‑CCGTGGTGTATCTCGAACCTTTAAA‑CCGTTATCAGAT‑TAATATGATCAACCCCCTCCCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGAtt                                                                                         <  2:887050/150‑1 (MQ=255)                                                             tCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGA‑TCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCg                                                                                     >  1:714398/1‑151 (MQ=255)                                                              cccTGCGCGTACTGGAACAGGTCCCCGCGGGTCCGGGAA‑CCGTGGTGTATCTCGACCCTTTAAC‑CCTCTAACAGGA‑TCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCgg                                                                                    <  2:2255/151‑1 (MQ=255)                                                                   ctcttACGGGAAAGGGCCGCCGCGTGT‑CCGGAACCCGGGGTGTATCTCAACCCTTTAAC‑CCGCTATAAGGA‑TAAAATAATAACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCg                                                                                  <  2:864036/144‑1 (MQ=255)                                                                      gTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGA‑TCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCt                                                                            >  1:71736/1‑151 (MQ=255)                                                                       tACTGGAACGGGCCGCCGCGCGTCCCGAAC‑CCGTGGTGTATCTCGAACCTTAAAC‑CCGCTAACAGAA‑TAATATGACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGa                                                                                 <  1:261820/145‑1 (MQ=255)                                                                           ggAACGGGCCGCCGCGCGTCCGTGAC‑CCGTGGTGTGCCTCGAACCTTTAAC‑CCCCTACTGGCA‑TAAAATGACCACCCCCCCCGCCGAGTCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATca                                                                        <  1:810236/150‑1 (MQ=255)                                                                           ggAACAGGTCGCCGCGCGTACCGGAC‑CCGTGGTGTATCTC‑AACCCTTAAACCCGCTATCAGGA‑TCATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATc                                                                         <  2:915621/149‑1 (MQ=255)                                                                           gaaaCGGGTCGCCGCGCGT‑CCGGAACCCGTGGTGTATCTCGACCCGTAAAC‑CCGCTATCAGGA‑TCATATGACCAACCCCCCCGCGGATGCCCGCCGTTACATCGTCGAAAACGATTTTAAACATGTACAGATTATCGGCGATTTCTATcac                                                                       <  2:302614/149‑1 (MQ=255)                                                                             aaCAGGCCGCCGCGCGTCCGGGAC‑CCGTGGTGTATCCCGACCCGTAAAC‑CCGCTACCAGTA‑TAATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATcaca                                                                      <  1:1170142/150‑1 (MQ=255)                                                                                aGGTCGCCGCGCGTACCGGAA‑CCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGA‑TCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGa                                                                   >  2:543213/1‑150 (MQ=255)                                                                                 ggTCGCCGCGCGTCCCGGAC‑CCGTGGTGTATCTCGACCCTTTAAC‑CCGCTATCAGGA‑TCATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATg                                                                    <  2:313690/148‑1 (MQ=255)                                                                                 ggTCGCCGCGCGTACCGGAC‑CCGTGGTGTATCTCGACCCTTTAAC‑CCGTTATCAGGA‑TCATATGATCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACAt                                                                     <  2:223308/147‑1 (MQ=255)                                                                                  gTCGCCGCGCGTACCGGAA‑CCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGA‑TCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAat                                                                 >  2:439821/1‑150 (MQ=255)                                                                                    cgccgcGCGTACCGGAA‑CCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGA‑TCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATc                                                              >  2:1138636/1‑151 (MQ=255)                                                                                            gTACCGGAA‑CCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGA‑TCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                         >  2:50756/1‑148 (MQ=255)                                                                                             gtcccgaacccGTGGTGTATCTCGAACCTTAAAC‑CCGCTATCAGGA‑TCATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGaagaag                                                       <  2:236302/141‑1 (MQ=255)                                                                                               ccGGAC‑CGGGGGTTTATCTCGACCCGTAAAC‑CCGCTATAAGAA‑TCATATGATAACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATaa                                                   <  2:372904/151‑1 (MQ=255)                                                                                               ccGGAC‑CCGTGGTTTCTCCCAACCCTTTAAC‑CCGTTATCGGAA‑TAATATAACAACCCCCCTCGCCGCGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaag                                                          <  1:403901/144‑1 (MQ=255)                                                                                               ccGAAA‑CGGTGGTGTATCCCGACCCTTAAAA‑CCGTTAGCAGAA‑ACATAGAACAAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGaagaag                                                       <  1:670462/147‑1 (MQ=255)                                                                                                 ggAC‑CCGTGTTGTATCTCGACCCTTAAAC‑CCGTTATAAGAA‑TCATATAACCACCCCCCCCGCCGTGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAAcc                                                 <  1:1058686/151‑1 (MQ=255)                                                                                                          ggtgttccgcaaaCTGTTAAA‑CCGTTATCAGGA‑TCATATGATAAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgc                                             <  2:105232/136‑1 (MQ=255)                                                                                                                  ctcaaCCCTTAAACCCGCTAT‑AGGATTCATATAACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAg                                         <  1:1050698/141‑1 (MQ=255)                                                                                                                    gAACCGTTAAA‑CCGCTATCAGGA‑TCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATg                               >  1:944426/1‑151 (MQ=255)                                                                                                                       ccGTTAAC‑CCGCTATCAGGA‑TAATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACa                            <  2:962578/151‑1 (MQ=255)                                                                                                                        cGTTAAA‑CCGCTATCAGGA‑TCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACa                            >  1:111225/1‑150 (MQ=255)                                                                                                                             aa‑CCGCTATCAGGA‑TCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgc                       >  1:449204/1‑150 (MQ=255)                                                                                                                                ccGCTATCAGGA‑TCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGAc                    >  2:785843/1‑151 (MQ=255)                                                                                                                                 cGCTATCAGGA‑TCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGAcc                   >  2:48773/1‑151 (MQ=255)                                                                                                                                  gCTATCAGGA‑TCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCt                  >  2:539379/1‑151 (MQ=255)                                                                                                                                   cTATCAGGA‑TCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTg                 >  2:643640/1‑151 (MQ=255)                                                                                                                                         ggA‑TCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGt           >  2:1115810/1‑151 (MQ=255)                                                                                                                                         ggA‑TAATATGATAACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCt                  <  2:975968/144‑1 (MQ=255)                                                                                                                                           a‑TCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTc          >  1:985989/1‑150 (MQ=255)                                                                                                                                              cATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAt        >  2:906695/1‑150 (MQ=255)                                                                                                                                                tatGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtgtg     >  2:1040451/1‑151 (MQ=255)                                                                                                                                                 atGACCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTc          <  2:739526/145‑1 (MQ=255)                                                                                                                                                  tGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGc    >  2:920206/1‑150 (MQ=255)                                                                                                                                                     caaaCCCCCCCGCTGGTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCat  <  2:536700/147‑1 (MQ=255)                                                                                                                                                          |                                                                                                                                                CGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CCGTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGA‑TCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCAT  >  NC_000913/1376355‑1376648

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