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

CGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGC‑GTACCGGAACC‑GTGGTGTATCTCGAACC‑GTTAAACCGCTATCAGGATC‑ATA‑TGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTG  >  NC_000913/1376355‑1376652                                                                                                                                                            |                                                                                                                                                    cGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGC‑GTACCGGAACC‑GTGGTGTATCTCGAACC‑GTTAAACCGCTATCAGGATC‑ATA‑TGATCAaca                                                                                                                                                     <  2:424838/151‑1 (MQ=255)     tACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGC‑GTACCGGAACC‑GTGGTGTATCTCGAACC‑GTTAAACCGCTATCAGGATC‑ATA‑TGATCAACACCCt                                                                                                                                                 >  1:848681/1‑151 (MQ=255)       ccgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGC‑GTACCGGAACC‑GTGGTGTATCTCGAACC‑GTTAAACCGCTATCAGGATC‑ATA‑TGATCAACACCCTc                                                                                                                                                <  1:398671/150‑1 (MQ=255)        cgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGC‑GTACCGGAACC‑GTGGTGTATCTCGAACC‑GTTAAACCGCTATCAGGATC‑ATA‑TGATCAACACCCTCGc                                                                                                                                              >  2:199499/1‑151 (MQ=255)          ccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGC‑GTACCGGAACC‑GTGGTGTATCTCGAACC‑GTTAAACCGCTATCAGGATC‑ATA‑TGATCAACACCCTCGCCg                                                                                                                                            <  2:252988/151‑1 (MQ=255)               gACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGC‑GTACCGGAACC‑GTGGTGTATCTCGAACC‑GTTAAACCGCTATCAGGATC‑ATA‑TGATCAACACCCTCGCCGATGcc                                                                                                                                       <  2:380945/151‑1 (MQ=255)                    cgccgcGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGC‑GTACCGGAACC‑GTGGTGTATCTCGACCC‑GTTAAACCGCTATCAGGATC‑ATA‑TGATCAACCCCCTCGCCGATGCccgcc                                                                                                                                   <  1:77425/150‑1 (MQ=255)                    cgccgcGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGC‑GTACCGGAACC‑GTGGTGTATCTCGAACC‑GTTAAACCGCTATCAGGATC‑ATA‑TGATCAACACCCTCGCCGATGCccgcc                                                                                                                                   <  2:966215/150‑1 (MQ=255)                              ccTGGACGGCGACCGAAAAAGGGTGAGTGATTCCCGGCGCGTACTGGAACAGGTCGCCGCGC‑GTACCGGACCG‑GTGGTGTATCTCGACCC‑GTTAACCCGCTATCAGGATC‑ATA‑TGATCAACCCCCCCGCCGATGCCCGCCGTTACAtcgt                                                                                                                         <  2:1141958/150‑1 (MQ=255)                                    cGGCGCCCGAAAAAGGGGGAGGGATCCCCTGCGCGTACGGGAAGGGGTCGCCGCGC‑GTCCCGGAACC‑GTGGTGTATCTCGACCC‑GTAAACCCGCTATAAGGATC‑ATA‑TAATCAACCCCCCCCCCGATGCCCGCCGTTACATCGTCGaaaa                                                                                                                   <  2:633012/150‑1 (MQ=255)                                     ggCGACCGAAAAATGGTGAGTGATCCCCTGCGCGTACTGAAACAGGTCGCCGCGC‑GTCCCGGACCC‑GTGGTGTATCTCGAACC‑GTTAAACCGCTATCAGGATC‑ATA‑TGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAAc                                                                                                                  <  2:519563/150‑1 (MQ=255)                                      gCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGC‑GTACCGGAACC‑GTGGTGTATCTCGAACC‑GTTAAACCGCTATCAGGATC‑ATA‑TGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACg                                                                                                                 >  2:353593/1‑150 (MQ=255)                                        gACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGC‑GTACCGGAACC‑GTGGTGTATCTCGAACC‑GTTAAACCGCTATCAGGATC‑ATA‑TGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATc                                                                                                              >  1:239620/1‑151 (MQ=255)                                        gACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGC‑GTACCGGAACC‑GTGGTGTATCTCGAACC‑GTTAAACCGCTATCAGGATC‑ATA‑TGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATc                                                                                                              >  2:722871/1‑151 (MQ=255)                                          ccGAAAATTGGTGAGTGTTCCCCTGCGC‑TTCTGGGACAGGGC‑CCCCGCGGGAACGGGACC‑GGGGTTGTTCCCCAACCTGTTAACCGGCT‑TAAGGTTCAATA‑TGACAACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCt                                                                                                             <  1:43138/150‑1 (MQ=255)                                           cGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGC‑GTACCGGAACC‑GTGGTGTATCTCGAACC‑GTTAAACCGCTATCAGGATC‑ATA‑TGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCtt                                                                                                            >  2:704954/1‑150 (MQ=255)                                             cAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGC‑GTACCGGAACC‑GTGGTGTATCTCGAACC‑GTTAAACCGCTATCAGGATC‑ATA‑TGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTaa                                                                                                          >  1:644090/1‑150 (MQ=255)                                             aaaaaTGGGGAGTGATTCCCGGCGCGTAGTGAAACGGGTCGCCGCGC‑GTCCCGAAACC‑GGGGTATATCTCGAACC‑GTTAACCCGCTATCAGGATT‑ATA‑TAATCAACCCCCGCGCCGGTGCCCGCCGTTACATCGTCGAAAACGATCTTaa                                                                                                          <  2:319910/149‑1 (MQ=255)                                              aaaaTGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGC‑GTACCGGAACC‑GTGGTGTATCTCGAACC‑GTTAAACCGCTATCAGGATC‑ATA‑TGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAAc                                                                                                        >  2:755344/1‑151 (MQ=255)                                                aaTGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGC‑GTACCGGAACC‑GTGGTGTATCTCGAACC‑GTTAAACCGCTATCAGGATC‑ATA‑TGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACAt                                                                                                      >  2:901540/1‑151 (MQ=255)                                                  tGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGC‑GTACCGGAACC‑GTGGTGTATCTCGAACC‑GTTAAACCGCTATCAGGATC‑ATA‑TGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATg                                                                                                     >  2:290455/1‑150 (MQ=255)                                                     tgagtgaTTCCCTGCGCGTACTGGAACAGGTCGCCGCGC‑GTACCGGAACC‑GTGGTGTATCTCGAACC‑GTTAAACCGCTATCAGGATC‑ATA‑TGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTAc                                                                                                  >  2:242312/1‑150 (MQ=255)                                                        gtgaTTCCCTGCGCGTACTGGAACAGGTCGCCGCGC‑GTACCGGAACC‑GTGGTGTATCTCGAACC‑GTTAAACCGCTATCAGGATC‑ATA‑TGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTAc                                                                                                  >  2:1021693/1‑147 (MQ=255)                                                        gtgaTTCCCGGCGCGTACGGGAACAGGTCGCCGCGC‑GTACCGGACCC‑GGGGTGTATCTCGACCC‑GTTAACCCGCTATCAGGATC‑ATA‑TGAACAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAg                                                                                                <  2:695004/149‑1 (MQ=255)                                                            ttCCCGGCGCGTCTTGGAAGAGGCCCCCGGGC‑GT‑CCGGAACCGGTGGTGTATCTCGACCC‑GTAAACCCGCTATCAGGATC‑ATA‑GGAACAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGAtt                                                                                             <  2:305470/148‑1 (MQ=255)                                                             tCCCTGCGCGTACTGGAACAGGTCGCCGCGC‑GTACCGGAACC‑GTGGTGTATCTCGAACC‑GTTAAACCGCTATCAGGATC‑ATA‑TGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCg                                                                                         >  2:947289/1‑151 (MQ=255)                                                               ccTGCGCGTACTGGAACAGGTCGCCGCGC‑GTACCGGACCC‑GTGGTTTATCTCAACCC‑GTTAACCCGCTATCAGGATC‑ATA‑TGATAAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGc                                                                                       <  2:820527/151‑1 (MQ=255)                                                                   cgcgTACTGGAACAGGTCGCCGCGC‑GTACCGGAACC‑GTGGTGTATCTCGAACC‑GTTAAACCGCTATCAGGATC‑ATA‑TGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGAtt                                                                                   >  2:675322/1‑151 (MQ=255)                                                                     cgTACGGGAACGGGTCGCCGCGC‑ATACCGGAACC‑GTGGTGTATCGCGACCC‑GTTAACCCGTTATCAGGATA‑ATA‑TGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGa                                                                                     <  2:538316/147‑1 (MQ=255)                                                                        aCTGGAACAGGTCGCCGCGC‑GTCCCGGACCC‑GTGGTGTATCTCGACCC‑TTTAACCCGCTATCAGGATA‑ATA‑TGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTAt                                                                              <  2:1129712/151‑1 (MQ=255)                                                                         cgggAACGGGTCGCCGCGC‑GTACCGGAACG‑GTGGTTTATCGCGAACC‑TTAAACCCGCTATCAGGATC‑ATA‑TGACAACCCCCCCCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTa                                                                               <  2:736376/147‑1 (MQ=255)                                                                         cTGGAACAGGTCGCCGCGC‑GTACCGGAACC‑GTGGTGTATCTCGAACC‑GTTAAACCGCTATCAGGATC‑ATA‑TGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATc                                                                             >  1:1063763/1‑151 (MQ=255)                                                                           ggAACGGGTCCCCGCCC‑GTCCCGAACCG‑GTGGTTTCTCTCAACCC‑GTAAACCCGCTTTAAGGATC‑ATA‑TAACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATcac                                                                           <  1:746062/151‑1 (MQ=255)                                                                           ggAACAGGTCGCCGCGC‑GTCCGGGAACC‑GTGGTTTATCTCGACCC‑GTTAACCCGTTATAAGGATC‑ATA‑TGACCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATcac                                                                           <  2:1004455/151‑1 (MQ=255)                                                                           ggAACAGGTCGCCGCGC‑GTACCGGAACC‑GTGGTGTATCTCGAACC‑GTTAAACCGCTATCAGGATC‑ATA‑TGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATcac                                                                           >  2:896972/1‑151 (MQ=255)                                                                             aactgcgcgccgcGG‑GTACCGGAACC‑GTGGTGTCTCTCCAACT‑TTAAACCCCCTATGGGGA‑C‑ATATTTACCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTAt                                                                              <  1:823716/139‑1 (MQ=255)                                                                                 ggTCGCCGCGC‑GTACCGGAACC‑GGGGTGTATCTCGAACC‑TTAAAACCGCTATCAGGATA‑ATA‑TGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATg                                                                        <  2:248246/148‑1 (MQ=255)                                                                                      ccgcGC‑GTCCGGAAACC‑GTGGTGTATCTCGAACC‑GTTAACCCTTTATCAGGATA‑ATA‑TGCTCACCCCCCTCGCCATGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATc                                                                  <  1:769438/149‑1 (MQ=255)                                                                                        gcgc‑gTACGGGAACC‑GTGGTTTATCGCAAACC‑TTTAAACCGCAATCGGGATA‑ATT‑TGATCAACCCCCGCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAata                                                                    <  1:497409/145‑1 (MQ=255)                                                                                          gc‑gTCCCGAAACC‑GTGGTGTATCTCGAACC‑GTTAACCCGTTATAAGGATA‑ATA‑TGACCACCCCCCTCCCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaagaa                                                            <  2:945706/151‑1 (MQ=255)                                                                                              tACCGGAACC‑GTGGTGTATCTCGAACC‑GTTAAACCGCTATCAGGATC‑ATA‑TGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                             >  1:869240/1‑147 (MQ=255)                                                                                              gtccgggaccgggggtgTATCTCGACCC‑GTTAACCCGTTATAAGAATA‑ATA‑TAACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAagaaga                                                          <  1:648592/138‑1 (MQ=255)                                                                                                ccGGACCC‑GGGGTGTCTCCCGACCC‑GTAAACCCGCTACCAGAATA‑ATA‑TAATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAagaaga                                                          <  1:814806/148‑1 (MQ=255)                                                                                                                  cccGACCC‑TTTAACCCGCTATCAGGATC‑AAA‑TGATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTg                                       <  1:41037/148‑1 (MQ=255)                                                                                                                     gAACC‑GTTAAACCGCTATCAGGATC‑ATA‑TGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATg                                   >  1:65843/1‑151 (MQ=255)                                                                                                                      aaCC‑TTAAAACCGCTATCAGGATC‑ATA‑TGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCAt                                    <  2:573889/149‑1 (MQ=255)                                                                                                                      aaCC‑GTTAAACCGCTATCAGGATC‑ATA‑TGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATg                                   >  2:878330/1‑150 (MQ=255)                                                                                                                        cc‑GTTAAACCGCTATCAGGATC‑ATA‑TGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGAc                                 >  2:52106/1‑150 (MQ=255)                                                                                                                           gTTAAACCGCTATCAGGATC‑ATA‑TGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAAc                              >  1:250140/1‑151 (MQ=255)                                                                                                                            ttAAACCGCTATCAGGATC‑ATA‑TGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAAcc                             >  1:1039628/1‑151 (MQ=255)                                                                                                                            ttAAACCGCTATCAGGATC‑ATA‑TGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAAcc                             >  1:453584/1‑151 (MQ=255)                                                                                                                             tAAACCGCTATCAGGATC‑ATA‑TGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcg                            >  1:205771/1‑151 (MQ=255)                                                                                                                              aaaCCGCTATCAGGATC‑ATA‑TGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgc                           >  2:167779/1‑151 (MQ=255)                                                                                                                              aaaCCGCTATCAGGATC‑ATA‑TGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgc                           >  2:824777/1‑151 (MQ=255)                                                                                                                              aaaCCGCTATCAGGATC‑ATA‑TGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgc                           >  1:917308/1‑151 (MQ=255)                                                                                                                                 ccGCTATCAGGATC‑ATA‑TGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGAc                        >  2:213196/1‑151 (MQ=255)                                                                                                                                  cGCTATCAGGATC‑ATA‑TGAACAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGAcc                       <  2:24330/151‑1 (MQ=255)                                                                                                                                      acaaGGATA‑ATA‑TGACCACCCCCCCCGCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCt                      <  1:1058580/145‑1 (MQ=255)                                                                                                                                       tCAGGATC‑ATA‑TGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTc                  >  2:590524/1‑151 (MQ=255)                                                                                                                                         aGGATC‑ATA‑TGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCgg                >  1:752377/1‑151 (MQ=255)                                                                                                                                          ggATC‑ATA‑TAAAAACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCg                 <  2:114397/149‑1 (MQ=255)                                                                                                                                            aataata‑tGACCAACCCCCTCGCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTc                  <  1:1021667/143‑1 (MQ=255)                                                                                                                                                    tGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGc        <  2:1146422/150‑1 (MQ=255)                                                                                                                                                    tGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGc        >  1:562643/1‑150 (MQ=255)                                                                                                                                                      aTCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCat      >  1:1122106/1‑150 (MQ=255)                                                                                                                                                         aaCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTg  <  2:189023/151‑1 (MQ=255)                                                                                                                                                         aaCACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTg  >  1:974187/1‑151 (MQ=255)                                                                                                                                                            |                                                                                                                                                    CGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGC‑GTACCGGAACC‑GTGGTGTATCTCGAACC‑GTTAAACCGCTATCAGGATC‑ATA‑TGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTG  >  NC_000913/1376355‑1376652

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