BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	459,899	C→T	7.4%	intergenic (‑45/+34)	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,899	0	C	T	7.4%	318.6 / 10.3	108	intergenic (‑45/+34)	pqqB/pqqA	coenzyme PQQ synthesis protein B/coenzyme PQQ synthesis protein A
Reads supporting (aligned to +/- strand): ref base C (47/53); new base T (5/3); total (52/56)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

Only 100 of 109 total aligned reads displayed.

GTGCCATCACGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACC‑GGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGCTA  >  NC_002947/459763‑460044

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

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

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

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

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

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

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

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

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

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

          ggTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACC‑GGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAAcc                                                                                                                              >  2:207900/1‑149 (MQ=255)

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

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

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

                ccTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACC‑GGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCg                                                                                                                                                  >  7:260728/1‑123 (MQ=255)

                ccTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACC‑GGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCg                                                                                                                                                  <  8:260728/123‑1 (MQ=255)

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

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

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

                               agttgcagttCCACTGCGGGAACCCGCCACC‑GGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTa                                                                                                        <  8:254686/149‑1 (MQ=255)

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

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

                                   gcagttCCACTGCGGGAACCCGCCACCGGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGc                                                                                                      >  1:22978/1‑148 (MQ=255)

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

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

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

                                           aCTGCGGGAACCCGCCACC‑GGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACATg                                                                                            >  1:158285/1‑149 (MQ=255)

                                           aCTGCGGGAACCCGCCACC‑GGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACATg                                                                                            >  2:375030/1‑149 (MQ=255)

                                            cTGCGGGAACCCGCCACC‑GGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACATgg                                                                                           >  2:172941/1‑149 (MQ=255)

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

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

                                                 ggAACCCGCCACC‑GGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACATGGt                                                                                          <  7:265250/145‑1 (MQ=255)

                                                 ggAACCCGCCACC‑GGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACATGGt                                                                                          >  8:265250/1‑145 (MQ=255)

                                                 ggAACCCGCCACC‑GGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACATGGTGACt                                                                                      <  1:229996/149‑1 (MQ=255)

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

                                                     cccGCCACC‑GGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCg                                                                                   >  7:24624/1‑148 (MQ=255)

                                                      ccGCCACC‑GGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCg                                                                                   >  4:44525/1‑147 (MQ=255)

                                                      ccGCCACC‑GGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGaa                                                                                 <  1:115319/149‑1 (MQ=255)

                                                        gCCACC‑GGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACATGGTGCCTTCGAAGc                                                                               >  1:218419/1‑149 (MQ=255)

                                                        gCCACC‑GGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGc                                                                               <  1:207900/149‑1 (MQ=255)

                                                          cACC‑GGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGcc                                                                              >  2:329459/1‑148 (MQ=255)

                                                          cACC‑GGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGcc                                                                              >  5:189176/1‑148 (MQ=255)

                                                           aCC‑GGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGa                                                                            >  2:43624/1‑149 (MQ=255)

                                                           aCC‑GGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTATATGCAGAAACC‑GGATCTAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTCGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGa                                                                            <  2:187503/149‑1 (MQ=255)

                                                             c‑ggcggcggAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACATGGTGACtt                                                                                     >  7:4994/1‑138 (MQ=255)

                                                             c‑ggcggcggAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACATGGTGACtt                                                                                     <  8:4994/138‑1 (MQ=255)

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

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

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

                                                                   gcggAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAg                                                                     >  6:246397/1‑149 (MQ=255)

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

                                                                     ggAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGt                                                                   >  1:334649/1‑149 (MQ=255)

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

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

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

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

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

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

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

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

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

                                                                               aCCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATg                                                          <  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‑CTGGCTTAGCGGTTAGCGAAGTACATg                                                                                            <  2:203006/112‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                                                                         tGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTcc                               <  5:246397/149‑1 (MQ=255)

                                                                                                            aGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTc                            >  2:303363/1‑149 (MQ=255)

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

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

                                                                                                                    ggATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGa                     <  1:172941/149‑1 (MQ=255)

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

                                                                                                                                 acg‑‑‑‑‑‑‑‑‑‑‑aaaacgCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGggttgggttg          <  1:43624/149‑1 (MQ=255)

                                                                                                                                               aaaacgCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTc       <  2:263675/149‑1 (MQ=255)

                                                                                                                                                 cacgCCTCGGCCTGCCGGGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCa      <  3:220268/147‑1 (MQ=255)

                                                                                                                                                  acgCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGc    >  8:42851/1‑149 (MQ=255)

                                                                                                                                                   cggcccGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGGCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGct   <  1:174826/146‑1 (MQ=255)

                                                                                                                                                   cgCCTCGGCCTGCCGGGGCGTGAAACCCTAAGCGCTGGCTTAGCGGTTAGCGAAGTACATGGTGCCTTCGAAGCCGATACGCAGGTCAGTGTAGGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGCGGGTTGGGTGGTTCAGc    >  1:51908/1‑149 (MQ=255)

                                                                                                                                                   cgCCTCGGCCTGCCGGGGCGTGAAACCCTAAGCGCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGTTTGTTCAGc    >  1:46121/1‑149 (MQ=255)

                                                                                                                                                   cgCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTgg                                                 >  4:366235/1‑103 (MQ=255)

                                                                                                                                                   cgCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTgg                                                 <  3:366235/103‑1 (MQ=255)

                                                                                                                                                    gCCTCGGCCTGCCGGGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAgggttggg             >  3:347919/1‑138 (MQ=255)

                                                                                                                                                    gCCTCGGCCTGCCGGGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAgggttggg             <  4:347919/138‑1 (MQ=255)

                                                                                                                                                    gCCTCGGCCTGCCGGGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGcta  >  2:133158/1‑149 (MQ=255)

                                                                                                                                                    gCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGcta  >  1:271222/1‑149 (MQ=255)

GTGCCATCACGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACC‑GGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACC‑GGATCGAAAAA‑‑ACG‑‑‑‑‑‑‑‑‑‑‑AAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGC‑CTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGCTA  >  NC_002947/459763‑460044

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