BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,560,455	+G	intergenic (‑2/+1)	glpR ← / ← glpR	pseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon/pseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,560,455	1	.	G	99.2%	406.3 / ‑3.7	118	intergenic (‑2/+1)	glpR/glpR	pseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon/pseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon
Reads supporting (aligned to +/- strand): ref base . (1/0); new base G (43/74); total (44/74)
Fisher's exact test for biased strand distribution p-value = 3.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

Only 100 of 118 total aligned reads displayed.

TATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGA‑ACTGGAAGGCAGCGCCGCACCGCCATGATGG‑C‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTT  >  NC_000913/3560312‑3560598

tatCGATAAACAGCGTCGAGCCTTGGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGA‑ACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGAt                                                                                                                                           <  1:1192260/151‑1 (MQ=255)

tatCAATAAACAGCGTCGAGCCATGGGGGTTTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGA‑ACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGAt                                                                                                                                           <  1:67429/151‑1 (MQ=255)

taccGATAAACAGCGTCGGCCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGA‑ACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGAt                                                                                                                                           <  1:927939/148‑1 (MQ=255)

 ttagATAAGCCGCTTCGGCCCTTTGGGTTTTTGCTCCGCCTCTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCCCGGCGTGTTAACCGA‑ACTGGAAGGCAGCGCCGCCCCGCCATGATGG‑CGCAGGATc                                                                                                                                          <  1:693301/148‑1 (MQ=255)

 atCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGA‑ACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGAt                                                                                                                                           >  2:817486/1‑150 (MQ=255)

  tCGATAACCAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGA‑ACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCCGGATCa                                                                                                                                         >  1:1479739/1‑151 (MQ=255)

  tCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGA‑ACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCa                                                                                                                                         >  2:721939/1‑151 (MQ=255)

  tCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGA‑ACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCa                                                                                                                                         >  2:1356301/1‑151 (MQ=255)

  tCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGA‑ACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCa                                                                                                                                         >  2:1100408/1‑151 (MQ=255)

  tCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGA‑ACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCa                                                                                                                                         >  1:764926/1‑151 (MQ=255)

  tCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGA‑ACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCa                                                                                                                                         >  2:416081/1‑151 (MQ=255)

  tCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGA‑ACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCa                                                                                                                                         >  2:6565/1‑151 (MQ=255)

  tCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCATGCGATCGTGCCACGGCGTGTTAACCGA‑ACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCa                                                                                                                                         >  1:614701/1‑151 (MQ=255)

   cGATAAACAGCGCCGAGCCATTGGGGATTTGCCCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCATGGCGTGTTAACCGA‑ACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCat                                                                                                                                        >  1:335028/1‑150 (MQ=255)

    gATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGA‑ACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAgg                                                                                                                                       >  2:1508965/1‑151 (MQ=255)

             cGTCGGGCCTTGGGGGATTTGCCCCGCCACTTTGCGGGCGAGGCGCTCTTTTTCTTCGGTCGGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGA‑ACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTc                                                                                                                               <  1:671127/150‑1 (MQ=255)

             cGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGA‑ACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGcgc                                                                                                                               >  1:978543/1‑148 (MQ=255)

                 gggCCATTGGGGTTTTGCTCCGCGACTTTGCGGGCGATGCCCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCGCGGGGTGTTAACCGA‑ACGGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGcc                                                                                                                           <  1:1333997/148‑1 (MQ=255)

                     ctttGGGGATTTGCTCCGCCCCTTTGCGGGCGATGGGCTTTTTTTCTGCGGTCTGGGTGGCCTTGCGATGGTGCCACGGCGTGTTACCCGA‑ACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTc                                                                                                                      <  2:635343/149‑1 (MQ=255)

                          gggTTTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGA‑ACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTg                                                                                                                  <  1:1132421/150‑1 (MQ=255)

                          gggATTTGCCCCGCCCCTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCCCGGCGTGTTAACCGA‑ACTGGAAGGCAGCGCCGCCCCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTg                                                                                                                  <  2:484583/150‑1 (MQ=255)

                           ggtttgggcccGCCACTTGGCGGGCGATGCGCTCTTTTTCTTCGGTCGGGGTGGCCTTGCGATCGTGCCCCGGCGTGTTAACCGA‑CCTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGa                                                                                                                 <  1:873330/141‑1 (MQ=255)

                             aTTTGCTCCCCCACTTTGCGGGCGATGCGCTCTTTTTCCTCGGTCGGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGA‑GCTGGAGGGCGGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAgg                                                                                                               <  1:1367163/150‑1 (MQ=255)

                                 gCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGA‑ACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTcgcg                                                                                                          >  1:760769/1‑151 (MQ=255)

                                  cTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGA‑ACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcgg                                                                                                         >  2:1212451/1‑151 (MQ=255)

                                  cTCCGCCACTTTGCGGGCGACGCGCTCTTTTTCTTCGATCTGGGTGGCCTTGCGATCGTGCCCCGGCGTGTTAACCGC‑ACTGGAAGGCAGCGCTGCACCGCCATGCTGG‑CGCAGGATCAGGCTTTGCTCCCCCAGCTCATTGAGGGTgcgg                                                                                                         >  1:619468/1‑151 (MQ=255)

                                       ctaatttGCGGGCGATGCGTTCTTTTTCTTCGGTCTGGGTGGCCTTGGGATCGGGCCACGGCGTGTTAACCAA‑ATTGAAAGGCGGCGCCGCACCGCCATGATGG‑CGCAGGGTCAGGTTTTGCTCCGCCAGCTCAGTGAGGTCGGGGCGAAt                                                                                                    <  1:857678/147‑1 (MQ=255)

                                         ccTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGA‑ACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAg                                                                                                  <  1:174481/150‑1 (MQ=255)

                                         aCTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGA‑ACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAg                                                                                                  <  2:1007988/151‑1 (MQ=255)

                                          cTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGA‑ACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGt                                                                                                 >  2:1213311/1‑151 (MQ=255)

                                                   cGATGCGCTCTTTTTCTTCGGTCTGGGTGGCTTTGCGATCGTGCCACGGCGTGTTAACCGA‑ACTGGAAGGCCGCGCCGCACCGCCATGATGGAG‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCggg                                                                                          >  2:72003/1‑149 (MQ=255)

                                                   cGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGA‑ACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGc                                                                                         <  2:1269600/150‑1 (MQ=255)

                                                   cGATGCGCTCTTTTTCTTCGGTCGGGGTGGCCTTGCGATCGTGCCCCGGCGTGTTAACCGA‑ACTGGAAGGCAGCGCCGCTCCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGc                                                                                         <  2:197608/150‑1 (MQ=255)

                                                   cGATGCCCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGA‑ACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCgggn                                                                                         <  2:1251096/150‑2 (MQ=255)

                                                      tGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGA‑ACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAAGTCGCGGCGAATAGTCTGCGGGCTGAc                                                                                     <  2:1155391/151‑1 (MQ=255)

                                                          ctctTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGA‑ACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGgag                                                                                 >  2:1486910/1‑151 (MQ=255)

                                                            ctTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGA‑ACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCCGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGaa                                                                               >  2:158818/1‑151 (MQ=255)

                                                            ctTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGA‑ACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGaa                                                                               >  1:595196/1‑151 (MQ=255)

                                                            ctTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGA‑ACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGaa                                                                               >  1:1391664/1‑151 (MQ=255)

                                                            ctTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGA‑ACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGaa                                                                               >  1:831032/1‑151 (MQ=255)

                                                            ctTTTTCTTCGGTCTGGGTGGCCTTGAGATCGTGCCACGGCGTGTTAACCGA‑ACTGGAAGGCAGCGCCGGCCCGCCATGATGG‑CGAAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGaa                                                                               >  1:1022944/1‑151 (MQ=255)

                                                             tttttCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGA‑ACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGaaa                                                                              <  2:24136/151‑1 (MQ=255)

                                                                ttcttcGGTCGGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGA‑ACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGc                                                                           <  1:409910/151‑1 (MQ=255)

                                                                  cttcGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGA‑ACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGct                                                                          >  1:90442/1‑150 (MQ=255)

                                                                       gTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGA‑ACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCa                                                                    >  2:366496/1‑151 (MQ=255)

                                                                         cgggggggCCTTGCGATCGTGCCACGGCGTGTTAACCGA‑ACTGGAAGGCGGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAg                                                                   <  1:70383/148‑1 (MQ=255)

                                                                         cTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGA‑ACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAg                                                                   >  2:681921/1‑150 (MQ=255)

                                                                           gggTGGCCTTGCGATCGGGCCACGGCGTGTTAACCGA‑ACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGGTTTGCTCCGCCAGCTCATTGAGGTCGAGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACAAGcg                                                                 >  1:958081/1‑149 (MQ=255)

                                                                            ggTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGA‑ACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGctc                                                                <  1:1441507/150‑1 (MQ=255)

                                                                                 ccTTGCGATGGTGCCACGGCGTGTTAACGGA‑ACTGGAGGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCgg                                                           <  1:1471929/150‑1 (MQ=255)

                                                                                 ccTTGCGATCGTGCCACGGCGTGTTAACCGA‑ACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGt                                                          >  1:329340/1‑151 (MQ=255)

                                                                                 ccTTGCGATCGTGCCACGGCGTGTTAACCGA‑ACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGt                                                          <  1:1341341/151‑1 (MQ=255)

                                                                                   ttGCGATCGTGCCCCGGCGTGTTAACCGA‑ACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTAc                                                        <  2:1078576/151‑1 (MQ=255)

                                                                                   ttGCGATCGTGCCACGGCGTGTTAACCGA‑ACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTAc                                                        <  2:649855/151‑1 (MQ=255)

                                                                                   ttGCGAGAGTGCGACGGCGTGTTAACCGA‑GCTGGAGGGGAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTAc                                                        <  1:51743/151‑1 (MQ=255)

                                                                                      cGATCGTGCCACGGCGTGTTAACCGA‑ACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGa                                                     >  1:74344/1‑151 (MQ=255)

                                                                                      cGATCATGCCACGGCGTGTAACCGGA‑ACGGGAAGGCAGCGCGGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGAGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTg                                                      <  1:967146/150‑1 (MQ=255)

                                                                                        aTCGTGCCACGGCGTGTTAACCGA‑ACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACa                                                   <  1:6554/151‑1 (MQ=255)

                                                                                        aTCGGGCCACGGCGTGTTAACCGA‑ACTGGAGGGCAGCGCCGCGCCGCCATGATGG‑CGCAGGATCAGGGTTTGCTCCGCCAGCTCATTGAGGTCGCGGCAAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACa                                                   <  1:1169835/151‑1 (MQ=255)

                                                                                              ccccGGCGTGTTAACCGA‑ACTGGAAGGCAGCGCCGCACCGCCTTGATGG‑CGCAGGATCACGTTTTGCGCCGCCAGCTCATCGAGGCCTCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCGCGACCAGCTCTTCGGTACtttcatggcc                                              <  1:158771/147‑10 (MQ=255)

                                                                                              ccACGGCGTGTTACCCGA‑CCTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAccc                                             <  1:774382/151‑1 (MQ=255)

                                                                                              ccACGGCGTGTTAACCGA‑ACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAcc                                              <  2:910562/150‑1 (MQ=255)

                                                                                              ccACGACGTGTAAACCGA‑ACTGGAGGGCGACGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGGCGGAGAAATGCTCTACCAGCTCTTCGGTACTGGCATAAccc                                             <  1:785315/151‑1 (MQ=255)

                                                                                              acACGGCGGGTT‑ACCGAGTCTGGAAGGAAGCGCGCCGCCGCTATGATGG‑CGCAGGATCAGGATTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAAGGCTCTCCCAGCTCTTCGGTACTGACATAAcc                                              <  1:848046/149‑1 (MQ=255)

                                                                                                    cGTGTTAACCGA‑ACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCctgctg                                        <  2:1369621/150‑1 (MQ=255)

                                                                                                     gtattAACCGA‑ACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGtt                                      >  2:593121/4‑151 (MQ=255)

                                                                                                      tgtTAACCGA‑ACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGtt                                      <  2:446477/150‑1 (MQ=255)

                                                                                                          tagcGA‑ACTTGACGGGAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAAcc                                 <  1:1002386/148‑1 (MQ=255)

                                                                                                          acccGA‑CCTGGAAGGCGGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCGGACGGAGAAATGCTCTACCAGCGCTTCGGTACTGACATAACCCTGCTGTTTAAcc                                 <  2:657416/149‑1 (MQ=255)

                                                                                                          acccGA‑ACTGGAAGGCAGCGCCGCGCCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAAcc                                 <  1:1355623/149‑1 (MQ=255)

                                                                                                          aaCCGG‑ACGGGAGGCGAGGGCCGCACCGCCATGATGG‑CGAATGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGTCGAATAGTCTGCGGGCTGACAGAGAAATGCTCTACCAGCTCTTCCGTACTGACTTAGCCCTGCTGTTTAAcc                                 <  1:1442043/151‑1 (MQ=255)

                                                                                                          aaCCGA‑CCTGAAAGCCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAAcc                                 <  2:20651/151‑1 (MQ=255)

                                                                                                          aaCCGA‑ACTGGAGGGCAGCGCCGCCCCGCCATTATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAAcc                                 <  2:745913/151‑1 (MQ=255)

                                                                                                          aaCCGA‑ACTGGAGGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAAcc                                 <  1:1253902/151‑1 (MQ=255)

                                                                                                          aaCCGA‑ACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAAcc                                 <  1:545738/151‑1 (MQ=255)

                                                                                                          aaCCGA‑ACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAAcc                                 <  1:738996/151‑1 (MQ=255)

                                                                                                          aaCCGA‑ACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAAcc                                 <  1:415949/151‑1 (MQ=255)

                                                                                                          aaCCGA‑ACTGAAAGGCAGCGCCGCCCCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAAcc                                 <  2:212524/151‑1 (MQ=255)

                                                                                                               a‑aCTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGtt                             >  2:268175/1‑150 (MQ=255)

                                                                                                               a‑aCTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGtt                             <  2:669783/150‑1 (MQ=255)

                                                                                                                  cTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGa                          >  2:84056/1‑151 (MQ=255)

                                                                                                                  cTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGa                          >  2:762069/1‑151 (MQ=255)

                                                                                                                  cTGGAAGGCAGCCCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACCTAACCCTGCTGTTTAACCAgcttga                          >  1:14235/1‑146 (MQ=255)

                                                                                                                      aaGGCGGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGataa                       <  2:780136/150‑1 (MQ=255)

                                                                                                                      aaGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGataat                      <  2:246412/151‑1 (MQ=255)

                                                                                                                         gaaacgacgcACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGATCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATAc                    <  1:555718/143‑1 (MQ=255)

                                                                                                                         gCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATAc                    <  1:217400/150‑1 (MQ=255)

                                                                                                                         gCAGCCGCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATAc                    <  1:1428225/150‑1 (MQ=255)

                                                                                                                         ccAGCGCCGCACCGCCATGATGG‑CGCAGGATCTGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAGCCCTGCTGTTTACCCCGTTCGATGATAc                    <  1:1487764/149‑1 (MQ=255)

                                                                                                                             cgccgcCCCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCg                  <  2:1383930/148‑1 (MQ=255)

                                                                                                                             cgccgcACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCg                  <  2:614995/148‑1 (MQ=255)

                                                                                                                               ccgcACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTtgtg             <  2:737416/151‑1 (MQ=255)

                                                                                                                                cgcACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGa            <  2:91583/151‑1 (MQ=255)

                                                                                                                                cgcACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGa            <  1:882287/151‑1 (MQ=255)

                                                                                                                                cgcACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGa            <  1:224339/151‑1 (MQ=255)

                                                                                                                                cgcACCGCCATGATGG‑CGCAGGATCAGGTTTTGCGCCGCCATCTCATTGAGGTCGCGGCGAATAGTCGGCGGTCTGACGGAGAAATGCTCTACCAGCTCTGCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGGTTTGa            <  1:320017/151‑1 (MQ=255)

                                                                                                                                  ccccGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACg          <  1:470964/149‑1 (MQ=255)

                                                                                                                                     cGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTtg       <  2:891769/151‑1 (MQ=255)

                                                                                                                                      gCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTtgt      >  1:193200/1‑151 (MQ=255)

                                                                                                                                          tgatgG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGGCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGATGtgttt  >  2:248740/1‑151 (MQ=255)

TATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGA‑ACTGGAAGGCAGCGCCGCACCGCCATGATGG‑C‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTT  >  NC_000913/3560312‑3560598

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