BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,397,888	Δ1 bp	11.9%	intergenic (+103/‑83)	PP_1214 → / → ruvC	transcriptional regulator/crossover junction endodeoxyribonuclease RuvC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,397,888	0	T	.	11.9%	398.5 / 52.8	135	intergenic (+103/‑83)	PP_1214/ruvC	transcriptional regulator/crossover junction endodeoxyribonuclease RuvC
Reads supporting (aligned to +/- strand): ref base T (80/33); new base . (11/5); total (93/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.46e-01

Only 100 of 148 total aligned reads displayed.

ATTCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAA‑A‑TGTGGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGCCAGAC  >  NC_002947/1397739‑1398035

aTTCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAA‑A‑tgt                                                                                                                                                       <  2:2477054/150‑1 (MQ=255)

aTTCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCGTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAA‑A‑tgt                                                                                                                                                       <  2:5038758/150‑1 (MQ=255)

 ttCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAA‑A‑tgtg                                                                                                                                                      <  1:1582474/150‑1 (MQ=255)

 ttCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAA‑A‑tgtg                                                                                                                                                      <  1:4791163/150‑1 (MQ=255)

    cAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAA‑A‑Tgtggt‑g                                                                                                                                                  <  2:4594265/150‑1 (MQ=255)

      aTGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAA‑A‑TGTGGT‑GCa                                                                                                                                                <  2:3673309/150‑1 (MQ=255)

       tGCCCATATTTCCGACGAGATCATGGAAAACCTCGGCTTAGGTGATTTGGCATTTTGTGCCAACTCTGAGAGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGTTTAAACCCGCTACTACAGGTATTTCGCC‑CATTTCAA‑A‑TGTGGT‑GCAg                                                                                                                                               <  1:765026/150‑1 (MQ=255)

         cccAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAA‑A‑TGTGGT‑GCAGTa                                                                                                                                             <  1:3673228/150‑1 (MQ=255)

         cccAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGGAAACCCGCGCCGACAGGGATGACCCC‑CCATTCAA‑A‑GGTGGT‑Gccgga                                                                                                                                             >  2:2502923/1‑146 (MQ=255)

          ccAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAA‑A‑TGTGGT‑GCAGTAc                                                                                                                                            <  1:5019910/150‑1 (MQ=255)

          ccAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGATTACGCC‑CCTTTCAC‑A‑TGTGGT‑GCAGTAc                                                                                                                                            >  2:5470955/1‑150 (MQ=255)

           cAGATATCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAATACCGCTCCGAACGGTAATTAGCC‑CCATTCAA‑A‑TGGGGT‑GCAGTAcc                                                                                                                                           >  2:4303954/1‑150 (MQ=255)

            aGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAA‑A‑TGTGGT‑GCAGTACCt                                                                                                                                          <  2:5612935/150‑1 (MQ=255)

             gagttCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAA‑A‑TGTGGT‑GCAGTACCTg                                                                                                                                         <  1:3495757/147‑1 (MQ=255)

               tttCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTACTGCACCACATTTGAA‑A‑TG‑GGC‑GTAATACCTGTa                                                                                                                                       <  2:1897019/150‑1 (MQ=21)

               tttCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTACTGCACCACATTTGAA‑A‑TG‑GGC‑GTAATACCTGTa                                                                                                                                       <  1:3880848/150‑1 (MQ=21)

                ttCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAA‑A‑TGTGGT‑GCAGTACCTGTAg                                                                                                                                      <  1:4592392/150‑1 (MQ=255)

                    gacgaGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAA‑A‑TGTGGT‑GCAGTACCTGTAGGAGc                                                                                                                                  <  1:5462700/150‑1 (MQ=255)

                      cgaGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAA‑A‑TGTGGT‑GCAGTACCTGTAGGAGC‑gg                                                                                                                               <  2:3394817/150‑1 (MQ=255)

                      cgaGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAA‑A‑TGTGGT‑GCAGTACCTGTAGGAGC‑gg                                                                                                                               <  1:743076/150‑1 (MQ=255)

                      cgaGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAA‑A‑TGTGGT‑GCAGTACCTGTAGGAGC‑gg                                                                                                                               <  2:2782954/150‑1 (MQ=255)

                            cATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAA‑A‑TGTGGT‑GCAGTACCTGTAGGAGC‑GGGTTTAc                                                                                                                         <  2:1742333/150‑1 (MQ=255)

                            cATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAA‑A‑TGTGGT‑GCAGTACCTGTAGGAGC‑GGGTTTAc                                                                                                                         <  2:880129/150‑1 (MQ=255)

                             aTGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAA‑A‑TGTGGT‑GCAGTACCTGTAGGAGC‑GGGTTTAcc                                                                                                                        <  1:1696609/150‑1 (MQ=255)

                             aTGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAA‑A‑TGTGGT‑GCAGTACCTGTAGGAGC‑GGGTTTAcc                                                                                                                        <  1:4719153/150‑1 (MQ=255)

                             aTGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAA‑A‑TGTGGT‑GCAGTACCTGTAGGAGC‑GGGTTTAcc                                                                                                                        <  1:2966082/150‑1 (MQ=255)

                               ggAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCAACAGGTATTACGCC‑CATTTCAA‑A‑TGTGGT‑GCAGTACCGGGAGGAGC‑GGGTTTACCcg                                                                                                                      >  2:4215637/1‑150 (MQ=255)

                               ggAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAAACCGCTCCTACAGGGATTACGCC‑CAGTTAAA‑A‑TGTGGT‑GCAGTACCGGGAGGAGC‑GGGTTTAACcg                                                                                                                      >  2:4215517/1‑150 (MQ=38)

                                gAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAAACCGCTCCTACAGGTATTACGCC‑CAGTTAAA‑A‑AGTGGT‑GCAGTACCAGTAGGAGC‑GGGGTTAACcgc                                                                                                                     >  2:1273044/1‑150 (MQ=38)

                                         ggCTATGGTTATTTGGAATTTTGTGCCAACTCTGAGGGTTTTATAGCGCCTGTGCCGGCCCTTTCGCGGGTACACCCGCTCCTACCGGTATTACGCC‑CATTTAAA‑A‑GGTGGT‑GCAGTCCATGTAGGAGC‑GGGTTTACCCGCGAAAGGGcc                                                                                                            <  2:2945823/150‑1 (MQ=25)

                                            tAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAA‑A‑tgt                                                                                                                                                       >  2:3902378/1‑106 (MQ=255)

                                                                   cAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTTAACCCGCGCCAACACGGATTACCCC‑CAATTACA‑A‑AGTGGG‑GCAGTAACCGTAGGTGC‑CGGGTTAACCGCCCAAGGGCCCGCGCGGTTTTTAGGGATGcgaaaat                                                                                  >  2:287406/1‑144 (MQ=255)

                                                                   cAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCAACACGTAATACCCC‑CATTTCAA‑A‑TGTGGT‑GCAGTACCGGTAGGAGC‑GGGGTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGAAGCATCAtt                                                                                  >  2:1264932/1‑150 (MQ=37)

                                                                        ctGAGCGTTTGATAGCGCCGGGGCCGGCCCTGTCGGGGTTAAACCCGCTCCTACCGGT‑TTACCCC‑CCTTTCAA‑ATTGGGGT‑GCCGTACGTTTTGGAGC‑GGTTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTc                                                                             <  2:4973269/150‑1 (MQ=255)

                                                                          gAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CGTTTCAA‑A‑TGTGGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGGATCATTTTTTCGa                                                                           >  2:2242375/1‑150 (MQ=255)

                                                                           aGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCCCACCACAAGGTATTACGCC‑CCTTTCAA‑A‑TGTGGG‑GCAGTAACGGGAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCCGTTTGTAGAGATGCCTAATTTTTTCGAc                                                                          >  2:2948915/1‑150 (MQ=37)

                                                                             cGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCCCTCCTACAGGTATTACGCC‑CCAGTCAA‑A‑TGTGGT‑GCACTAACGGTAGGAGC‑CGGGTTAACCGCGAAAGGGCCCGCGCACGTTGTAGAGATGCAACAAGTTTTCGACgg                                                                        >  2:3730950/1‑150 (MQ=25)

                                                                               ttttATCGCGCGTGTGCCGGCCCTTGCGGGTTTAACCCGGCCCCCACCGGTATACCGCC‑CTTTCCAA‑T‑TGTGGT‑GCTTTACCTTTAGGCGG‑GGTTTAACCCGCGAATGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGc                                                                      <  2:4354331/150‑1 (MQ=255)

                                                                                 tGATAGCGCCTGTGCCGGCCCTTTCGCGGGTTAAACCCCTACGACACGTATTAAGCC‑CAATTAAA‑A‑AGTGGT‑GCAGTACAGGTAGGTGCGGGGTTAACCCCCGAACGGGCCGG‑GCCGGGTGTAGAGATTCGAGATTTTTTTGGCgggggt                                                                    >  2:3286867/1‑147 (MQ=255)

                                                                                  gATCGCGCCTGTGCCGGCCCTTTCGCGGTTAAACCCGCCCCTACTGTTATTACGCC‑CTTTTCAA‑A‑TGTGGT‑GCTGTACCTGTTGGAGG‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCttt                                                                   <  2:1846021/150‑1 (MQ=37)

                                                                                    tcgcgcGTGGGCCGGCCCTTGGGCGGGTACCCCGGCTCCTCCTGTTTTCCCGCC‑AATTGCAA‑T‑GGTGGT‑GCCGTCCGTGTTGGGGG‑GGTTTTACCCGGTAACGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTAt                                                                 <  2:2567521/148‑1 (MQ=255)

                                                                                    tcgcgcCGGGGCCGGCCGTTGCGGGGGTAACCCGGCGCGTCCCGGTATTCCGCC‑CTTTTCAA‑T‑GGTGGT‑GCTGTCCGGGTCGGGGG‑GGTTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTAt                                                                 <  2:5615009/148‑1 (MQ=255)

                                                                                        gcgtgggccgccCCTGTCGGGTCTAACCCGGCTCCACCCGGTATACCCCC‑CATTTC‑A‑A‑TGGGGTGGCCTCACCTTTAGGGGG‑GGTTTCACCCGCGAATGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGAct                                                             <  2:1606156/144‑1 (MQ=255)

                                                                                               ccGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAA‑A‑tgtg                                                                                                                                                      >  2:2775230/1‑56 (MQ=255)

                                                                                               ccGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAA‑A‑tgtg                                                                                                                                                      <  1:2775230/56‑1 (MQ=255)

                                                                                                 gccccTTGCGCGGTTAACCCCGCGCCTACCGTTATACCCCA‑CATTTCAA‑T‑GGTGGT‑GCCGACCCGGTTGGAGC‑GGTTTACCCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTa                                                    <  2:5281642/148‑1 (MQ=255)

                                                                                                  gCCCTTGCGGGGGTAAACCCGGCTCTAACGGTTTTACCCC‑CATTTAAATT‑TGGGGG‑GCCGGCCGTGTGGGAGG‑GTTTTCACCCGCGAACGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTa                                                    <  2:1860829/150‑1 (MQ=255)

                                                                                                    ccTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CCTTTCAA‑A‑TGTGGT‑GCAGTACCTGTGGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCCGTTTTTAGATATGCCTCATTTTTTCGACGGGGTTTATGACTCGGATTATAGGt                                                 >  2:756205/1‑150 (MQ=25)

                                                                                                             ggTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAA‑A‑TGTGGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTGGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGCTTCTAGGTATCGAcccc                                        >  2:4773723/1‑150 (MQ=255)

                                                                                                              gTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAA‑A‑TGTGGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCg                                       >  1:1693442/1‑150 (MQ=255)

                                                                                                               tAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAA‑A‑TGTGGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCg                                                                            >  1:2507321/1‑112 (MQ=255)

                                                                                                               tAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAA‑A‑TGTGGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCgg                                      >  2:175295/1‑150 (MQ=255)

                                                                                                                  aCCCGCTCCTACAGGTATTACGCC‑CATTTCAA‑A‑TGTGGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTc                                   >  1:3183683/1‑150 (MQ=255)

                                                                                                                  aCCCGCTCCTACAGGTATTACGCC‑CATTTCAA‑A‑TGTGGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTc                                   >  2:4487617/1‑150 (MQ=255)

                                                                                                                    tcGCTCCTACAGGTACTGCACCACATTTGAA‑A‑TG‑GGC‑GTAATACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTcgc                                 <  1:4154527/149‑1 (MQ=18)

                                                                                                                    ccGCTCCTACAGGTATTACGCC‑CATTTCAA‑A‑TGTGGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTcgc                                 >  1:3409097/1‑150 (MQ=255)

                                                                                                                    ccGCTCCTACAGGTACTGCACCACATTTGAA‑A‑TG‑GGC‑GTAATACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTcgc                                 >  2:4755345/1‑150 (MQ=18)

                                                                                                                     cGCTCCTACAGGTATTACGCC‑CATTTCAA‑A‑TGTGGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTcgcg                                >  1:2636855/1‑150 (MQ=255)

                                                                                                                      gCTCCTACAGGTACTGCACCACATTTGAA‑A‑TG‑GGC‑GTAATACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCGAGGTATCGACCCCGGCTCgcgc                               >  2:2909495/1‑150 (MQ=18)

                                                                                                                          cTACAGTTATTCCGCC‑AATTTAAA‑T‑TGTGGT‑GCAGCCCCTTTTGGAGG‑GGTTTAACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCa                           <  2:2165777/150‑1 (MQ=255)

                                                                                                                          cTACAGGTATTACGCC‑CATTTCAA‑A‑TGTGGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCa                           >  2:879690/1‑150 (MQ=255)

                                                                                                                          cTACAGGTATTACGCC‑CATTTCAA‑A‑TGTGGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCa                           >  2:4147370/1‑150 (MQ=255)

                                                                                                                           tACAGGTATTACGCC‑CATTTCAA‑A‑TGTGGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGAt                                                        >  1:2510497/1‑120 (MQ=255)

                                                                                                                           tACAGGTATTACGCC‑CATTTCAA‑A‑TGTGGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGAt                                                        <  2:2510497/120‑1 (MQ=255)

                                                                                                                           tACAGGTATTACGCC‑CATTTCAA‑A‑TGTGGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCAc                          >  1:2147393/1‑150 (MQ=255)

                                                                                                                           tACAGGTACTGCACCACATTTGAA‑A‑TG‑GGC‑GTAATACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCAc                          >  1:2292607/1‑150 (MQ=21)

                                                                                                                           tACAGGTACTGCACCACATTTGAA‑A‑TG‑GGC‑GTAATACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCAc                          >  1:5450429/1‑150 (MQ=21)

                                                                                                                           tACAGGTACTGCACCACATTTGAA‑A‑TG‑GGC‑GTAATACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCAc                          >  1:3965967/1‑150 (MQ=21)

                                                                                                                           tACAGGTACTGCACCACATTTGAA‑A‑TG‑GGA‑GTAATACCTGTAGGAGC‑GGGTTTACCCGCGAATGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCAc                          >  1:3452668/1‑150 (MQ=21)

                                                                                                                            aCAGGTATTACGCC‑CATTTCAA‑A‑TGTGGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCAcc                         >  1:1049564/1‑150 (MQ=255)

                                                                                                                            aCAGGTATTACGCC‑CATTTCAA‑A‑TGTGGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCAcc                         >  2:2199431/1‑150 (MQ=255)

                                                                                                                             cAGGTATTACGCC‑CATTTCAA‑A‑TGTGGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATAACCg                        >  1:4950502/1‑150 (MQ=255)

                                                                                                                              aGGTATTACGCC‑CATTTCAA‑A‑TGTGGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCgg                       >  1:4515657/1‑150 (MQ=255)

                                                                                                                              aGGTATTACGCC‑CATTTCAA‑A‑TGTGGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCgg                       >  1:4515405/1‑150 (MQ=255)

                                                                                                                              aGGTATTACGCC‑CATTTCAA‑A‑TGTGGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCgg                       >  2:613422/1‑150 (MQ=255)

                                                                                                                              aGGTATTACGCC‑CATTTCAA‑A‑TGTGGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCgg                       >  2:2589632/1‑150 (MQ=255)

                                                                                                                                gTATTACGCC‑CATTTCAA‑A‑TGTGGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCt                     >  2:2399730/1‑150 (MQ=255)

                                                                                                                                   ttACGCC‑CATTTCAA‑A‑TGTGGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACGGGCTATg                  >  1:46050/1‑150 (MQ=255)

                                                                                                                                   ttACGCC‑CATTTCAA‑A‑TGTGGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATg                  >  1:3425597/1‑150 (MQ=255)

                                                                                                                                      cGCC‑CATTTCAA‑T‑TGTGGT‑GCAGTCCCTGTCGGTGC‑GGGTTCACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCg               <  2:4497820/150‑1 (MQ=255)

                                                                                                                                       gCC‑CATTTCAA‑A‑TGTGGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTCGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGt              >  1:4073859/1‑150 (MQ=255)

                                                                                                                                       gCC‑CATTTCAA‑A‑TGTGGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATTACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGt              >  2:4056573/1‑150 (MQ=255)

                                                                                                                                       gCC‑CATTTCAA‑A‑TGTGGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGt              >  1:1434410/1‑150 (MQ=255)

                                                                                                                                       gCC‑CATTTCAA‑A‑TGTGGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGt              >  2:4206448/1‑150 (MQ=255)

                                                                                                                                       gCC‑CATTTCAA‑A‑TGTGGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGt              >  1:2368319/1‑150 (MQ=255)

                                                                                                                                       gCC‑CATTTCAA‑A‑TGTGGT‑GCAGTACCGGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGt              >  2:2578688/1‑150 (MQ=255)

                                                                                                                                       accacATTTGAA‑A‑TG‑GGC‑GTAATACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGt              >  1:576119/5‑150 (MQ=37)

                                                                                                                                        cc‑cATTTCAA‑A‑TGTGGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTg             <  1:2025211/150‑1 (MQ=255)

                                                                                                                                        cc‑cATTTCAA‑A‑TGTGGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTg             >  1:4702389/1‑150 (MQ=255)

                                                                                                                                        cc‑cATATCAA‑A‑AGGCGA‑GCAGTACCTGGAGGAGC‑CGGTTTACCCGCGAAAGGGCCGGCGCAGAGTGTAGAGATGCATAATATTTTAGACGCGCATTCTGACTCAGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTg             >  2:2561693/1‑150 (MQ=255)

                                                                                                                                         c‑cATTTCAA‑A‑TGTGGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTgg            >  2:4613762/1‑150 (MQ=255)

                                                                                                                                         c‑cATTTCAA‑A‑TGTGGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTgg            >  2:3669727/1‑150 (MQ=255)

                                                                                                                                         c‑cATTTCAA‑A‑TGTGGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTgg            >  1:663356/1‑150 (MQ=255)

                                                                                                                                             tttCAA‑A‑TGTGGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTAc         >  1:1994907/1‑150 (MQ=255)

                                                                                                                                             tttCAA‑A‑TGTGGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTAc         >  2:2564126/1‑150 (MQ=255)

                                                                                                                                             tttCAA‑A‑TGTGGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTAc         >  1:1935811/1‑150 (MQ=255)

                                                                                                                                              ttCAA‑A‑TGTGGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACg        >  2:1839605/1‑150 (MQ=255)

                                                                                                                                                 aa‑a‑TGTGGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGCCa     >  1:2886836/1‑150 (MQ=255)

                                                                                                                                                 aa‑a‑TGTGGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGCCa     >  1:2886470/1‑150 (MQ=255)

                                                                                                                                                  a‑a‑TGTGGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGCCAg    >  1:4768655/1‑150 (MQ=255)

ATTCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAA‑A‑TGTGGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGCCAGAC  >  NC_002947/1397739‑1398035

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