BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,397,895	G→A	12.1%	intergenic (+110/‑76)	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,895	0	G	A	12.1%	357.4 / 32.3	136	intergenic (+110/‑76)	PP_1214/ruvC	transcriptional regulator/crossover junction endodeoxyribonuclease RuvC
Reads supporting (aligned to +/- strand): ref base G (85/31); new base A (11/5); total (97/39)
Fisher's exact test for biased strand distribution p-value = 7.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.74e-01

Only 100 of 140 total aligned reads displayed.

TGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAAA‑TGT‑GGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGCCAGACCGCCCGTGG  >  NC_002947/1397746‑1398044

   ccAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTTAACCCGCTCCAACAGGTATTACGCC‑CCATTCAAA‑TGT‑GGT‑GCCGTAc                                                                                                                                                     >  2:3463578/1‑150 (MQ=255)

   ccAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAAA‑TGT‑GGT‑GCAGTAc                                                                                                                                                     <  1:5019910/150‑1 (MQ=255)

    cAGATATCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAATACCGCTCCGAACGGTAATTAGCC‑CCATTCAAA‑TGG‑GGT‑GCAGTAcc                                                                                                                                                    >  2:4303954/1‑150 (MQ=255)

     aGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAAA‑TGT‑GGT‑GCAGTACCt                                                                                                                                                   <  2:5612935/150‑1 (MQ=255)

      gagttCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAAA‑TGT‑GGT‑GCAGTACCTg                                                                                                                                                  <  1:3495757/147‑1 (MQ=255)

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

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

        tttCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTACTGCACCACATTTGAAA‑TG‑‑GGC‑GTAATACCTGTa                                                                                                                                                <  2:3255665/150‑1 (MQ=21)

         ttCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAAA‑TGT‑GGT‑GCAGTACCTGTAg                                                                                                                                               <  1:4592392/150‑1 (MQ=255)

         ttCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAAA‑TGT‑GGT‑GCAGTACCTGTAg                                                                                                                                               <  1:345264/150‑1 (MQ=255)

             gacgaGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAAA‑TGT‑GGT‑GCAGTACCTGTAGGAGc                                                                                                                                           <  1:5462700/150‑1 (MQ=255)

             gacgaGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTACTGCACCACATTTGAAA‑TG‑‑GGC‑GTAATACCTGTAGGAGc                                                                                                                                           <  2:4138137/150‑1 (MQ=18)

               cgaGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAAA‑TGT‑GGT‑GCAGTACCTGTAGGAGC‑gg                                                                                                                                        <  2:1841702/150‑1 (MQ=255)

               cgaGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAAA‑TGT‑GGT‑GCAGTACCTGTAGGAGC‑gg                                                                                                                                        <  2:2782954/150‑1 (MQ=255)

               cgaGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAAA‑TGT‑GGT‑GCAGTACCTGTAGGAGC‑gg                                                                                                                                        <  1:743076/150‑1 (MQ=255)

               cgaGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAAA‑TGT‑GGT‑GCAGTACCTGTAGGAGC‑gg                                                                                                                                        <  2:5089330/150‑1 (MQ=255)

                     cATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAAA‑TGT‑GGT‑GCAGTACCTGTAGGAGC‑GGGTTTAc                                                                                                                                  <  2:880129/150‑1 (MQ=255)

                      aTGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAAA‑TGT‑GGT‑GCAGTACCTGTAGGAGC‑GGGTTTAcc                                                                                                                                 <  1:4719280/150‑1 (MQ=255)

                      aTGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAAA‑TGT‑GGT‑GCAGTACCTGTAGGAGC‑GGGTTTAcc                                                                                                                                 <  1:2966082/150‑1 (MQ=255)

                        ggAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCAACAGGTATTACGCC‑CATTTCAAA‑TGT‑GGT‑GCAGTACCGGGAGGAGC‑GGGTTTACCcg                                                                                                                               >  2:4215637/1‑150 (MQ=255)

                        ggAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAAACCGCTCCTACAGGGATTACGCC‑CAGTTAAAA‑TGT‑GGT‑GCAGTACCGGGAGGAGC‑GGGTTTAACcg                                                                                                                               >  2:4215517/1‑150 (MQ=38)

                              ccTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAAA‑TGT‑GGT‑GCAGTAACTGTAGGAGC‑GGGTTTACCCGCGAAAg                                                                                                                         >  2:4210017/1‑150 (MQ=255)

                                  ggCTATGGTTATTTGGAATTTTGTGCCAACTCTGAGGGTTTTATAGCGCCTGTGCCGGCCCTTTCGCGGGTACACCCGCTCCTACCGGTATTACGCC‑CATTTAAAA‑GGT‑GGT‑GCAGTCCATGTAGGAGC‑GGGTTTACCCGCGAAAGGGcc                                                                                                                     <  2:2945823/150‑1 (MQ=25)

                                                            cAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTTAACCCGCGCCAACACGGATTACCCC‑CAATTACAA‑AGT‑GGG‑GCAGTAACCGTAGGTGC‑CGGGTTAACCGCCCAAGGGCCCGCGCGGTTTTTAGGGATGcgaaaat                                                                                           >  2:287406/1‑144 (MQ=255)

                                                            cAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCAACACGTAATACCCC‑CATTTCAAA‑TGT‑GGT‑GCAGTACCGGTAGGAGC‑GGGGTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGAAGCATCAtt                                                                                           >  2:1264932/1‑150 (MQ=37)

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

                                                                      cGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCCCTCCTACAGGTATTACGCC‑CCAGTCAAA‑TGT‑GGT‑GCACTAACGGTAGGAGC‑CGGGTTAACCGCGAAAGGGCCCGCGCACGTTGTAGAGATGCAACAAGTTTTCGACgg                                                                                 >  2:3730950/1‑150 (MQ=25)

                                                                        ttttATCGCGCGTGTGCCGGCCCTTGCGGGTTTAACCCGGCCCCCACCGGTATACCGCC‑CTTTCCAAT‑TGT‑GGT‑GCTTTACCTTTAGGCGG‑GGTTTAACCCGCGAATGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGc                                                                               <  2:4354331/150‑1 (MQ=255)

                                                                          tGATAGCGCCTGTGCCGGCCCTTTCGCGGGTTAAACCCCTACGACACGTATTAAGCC‑CAATTAAAA‑AGT‑GGT‑GCAGTACAGGTAGGTGCGGGGTTAACCCCCGAACGGGCCGG‑GCCGGGTGTAGAGATTCGAGATTTTTTTGGCgggggt                                                                             >  2:3286867/1‑147 (MQ=255)

                                                                           gATCGCGCCTGTGCCGGCCCTTTCGCGGTTAAACCCGCCCCTACTGTTATTACGCC‑CTTTTCAAA‑TGT‑GGT‑GCTGTACCTGTTGGAGG‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCttt                                                                            <  2:1846021/150‑1 (MQ=37)

                                                                             tcgcgcGTGGGCCGGCCCTTGGGCGGGTACCCCGGCTCCTCCTGTTTTCCCGCC‑AATTGCAAT‑GGT‑GGT‑GCCGTCCGTGTTGGGGG‑GGTTTTACCCGGTAACGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTAt                                                                          <  2:2567521/148‑1 (MQ=255)

                                                                                 gcgtgggccgccCCTGTCGGGTCTAACCCGGCTCCACCCGGTATACCCCC‑CATTTC‑AA‑TGG‑GGTGGCCTCACCTTTAGGGGG‑GGTTTCACCCGCGAATGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGAct                                                                      <  2:1606156/144‑1 (MQ=255)

                                                                                   cTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTAAAGGTATTACGCC‑CATTTCAAA‑TGT‑GGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCCGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGActct                                                                    >  2:5545014/1‑150 (MQ=255)

                                                                                          gccccTTGCGCGGTTAACCCCGCGCCTACCGTTATACCCCA‑CATTTCAAT‑GGT‑GGT‑GCCGACCCGGTTGGAGC‑GGTTTACCCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTa                                                             <  2:5281642/148‑1 (MQ=255)

                                                                                             ccTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CCTTTCAAA‑TGT‑GGT‑GCAGTACCTGTGGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCCGTTTTTAGATATGCCTCATTTTTTCGACGGGGTTTATGACTCGGATTATAGGt                                                          >  2:756205/1‑150 (MQ=25)

                                                                                                      ggTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAAA‑TGT‑GGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTGGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGCTTCTAGGTATCGAcccc                                                 >  2:4773723/1‑150 (MQ=255)

                                                                                                       gTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAAA‑TGT‑GGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCg                                                >  1:1693442/1‑150 (MQ=255)

                                                                                                        tAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAAA‑TGT‑GGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCg                                                                                     <  2:2507321/112‑1 (MQ=255)

                                                                                                        tAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAAA‑TGT‑GGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCgg                                               >  2:175295/1‑150 (MQ=255)

                                                                                                        tAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAAA‑TGT‑GCT‑GCAGAACCGGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGGTTGTAGAGATGCATCATTTTTACGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCgg                                               >  2:175264/1‑150 (MQ=38)

                                                                                                           aCCCGCTCCTACAGGTATTACGCC‑CATTTCAAA‑TGT‑GGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTc                                            >  1:3183683/1‑150 (MQ=255)

                                                                                                           aCCCGCTCCTACAGGTATTACGCC‑CATTTCAAA‑TGT‑GGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTc                                            >  2:4487617/1‑150 (MQ=255)

                                                                                                             ccGCTCCTACAGGTATTACGCC‑CATTTCAAA‑TGT‑GGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTcgc                                          >  1:3409097/1‑150 (MQ=255)

                                                                                                             ccGCTCCTACAGGTATTACGCC‑CATTTCAAA‑TGT‑GGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTcgc                                          >  1:3549983/1‑150 (MQ=255)

                                                                                                             ccGCTCCTACAGGTACTGCACCACATTTGAAA‑TG‑‑GGC‑GTAATACCTGTAGGAGC‑GGGTTTACCCGCGAACGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTcgc                                          >  1:1882002/1‑150 (MQ=18)

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

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

                                                                                                               gCTCCTACAGGTACTGCACCACATTTGAAA‑TG‑‑GGC‑GTAATACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCACAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCgcgc                                        >  1:4518027/1‑150 (MQ=18)

                                                                                                                 tCCTACAGGTACTGCACCACATTTGAAA‑TG‑‑GGC‑GTAATACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGAAGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCAt                                      >  2:2563398/1‑150 (MQ=18)

                                                                                                                  ccTACAGGTATTACGCC‑CATTTCAAA‑TGT‑GGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATc                                     >  1:258441/1‑150 (MQ=255)

                                                                                                                   cTACAGGTATTACGCC‑CATTTCAAA‑TGT‑GGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCa                                    >  2:879690/1‑150 (MQ=255)

                                                                                                                   cTACAGGTATTACGCC‑CATTTCAAA‑TGT‑GGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCa                                    >  2:4147370/1‑150 (MQ=255)

                                                                                                                   cTACAGGTATTACGCC‑CATTTCAAA‑TGT‑GGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCa                                    >  2:1294299/1‑150 (MQ=255)

                                                                                                                    tACAGGTATTACGCC‑CATTTCAAA‑TGT‑GGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGAt                                                                 >  1:2510497/1‑120 (MQ=255)

                                                                                                                    tACAGGTATTACGCC‑CATTTCAAA‑TGT‑GGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCAc                                   >  1:2147393/1‑150 (MQ=255)

                                                                                                                    tACAGGTATTACGCC‑CATTTCAAA‑TGT‑GGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCAc                                   >  2:4519346/1‑150 (MQ=255)

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

                                                                                                                    tACAGGTACTGCACCACATTTGAAA‑TG‑‑GGC‑GTAATACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCAc                                   >  2:2658925/1‑150 (MQ=21)

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

                                                                                                                     aCAGGTATTACGCC‑CATTTCAAA‑TGT‑GGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCAcc                                  >  1:1049564/1‑150 (MQ=255)

                                                                                                                     aCAGGTATTACGCC‑CATTTCAAA‑TGT‑GGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCAcc                                  >  2:2199431/1‑150 (MQ=255)

                                                                                                                     aCAGGTATTACGCC‑CATTTCAAA‑TGT‑GGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCAcc                                  >  2:2199441/1‑150 (MQ=255)

                                                                                                                       aGGTATTACGCC‑CATTTCAAA‑TGT‑GGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCgg                                >  2:613422/1‑150 (MQ=255)

                                                                                                                       aGGTATTACGCC‑CATTTCAAA‑TGT‑GGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCgg                                >  1:2450205/1‑150 (MQ=255)

                                                                                                                       aGGTATTACGCC‑CATTTCAAA‑TGT‑GGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCgg                                >  1:4515405/1‑150 (MQ=255)

                                                                                                                        ggTATTACGCC‑CATTTCAAA‑TGT‑GGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGc                               >  1:280952/1‑150 (MQ=255)

                                                                                                                         gTATTACGCC‑CATTTCAAA‑TGT‑GGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCt                              >  2:2399730/1‑150 (MQ=255)

                                                                                                                         gTATTACGCC‑CATTTCAAA‑TGT‑GGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCt                              >  1:4950151/1‑150 (MQ=255)

                                                                                                                            ttACGCC‑CATTTCAAA‑TGT‑GGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATg                           >  1:727019/1‑150 (MQ=255)

                                                                                                                            ttACGCC‑CATTTCAAA‑TGT‑GGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATg                           >  1:3425597/1‑150 (MQ=255)

                                                                                                                            attcccc‑cATTTCCAA‑CGGGGGT‑GCATCACCTTTTGGCGC‑GTTTTAACCCGCGAATGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTAt                            <  1:460543/145‑1 (MQ=255)

                                                                                                                               cGCC‑CATTTCAAT‑TGT‑GGT‑GCAGTCCCTGTCGGTGC‑GGGTTCACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCg                        <  2:4497820/150‑1 (MQ=255)

                                                                                                                                gCC‑CATTTCAAA‑TGT‑GGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTCGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGt                       >  1:4073859/1‑150 (MQ=255)

                                                                                                                                gCC‑CATTTCAAA‑TGT‑GGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGt                       >  1:2368319/1‑150 (MQ=255)

                                                                                                                                gCC‑CATTTCAAA‑TGT‑GGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGt                       >  1:1434410/1‑150 (MQ=255)

                                                                                                                                gCC‑CATTTCAAA‑TGT‑GGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGt                       >  2:2578665/1‑150 (MQ=255)

                                                                                                                                gCC‑CATTTCAAA‑TGT‑GGT‑GCAGTACCGGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGt                       >  2:2578688/1‑150 (MQ=255)

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

                                                                                                                                 cc‑cATTTCAAA‑TGT‑GGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTg                      <  1:2025211/150‑1 (MQ=255)

                                                                                                                                 cc‑cATTTCAAA‑TGT‑GGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTg                      >  1:4702389/1‑150 (MQ=255)

                                                                                                                                 cc‑cATATCAAA‑AGG‑CGA‑GCAGTACCTGGAGGAGC‑CGGTTTACCCGCGAAAGGGCCGGCGCAGAGTGTAGAGATGCATAATATTTTAGACGCGCATTCTGACTCAGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTg                      >  2:2561693/1‑150 (MQ=255)

                                                                                                                                  c‑cATTTCAAA‑TGT‑GGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTgg                     >  2:3669727/1‑150 (MQ=255)

                                                                                                                                  c‑cATTTCAAA‑TGT‑GGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTgg                     >  1:5506699/1‑150 (MQ=255)

                                                                                                                                  c‑cATTTCAAA‑TGT‑GGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTgg                     >  1:2024596/1‑150 (MQ=255)

                                                                                                                                      tttCAAA‑TGT‑GGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTAc                  >  2:2457927/1‑150 (MQ=255)

                                                                                                                                      tttCAAA‑TGT‑GGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTAc                  >  1:1994907/1‑150 (MQ=255)

                                                                                                                                      tttCAAA‑TGT‑GGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTAc                  >  2:5618778/1‑150 (MQ=255)

                                                                                                                                       ttCAAA‑TGT‑GGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACg                 >  2:1027137/1‑150 (MQ=255)

                                                                                                                                          aaa‑TGT‑GGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGCCa              <  1:4190953/150‑1 (MQ=255)

                                                                                                                                          aaa‑TGT‑GGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGCCa              >  1:2886470/1‑150 (MQ=255)

                                                                                                                                          aaa‑TGT‑GGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGCCa              >  1:2886836/1‑150 (MQ=255)

                                                                                                                                          aaa‑TGT‑GGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGCCa              >  1:5344716/1‑150 (MQ=255)

                                                                                                                                           aa‑TGT‑GGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGCCAg             >  1:4768655/1‑150 (MQ=255)

                                                                                                                                              tgt‑gGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGCCAGAc           >  2:1975060/1‑150 (MQ=255)

                                                                                                                                                   gt‑gCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGCCAGACCGcc       >  2:2059345/1‑150 (MQ=255)

                                                                                                                                                    t‑gCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGCCAGACCGccc      >  1:940494/1‑150 (MQ=255)

                                                                                                                                                        aGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGCCAGACCGCCCGTg   >  2:1793772/1‑150 (MQ=255)

                                                                                                                                                         gTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGCCAGACCGCCCGTgg  >  1:4544075/1‑150 (MQ=255)

                                                                                                                                                         gTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGCCAGACCGCCCGTgg  >  1:2150562/1‑150 (MQ=255)

                                                                                                                                                         gTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGCCAGACCGCCCGTgg  >  1:1727532/1‑150 (MQ=255)

TGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAAA‑TGT‑GGT‑GCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGCCAGACCGCCCGTGG  >  NC_002947/1397746‑1398044

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