BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	459,910	A→G	7.0%	intergenic (‑56/+23)	pqqB ← / ← pqqA	coenzyme PQQ synthesis protein B/coenzyme PQQ synthesis protein A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	459,910	0	A	G	7.0%	335.2 / 7.6	112	intergenic (‑56/+23)	pqqB/pqqA	coenzyme PQQ synthesis protein B/coenzyme PQQ synthesis protein A
Reads supporting (aligned to +/- strand): ref base A (49/55); new base G (5/3); total (54/58)
Fisher's exact test for biased strand distribution p-value = 4.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 112 total aligned reads displayed.

GTGCCATCACGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGCCGATACGCAGGTCAGTGTA‑‑TGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGCTACTACTTG  >  NC_002947/459763‑460051

gTGCCATCACGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAg                                                                                                                                                  <  5:40712/149‑1 (MQ=255)

 tGCCATCACGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAgg                                                                                                                                                 <  3:271726/149‑1 (MQ=255)

    cATCACGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGt                                                                                                                                              <  2:289269/149‑1 (MQ=255)

    cATCACGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGt                                                                                                                                              <  2:284487/149‑1 (MQ=255)

    cATCACGGTACCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGt                                                                                                                                              <  2:102330/149‑1 (MQ=255)

     aTCACGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGt                                                                                                                                              >  8:320672/1‑148 (MQ=255)

      tCACGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTg                                                                                                                                             >  2:115319/1‑148 (MQ=255)

        aCGGTAGCGCTTGCAGTAGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGaaa                                                                                                                                          <  5:308432/149‑1 (MQ=255)

        aCGGTAGCCCTTTCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGATCCTGGATGTACATGGCCACTCTCTTTGCAGAAACC‑GGTTCGTAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCTCGGCCTGCCGGGGCGttttc                                                                                                                                          <  5:268324/149‑5 (MQ=255)

        aCGGTAGCCCTTGCAGTTCACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGCATGTACATGGCCACTCTCTATGCAGAAACC‑GGAGCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑TAAACGCCCCGGCAGGCCGAGGCGTGaaa                                                                                                                                          <  1:215295/149‑1 (MQ=255)

          ggTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAAcc                                                                                                                                        >  7:254686/1‑149 (MQ=255)

          ggTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAAcc                                                                                                                                        <  1:281262/149‑1 (MQ=255)

           gTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAAcc                                                                                                                                        >  2:174826/1‑148 (MQ=255)

               cccTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAg                                                                                                                                   >  8:102516/1‑149 (MQ=255)

                    gCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTg                                                                                                                              <  7:102516/148‑1 (MQ=255)

                        ttGACGCAGTTGCAGATCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTa                                                                                                                         <  5:93812/149‑1 (MQ=255)

                          gACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAg                                                                                                                        <  8:73697/148‑1 (MQ=255)

                                gttgcagttCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCTGAACCC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCg                                                                                                                      >  4:10138/1‑144 (MQ=255)

                                gttgcagttCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCTGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCg                                                                                                                      <  3:10138/144‑1 (MQ=255)

                                gttgcagttCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAg                                                                                                                 <  1:371982/149‑1 (MQ=255)

                                    cagttCCAGTGCGGGAACCCGCCACGGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGaa                                                                                                             <  4:49750/149‑1 (MQ=255)

                                    cagttCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGaa                                                                                                             <  4:337767/149‑1 (MQ=255)

                                     agttCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAg                                                                                                            >  1:243099/1‑149 (MQ=255)

                                       ttCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTa                                                                                                          >  2:229996/1‑149 (MQ=255)

                                           aCTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑g                                                                                                     >  2:375030/1‑149 (MQ=255)

                                           aCTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑g                                                                                                     >  1:158285/1‑149 (MQ=255)

                                            cTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑gg                                                                                                    >  2:172941/1‑149 (MQ=255)

                                            cTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CAGGCTTAgg                                                                                                                       >  7:101719/1‑130 (MQ=255)

                                            cTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CAGGCTTAgg                                                                                                                       <  8:101719/131‑2 (MQ=255)

                                                 ggAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGt                                                                                                   >  8:265250/1‑145 (MQ=255)

                                                 ggAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACt                                                                                               <  1:229996/149‑1 (MQ=255)

                                                  gAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTGCATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAAccc                                                                                                                                       <  8:290680/110‑1 (MQ=255)

                                                     cccGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCg                                                                                            >  7:24624/1‑148 (MQ=255)

                                                      ccGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCg                                                                                            >  4:44525/1‑147 (MQ=255)

                                                      ccGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGaa                                                                                          <  1:115319/149‑1 (MQ=255)

                                                        gCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGCCTTCGAAGc                                                                                        >  1:218419/1‑149 (MQ=255)

                                                        gCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGc                                                                                        <  1:207900/149‑1 (MQ=255)

                                                          cACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGcc                                                                                       >  5:189176/1‑148 (MQ=255)

                                                          cACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGcc                                                                                       >  2:329459/1‑148 (MQ=255)

                                                           aCCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGCCGa                                                                                     >  2:43624/1‑149 (MQ=255)

                                                           aCCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTATATGCAGAAACC‑GGATCTAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTCGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGCCGa                                                                                     <  2:187503/149‑1 (MQ=255)

                                                             cggcggcggAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACtt                                                                                              <  8:4994/138‑1 (MQ=255)

                                                             cggcggcggAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACtt                                                                                              >  7:4994/1‑138 (MQ=255)

                                                               gcggcggAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGCCGATACg                                                                                 >  2:226484/1‑149 (MQ=255)

                                                                cggcggAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGCCGATACGc                                                                                <  5:369150/149‑1 (MQ=255)

                                                                  gcggGGCCGAGGACCTGGATGTACATGGGCACGCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAAAGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCGTTGCGGGTAGCGAAGTACAT‑GGTGACTTGGAAGCCGATTCGc                                                                                >  6:369150/1‑147 (MQ=255)

                                                                    ggAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGCCGATACGCAGGt                                                                            <  1:329459/149‑1 (MQ=255)

                                                                      aGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAACTACGCTTCCCCTTCTACAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGcc                                                                                       <  1:121900/149‑1 (MQ=255)

                                                                      aGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGCCGATACGCAGGTCa                                                                          <  2:22978/149‑1 (MQ=255)

                                                                       gCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGCCGATACGCAGGt                                                                            >  1:263675/1‑147 (MQ=255)

                                                                       gCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGCCGATACGCAGGTCa                                                                          <  2:158285/148‑1 (MQ=255)

                                                                          gagaACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGCCGa                                                                                     >  3:191650/1‑134 (MQ=255)

                                                                          gagaACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGCCGa                                                                                     <  4:191650/134‑1 (MQ=255)

                                                                          gagaACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGCCGATACGCAGGTCAgt                                                                        >  6:3971/1‑147 (MQ=255)

                                                                          gagaACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCAGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGCCGATACGCAGGTCAgt                                                                        <  5:3971/147‑1 (MQ=255)

                                                                           agaACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGCCGATACGCAGGTCAGTGTa                                                                     >  7:301566/1‑149 (MQ=255)

                                                                              aCCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGCCGATACGCAGGTCAGTGTA‑‑Tg                                                                 <  1:375030/148‑1 (MQ=255)

                                                                                cTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑Ct                                                                                                                               >  8:357557/1‑87 (MQ=255)

                                                                                cTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑Ct                                                                                                                               <  7:357557/87‑1 (MQ=255)

                                                                                cTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑g                                                                                                     >  1:203006/1‑112 (MQ=255)

                                                                                cTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑g                                                                                                     <  2:203006/112‑1 (MQ=255)

                                                                                cTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGCCGATACGCAGGTCAGTGTA‑‑TGCAg                                                              >  3:303493/1‑149 (MQ=255)

                                                                                    aTGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGCCGATACGCAGGTCAGTGTA‑‑TGCAGGttt                                                          <  7:320672/149‑1 (MQ=255)

                                                                                         cATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGCCGATACGCAGGTCAGTGTA‑‑TGCAGGTTTGGTc                                                      <  6:189176/148‑1 (MQ=255)

                                                                                          aTGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGaa                                                                                                             >  2:254805/1‑95 (MQ=255)

                                                                                          aTGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGaa                                                                                                             <  1:254805/95‑1 (MQ=255)

                                                                                           tGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCCGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGCCGATACGCAGGTCAGTGTA‑‑TGCAGGTTTGGTCca                                                    >  6:332016/1‑148 (MQ=255)

                                                                                            ggCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGCCGATACGCAGGTCAGTGTA‑‑TGCAGGTTTGGt                                                       <  4:97801/144‑1 (MQ=255)

                                                                                            ggCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGCCGATACGCAGGTCAGTGTA‑‑TGCAGGTTTGGt                                                       >  3:97801/1‑144 (MQ=255)

                                                                                            ggCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGCCGATACGCAGGTCAGTGTA‑‑TGCAGGTTTGGTCcaca                                                  >  1:52700/1‑149 (MQ=255)

                                                                                             gCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGCCGATACGCAGGTCAGTGTA‑‑TGCAGGTTTGGTCcaca                                                  <  8:301566/148‑1 (MQ=255)

                                                                                              ccACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCCGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGCCGATACGCAGGTCAGTGTA‑‑TGCAGGTTTGGTCCACATg                                                <  5:332016/149‑1 (MQ=255)

                                                                                              ccACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGCCGATACGCAGGTCAGTGTA‑‑TGCAGGTTTGGTCCACATg                                                >  7:184228/1‑149 (MQ=255)

                                                                                                aCTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGCCGATACGCAGGTCAGTGTA‑‑TGCAGGTTTGGTCCACATgg                                               >  7:23742/1‑148 (MQ=255)

                                                                                                  tctctATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGCCGATACGCAGGTCAGTGTA‑‑TGCAGGTTTg                                                         <  1:5494/136‑1 (MQ=255)

                                                                                                  tctctATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGCCGATACGCAGGTCAGTGTA‑‑TGCAGGTTTg                                                         >  2:5494/1‑136 (MQ=255)

                                                                                                    tctATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGCCGATACGCAGGTCAGTGTA‑‑TGCAGGTTTGGTCCACATGGAAt                                            <  8:184228/147‑1 (MQ=255)

                                                                                                    tctATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGCCGATACGCAGGTCAGTGTA‑‑TGCAGGTTTGGTCCACATGGAAt                                            <  2:243099/147‑1 (MQ=255)

                                                                                                       aTGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGCCGATACGCAGGTCAGTGTA‑‑TGCAGGTTTGGTCCACATGGAATTACTc                                       >  3:253582/1‑149 (MQ=255)

                                                                                                        tGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGCCGATACGCAGGTCAGTGTA‑‑TGCAGGTTTGGTCCACATGGAATTACTcc                                      <  2:334649/149‑1 (MQ=255)

                                                                                                            gAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGCCGATACGCAGGTCAGTGTA‑‑TGCAGGTTTGGTCCACATGGAATTACTCCTTcc                                  <  3:44525/149‑1 (MQ=255)

                                                                                                                   ggATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGCCGATACGCAGGTCAGTGTA‑‑TGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGa                            <  1:172941/149‑1 (MQ=255)

                                                                                                                   ggATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGCCGATACGCAGGTCAGTGTA‑‑TGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGa                            <  5:190314/149‑1 (MQ=255)

                                                                                                                         aaaaa‑‑aCG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGCCGATACGCAGGTCAGTGTA‑‑TGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAgggttg                      >  7:333543/1‑149 (MQ=255)

                                                                                                                                acg‑‑‑‑‑‑‑‑‑‑‑aaaacgCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGCCGATACGCAGGTCAGTGTA‑‑TGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGggttgggttg                 <  1:43624/149‑1 (MQ=255)

                                                                                                                                              aaaacgCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGCCGATACGCAGGTCAGTGTA‑‑TGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTc              <  2:263675/149‑1 (MQ=255)

                                                                                                                                                 acgCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGCCGATACGCAGGTCAGTGTA‑‑TGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGc           >  8:42851/1‑149 (MQ=255)

                                                                                                                                                  cggcccGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGCCGATACGCAGGGCAGTGTA‑‑TGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGct          <  1:174826/146‑1 (MQ=255)

                                                                                                                                                  cgCCTCGGCCTGCCGGGGCGTGAAACCCTAAGCGCTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGCCTTCGAAGCCGATACGCAGGTCAGTGTA‑‑GGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGCGGGTTGGGTGGTTCAGc           >  1:51908/1‑149 (MQ=255)

                                                                                                                                                  cgCCTCGGCCTGCCGGGGCGTGAAACCCTAAGCGCTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGCCGATACGCAGGTCAGTGTA‑‑TGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGTTTGTTCAGc           >  1:46121/1‑149 (MQ=255)

                                                                                                                                                  cgCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGCCGATACGCAGGTCAGTGTA‑‑TGCAGGTTTgg                                                        >  4:366235/1‑103 (MQ=255)

                                                                                                                                                  cgCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGCCGATACGCAGGTCAGTGTA‑‑TGCAGGTTTgg                                                        <  3:366235/103‑1 (MQ=255)

                                                                                                                                                   gCCTCGGCCTGCCGGGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGCCGATACGCAGGTCAGTGTA‑‑TGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAgggttggg                    >  3:347919/1‑138 (MQ=255)

                                                                                                                                                   gCCTCGGCCTGCCGGGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGCCGATACGCAGGTCAGTGTA‑‑TGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAgggttggg                    <  4:347919/138‑1 (MQ=255)

                                                                                                                                                   gCCTCGGCCTGCCGGGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGCCGATACGCAGGTCAGTGTA‑‑TGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGcta         >  2:133158/1‑149 (MQ=255)

                                                                                                                                                   gCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGCCGATACGCAGGTCAGTGTA‑‑TGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGcta         >  1:271222/1‑149 (MQ=255)

                                                                                                                                                         gCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGCCGATACGCAGGTCAGTGTA‑‑TGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGCTACTACtt   <  2:52700/149‑1 (MQ=255)

                                                                                                                                                         gCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGCCGATACGCAGGTCAGTGTA‑‑TGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGCTACTACtt   >  1:357100/1‑149 (MQ=255)

                                                                                                                                                           aGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACATGGGTGACTTCGAAGCCGATACGCAGGTCAGTGTA‑‑TGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATCAGGGTTGGGTTGTTCAGATAGTAGTTg  >  1:94057/1‑149 (MQ=255)

                                                                                                                                                           aGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGCCGATACGCAGGTCAGTGTACCTGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGCTACTACtt   <  1:303363/149‑1 (MQ=255)

GTGCCATCACGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACAT‑GGTGACTTCGAAGCCGATACGCAGGTCAGTGTA‑‑TGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGCTACTACTTG  >  NC_002947/459763‑460051

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