BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	1,345,484	A→C	9.4%	intergenic (‑188/+180)	fabI ← / ← ycjD	enoyl‑[acyl‑carrier‑protein] reductase, NADH‑dependent/DUF559 family endonuclease‑related protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	1,345,484	0	A	C	9.4%	474.9 / 11.2	179	intergenic (‑188/+180)	fabI/ycjD	enoyl‑[acyl‑carrier‑protein] reductase, NADH‑dependent/DUF559 family endonuclease‑related protein
Reads supporting (aligned to +/- strand): ref base A (45/112); new base C (20/2); total (65/114)
Fisher's exact test for biased strand distribution p-value = 2.24e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 181 total aligned reads displayed.

AATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCAACAAACAGGCTCAT‑ACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACCACAACAG  >  CP009273/1345335‑1345633

aaTCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCa                                                                                                                                                          <  1:2508164/150‑1 (MQ=255)

       cGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGc                                                                                                                                                   <  2:2277164/150‑1 (MQ=255)

         aTTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTGCCCCCCCCCAAACCCACAACCCCACAGGGGCg                                                                                                                                                  >  2:1264617/1‑150 (MQ=255)

             gtcgtTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCGTTTGCCCCCACCCCAACCCAC‑AACCCCCCGGGGCGC‑gggg                                                                                                                                            >  1:2193739/1‑150 (MQ=255)

             gtcgtTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCAGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTACCCCCCCCCCAAACCAC‑AACCCAAAAGGGCCC‑gggg                                                                                                                                            >  2:891811/1‑150 (MQ=255)

                gtTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGAcc                                                                                                                                         <  1:3273117/150‑1 (MQ=255)

                   ttAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTa                                                                                                                                      <  1:4075083/150‑1 (MQ=255)

                    tAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTAt                                                                                                                                     <  2:222617/150‑1 (MQ=255)

                    tAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCCTTTCCCCCCCCCCCAACCCAC‑AACCCACAAGGGCGC‑GGGGCCCGtct                                                                                                                                     >  1:4446999/1‑148 (MQ=255)

                     aGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTAtt                                                                                                                                    <  1:1050117/150‑1 (MQ=255)

                     aGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTAtt                                                                                                                                    <  1:2624174/150‑1 (MQ=255)

                      gTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTTCCCGCCCCCCAAACCAC‑ACCCCACAGGGGGGG‑CGGGGCCGGAttt                                                                                                                                   >  2:3126720/1‑149 (MQ=255)

                      gTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATtg                                                                                                                                   <  1:446006/150‑1 (MQ=255)

                        aaaCAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATtgt‑g                                                                                                                                <  1:2951459/150‑1 (MQ=255)

                         aaCAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCGTTTCCCCGCACCCAAACCCTC‑GACCCAGCGGGGCGA‑GGGGCCCCTAGTGG‑Gc                                                                                                                               >  2:2923234/1‑150 (MQ=255)

                          aCAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCa                                                                                                                              <  2:394052/150‑1 (MQ=255)

                            aGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCACACCCTAACCCAC‑TCCCCCGAGGGGGGT‑GGGGTCCGGATTGT‑GCaaa                                                                                                                            >  2:643300/1‑150 (MQ=255)

                             gTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTGCCCGCCCCCGCAACCAC‑ACCCCCCAGGGGGGC‑TGGGCCCCTATTTT‑TCaaaa                                                                                                                           >  2:1569396/1‑149 (MQ=255)

                              tACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAata                                                                                                                          <  2:1801097/150‑1 (MQ=255)

                              tACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAata                                                                                                                          <  2:1117212/150‑1 (MQ=255)

                               aCGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAatat                                                                                                                         <  2:1514702/150‑1 (MQ=255)

                               aCGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCGTTTGCCCACACCCCAACCCAC‑ACCCCCCACGGGCGC‑GGGGCCCCGATTTT‑GCAAataa                                                                                                                         >  2:1987474/1‑149 (MQ=255)

                                 gAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACGCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCACATATTg                                                                                                                       >  1:4332631/1‑150 (MQ=255)

                                  aaCAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTTCCCACACCCAAACCCAC‑ACCCCCCAGGGGCGC‑GGGGGCCCTAGTGT‑GCCAataattt                                                                                                                      >  2:3587856/1‑146 (MQ=255)

                                   aCAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGtt                                                                                                                     <  1:942778/150‑1 (MQ=255)

                                    cAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTa                                                                                                                    <  1:3777883/150‑1 (MQ=255)

                                     aGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTAc                                                                                                                   <  1:3457712/150‑1 (MQ=255)

                                     aGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTAc                                                                                                                   <  1:3457426/150‑1 (MQ=255)

                                      gATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTACCCACACCCAACACCAC‑ACCCCAGAGGGGCGA‑GGGGTCCCGATTGT‑GGCAATATTGTTaac                                                                                                                  >  2:4003210/1‑148 (MQ=255)

                                        tAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCGTTTGACCGCAACCGAAACCAC‑AAACCCCAGGGGCGG‑GGGGGCCGATTTGT‑GGAAAAATtttttaac                                                                                                                 >  1:1833592/1‑146 (MQ=255)

                                        tAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTAcccc                                                                                                                <  1:527381/150‑1 (MQ=255)

                                        tAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCCCCACCAAAACCCC‑ACCCCACAGGGGCGA‑GGGGACCCTAATGTGGAAAAAATTTTTAccc                                                                                                                 >  2:2421196/1‑150 (MQ=255)

                                         aaaCGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCa                                                                                                               <  2:3683255/150‑1 (MQ=255)

                                           aCGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCGTTTCCCCACCCCCAAAACCAC‑ACCCCCGAGGGGCGG‑GGGGGCCCGATTGT‑GCAAATATTTTTAccccccc                                                                                                             >  2:1798320/1‑147 (MQ=255)

                                             ggTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGcaa                                                                                                           <  1:3154778/150‑1 (MQ=255)

                                              gTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTACCCACACCCAAACCCAC‑TCCCCCCAGGGGCGG‑GGGGACCCGATTGT‑GCAAAAAATGTTAACCCCCcaac                                                                                                          >  1:4145384/1‑150 (MQ=255)

                                                  ttatAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTGCCCCCCCCCCAAACCAC‑ACCCCCGAGGGGGGAGGGGGGCCGGTTTGT‑GGAAAATTTTTTACCCC‑CCaacaaac                                                                                                      >  1:4527180/1‑150 (MQ=255)

                                                  ttatAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCaacaaac                                                                                                      <  1:1321008/150‑1 (MQ=255)

                                                  ttatAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGAGACCGTATTGT‑GCAAATATTGTTACCCCAGCaacaaac                                                                                                      >  2:4132778/1‑150 (MQ=255)

                                                  ttatAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCACACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGAGACCGTATTGT‑GCAAATATTGTTACCCCAGCaacaaac                                                                                                      >  2:4132780/1‑150 (MQ=255)

                                                   tataATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCaacaaaca                                                                                                     <  1:3438452/150‑1 (MQ=255)

                                                   tataATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCaacaaaca                                                                                                     <  2:2418600/150‑1 (MQ=255)

                                                    ataATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCAACAAACAg                                                                                                    <  1:3500996/150‑1 (MQ=255)

                                                       aTCCACCTGGCGGTGAGTAGCGATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCAACAAACAGGCt                                                                                                 <  2:3338240/150‑1 (MQ=255)

                                                       aTCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCAACAAACAGGCt                                                                                                 <  2:376884/150‑1 (MQ=255)

                                                       aTCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCAACAAACAGGCt                                                                                                 <  2:377563/150‑1 (MQ=255)

                                                         cAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCGTTTCCCCACCCCCCAACCCCC‑ACCCCCCCGGGGCGC‑GGGGCCCGGATTGT‑GGAAATATTGTTACCCCACCAACCAACACGCGCa                                                                                               >  2:1471345/1‑150 (MQ=255)

                                                               ggCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCAACAAACAGGCTCAT‑ACAGc                                                                                        <  1:1029576/150‑1 (MQ=255)

                                                                gCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCAACAAACAGGCTCAT‑ACAGcc                                                                                       <  1:2471373/150‑1 (MQ=255)

                                                                  tgtgAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCACACCCAAACCCAC‑AACCCAGAGGGGCGA‑TGGGGCCCTATTTT‑GGAAAAAATGTTACCCCACCAAACAACACGGGCATAACAGccc                                                                                      >  2:300359/1‑150 (MQ=255)

                                                                   gtgAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTGCCCGCACCCAAACCCAC‑ACCCCCGACGGGGGG‑GGGGGCCGGATTGT‑GGAAAAATTGTTACCCCCCCAACAAAAACGCTCAA‑AAAGccccc                                                                                    >  2:208733/1‑149 (MQ=255)

                                                                   gtgAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCAACAAACAGGCTCAT‑ACAGCCCCt                                                                                    <  1:1806072/150‑1 (MQ=255)

                                                                   gtgAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCAACAAACAGGCTCAT‑ACAGCCCCt                                                                                    <  1:1420699/150‑1 (MQ=255)

                                                                    tgAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCAACAAACAGGCTCAT‑ACAGCCCCTa                                                                                   <  2:1856933/150‑1 (MQ=255)

                                                                      aGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCAACAAACAGGCTCAT‑ACAGCCCCTAAc                                                                                 <  1:17202/150‑1 (MQ=255)

                                                                      aGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTACCCGCACCCAAAACCAC‑TCCCCAGAGGGGCGA‑GGGGACCGGATTGT‑GCAAAAAATGTTAACCCCGCAACAAACAGGCTCAT‑ACAGCCCCAAAc                                                                                 >  2:754418/1‑150 (MQ=255)

                                                                           cTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCAACAAACAGGCTCAT‑ACAGCCCCTAACCCttt                                                                            <  1:1115509/150‑1 (MQ=255)

                                                                             ataGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTTCCCCCCCCCAAAACCAA‑AACCCCCAGGGGGGC‑GGGGGCCCGCTTGT‑GGCAATATTGTTTACCCCCCCACAAAAACGCGCAA‑AAACCCCCAAAcccagtca                                                                          >  2:1350636/1‑145 (MQ=255)

                                                                              taGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCAACAAACAGGCTCAT‑ACAGCCCCTAACCCTTTCAt                                                                         <  1:3307215/150‑1 (MQ=255)

                                                                                 ttGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCAACAAACAGGCTCAT‑ACAGCCCCTAACCCTTTCATGGc                                                                      <  2:3931356/150‑1 (MQ=255)

                                                                                 gtGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCAACAAACAGGCTCAT‑ACAGCCCCTAACCCTTTCATGGc                                                                      <  1:2118070/149‑1 (MQ=255)

                                                                                          ccGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCAACAAACAGGCTCAT‑ACAGCCCCTAACCCTTTCATGGCGATGGCTg                                                              <  2:2557568/150‑1 (MQ=255)

                                                                                          ccGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCAACAAACAGGCTCAT‑ACAGCCCCTAACCCTTTCATGGCGATGGCTg                                                              <  2:2555110/150‑1 (MQ=255)

                                                                                           cGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCAACAAACAGGCTCAT‑ACAGCCCCTAACCCTTTCATGGCGATGGCTggg                                                            <  2:2419783/150‑1 (MQ=255)

                                                                                              ccGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCAACAAACAGGCTCAT‑ACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACg                                                         <  1:372740/150‑1 (MQ=255)

                                                                                                  aGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCAACAAACAGGCTCAT‑ACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTc                                                     <  2:582513/150‑1 (MQ=255)

                                                                                                    cAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCCGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCAACAAACAGGCTCAT‑ACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAg                                                   <  1:3713854/150‑1 (MQ=255)

                                                                                                      ggCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCAACAAACAGGCTCAT‑ACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGAc                                                 <  2:855189/150‑1 (MQ=255)

                                                                                                       gCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCAACAAACAGGCTCAT‑ACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGAcc                                                <  1:2189781/150‑1 (MQ=255)

                                                                                                       gCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCAACAAACAGGCTCAT‑ACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGAcc                                                <  1:3513064/150‑1 (MQ=255)

                                                                                                        cTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCAACAAACAGGCTCAT‑ACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCt                                               <  1:2773219/150‑1 (MQ=255)

                                                                                                         tGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCAACAAACAGGCTCAT‑ACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCtt                                              <  1:1539429/150‑1 (MQ=255)

                                                                                                          gCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCAACAAACAGGCTCAT‑ACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTg                                             <  2:637531/150‑1 (MQ=255)

                                                                                                            ggCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCAACAAACAGGCTCAT‑ACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGcc                                           <  1:2664884/150‑1 (MQ=255)

                                                                                                               a‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCAACAAACAGGCTCAT‑ACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGaa                                        <  1:1529854/150‑1 (MQ=255)

                                                                                                               a‑GGGGGGGCGCTTGTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCAACAAACAGGCTCAT‑ACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGaa                                        <  1:2735920/150‑1 (MQ=255)

                                                                                                                     gggCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCAACAAACAGGCTCAT‑ACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATtc                                   <  2:1699913/150‑1 (MQ=255)

                                                                                                                      ggCGCTTTTCCCCTCACCCTAATCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCAACAAACAGGCTCAT‑ACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATtct                                  <  1:855495/150‑1 (MQ=255)

                                                                                                                        cgcTTTTCCCCACACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCAACAAACAGGCTCAT‑ACAGCCCCTCACCCGTTCATGGCGCTGGCTGGGACGGTTCAGTCCTTGCCGAATATTCTcc                                >  2:1777918/1‑150 (MQ=39)

                                                                                                                         gcTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCAACAAACAGGCTCAT‑ACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCa                               <  1:75764/150‑1 (MQ=255)

                                                                                                                          cTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCAACAAACAGGCTCAT‑ACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAg                              >  1:1706878/1‑150 (MQ=255)

                                                                                                                              tCCCCTCACCCAAACCCAC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCCACAAACAGGCTCAT‑ACAGCCCCAAACCCGTTCATTGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCCGCAcc                          >  2:787198/1‑150 (MQ=39)

                                                                                                                                cccTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCAACAAACAGGCTCAT‑ACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGt                        <  1:1577025/150‑1 (MQ=255)

                                                                                                                                 ccTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCAACAAACAGGCTCAT‑ACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGtc                       <  1:4239374/150‑1 (MQ=255)

                                                                                                                                 ccTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCAACAAACAGGCTCAT‑ACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGtc                       <  1:4239525/150‑1 (MQ=255)

                                                                                                                                  cTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCAACAAACAGGCTCAT‑ACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGtct                      >  2:3390859/1‑150 (MQ=255)

                                                                                                                                   tCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCAACAAACAGGCTCAT‑ACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGtctc                     >  2:1909127/1‑150 (MQ=255)

                                                                                                                                    cACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCAACAAACAGGCTCAT‑ACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTcc                    >  2:442088/1‑150 (MQ=255)

                                                                                                                                    cACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCAACAAACAGGCTCAT‑ACAGCCCCTAACCATTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTcc                    >  1:1684946/1‑150 (MQ=255)

                                                                                                                                      cccTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCAACAAACAGGCTCAT‑ACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCAt                  <  1:3750988/150‑1 (MQ=255)

                                                                                                                                       ccTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCAACAAACAGGCTCAT‑ACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATg                 >  1:2318931/1‑150 (MQ=255)

                                                                                                                                         tAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCAACAAACAGGCTCAT‑ACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGtt               <  1:2145133/150‑1 (MQ=255)

                                                                                                                                         tAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCAACAAACAGGCTCAT‑ACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGtt               <  2:1073569/150‑1 (MQ=255)

                                                                                                                                              ctc‑tCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCAACAAACAGGCTCCT‑ACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTcac          >  2:1743479/1‑150 (MQ=255)

                                                                                                                                              ctc‑tCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCAACAAACAGGCTCAT‑ACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTcac          >  2:4303646/1‑150 (MQ=255)

                                                                                                                                               tc‑tcCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCAACAAACAGGCTCAT‑ACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTcacc         >  1:4354168/1‑150 (MQ=255)

                                                                                                                                               tc‑tcCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCAACAAACAGGCTCAT‑ACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTcacc         <  1:1987474/150‑1 (MQ=255)

                                                                                                                                                    cccAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCCGCAACAAACAGGCTCAT‑ACAGCCCCTAACCCTTTCCTGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACCacaa     <  2:2158764/150‑1 (MQ=255)

                                                                                                                                                    cccAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCAACAAACAGGCTCAT‑ACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACCacaa     >  1:2153188/1‑150 (MQ=255)

                                                                                                                                                     ccAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCAACAAACAGGCTCAT‑ACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACCacaac    >  1:2938662/1‑150 (MQ=255)

AATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTC‑TCCCCAGAGGGGCGA‑GGGGACCGTATTGT‑GCAAATATTGTTACCCCAGCAACAAACAGGCTCAT‑ACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACCACAACAG  >  CP009273/1345335‑1345633

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 11 ≤` `ATCG/ATCG` `< 25 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`