BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,392,944	C→T	13.1%	intergenic (+108/‑78)	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,392,944	0	C	T	13.1%	512.5 / 55.1	201	intergenic (+108/‑78)	PP_1214/ruvC	transcriptional regulator/crossover junction endodeoxyribonuclease RuvC
Reads supporting (aligned to +/- strand): ref base C (101/69); new base T (9/16); total (114/87)
Fisher's exact test for biased strand distribution p-value = 3.20e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.82e-02

Only 100 of 206 total aligned reads displayed.

AATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTG‑CCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTAC‑GCC‑CA‑TTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGCCAGACCGCCCGT  >  NC_002947/1392795‑1393093

aaTGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTG‑CCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTAC‑GCC‑CA‑TTTCAAATGTGGTGc                                                                                                                                                        <  2:8846411/150‑1 (MQ=255)

 aTGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTG‑CCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTAC‑GCC‑CA‑TTTCAAATGTGGTGCa                                                                                                                                                       <  1:6313288/150‑1 (MQ=255)

 aTGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTG‑CCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTAC‑GCC‑CA‑TTTCAAATGTGGTGCa                                                                                                                                                       <  2:3713488/150‑1 (MQ=255)

 aTGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTG‑CCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTAC‑GCC‑CA‑TTTCAAATGTGGTGCa                                                                                                                                                       <  2:4485213/150‑1 (MQ=255)

     ccAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTG‑CCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTAC‑GCC‑CA‑TTTCAAATGTGGTGCAGTAc                                                                                                                                                   <  1:6999786/150‑1 (MQ=255)

      cAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTG‑CCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCATTTCGCGGGTAAACCCGCTCCTACAGGTATTAC‑GCC‑CA‑TTTCAAATGTGGTGCAGTAcc                                                                                                                                                  <  2:373847/150‑1 (MQ=255)

       aGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTG‑CCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTAC‑GCC‑CA‑TTTCAAATGTGGTGCAGTACCt                                                                                                                                                 <  1:4588784/150‑1 (MQ=255)

         aTTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTG‑CCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTACTGC‑ACCACA‑TTTGAAATG‑GGCGTAATACCTGt                                                                                                                                               <  2:6470613/150‑1 (MQ=21)

          tttCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTG‑CCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTAC‑GCC‑CA‑TTTCAAATGTGGTGCAGTACCTGTa                                                                                                                                              <  1:7542474/150‑1 (MQ=255)

          tttCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTG‑CCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTACTGC‑ACCACA‑TTTGAAATG‑GGCGTAATACCTGTa                                                                                                                                              <  2:6357651/150‑1 (MQ=21)

          tttCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTG‑CCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTACTGC‑ACCACA‑TTTGAAATG‑GGCGTAATACCTGTa                                                                                                                                              <  1:5464451/150‑1 (MQ=21)

           ttCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTG‑CCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTAC‑GCC‑CA‑TTTCAAATGTGGTGCAGTACCTGTAg                                                                                                                                             <  1:2615634/150‑1 (MQ=255)

           ttCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTG‑CCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTAC‑GCC‑CA‑TTTCAAATGTGGTGCAGTACCTGTAg                                                                                                                                             <  1:6540923/150‑1 (MQ=255)

           ttCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTG‑CCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTAC‑GCC‑CA‑TTTCAAATGTGGTGCAGTACCTGTAg                                                                                                                                             <  2:6805808/150‑1 (MQ=255)

            tCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTG‑CCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTAC‑GCC‑CA‑TTTCAAATGTGGTGCAGTACCTGTAgg                                                                                                                                            <  1:6254066/150‑1 (MQ=255)

               gacgaGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTG‑CCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTAC‑GCC‑CA‑TTTCAAATGTGGTGCAGTACCTGTAGGAGc                                                                                                                                         <  2:6441871/150‑1 (MQ=255)

               gacgaGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTG‑CCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTACTGC‑ACCACA‑TTTGAAATG‑GGCGTAATACCTGTAGGAGc                                                                                                                                         <  1:2944319/150‑1 (MQ=18)

               gacgaGAACATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTG‑CCAACTCTGAGCGTTTGATAGCGCCTGTGCCGCCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTAC‑GCC‑CA‑TTTCAAATGTGGTGCAGTACCTGTAGGAGc                                                                                                                                         <  2:2555391/150‑1 (MQ=255)

                 cgaGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTG‑CCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTAC‑GCC‑CA‑TTTCAAATGTGGTGCAGTACCTGTAGGAGCgg                                                                                                                                       <  1:1703740/150‑1 (MQ=255)

                 cgaGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTG‑CCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTAC‑GCC‑CA‑TTTCAAATGTGGTGCAGTACCTGTAGGAGCgg                                                                                                                                       <  1:767109/150‑1 (MQ=255)

                  gagaTCATGGAAAACCTCGGCTAAGGTGATTTGGCTTTTTGTG‑CCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTACTGC‑ACCACA‑TTTGAAATG‑GGCGTAATACCTGTAGGAGCggg                                                                                                                                      <  1:3869268/150‑1 (MQ=18)

                  gagaTCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTG‑CCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTAC‑GCC‑CA‑TTTCAAATGTGGTGCAGTACCTGTAGGAGCggg                                                                                                                                      <  1:7313638/150‑1 (MQ=255)

                  gagaTCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTG‑CCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTAC‑GCC‑CA‑TTTCAAATGTGGTGCAGTACCTGTAGGAGCggg                                                                                                                                      <  1:3703052/150‑1 (MQ=255)

                    gaTCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTG‑CCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTAC‑GCC‑CA‑TTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGtt                                                                                                                                    <  1:4035994/150‑1 (MQ=255)

                    gaTCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTG‑CCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTAC‑GCC‑CA‑GTTCAAATGTGGTGCAGTACCTGGAGGAGCGGGtt                                                                                                                                    >  2:9023992/1‑150 (MQ=255)

                        aTGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTG‑CCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTAC‑GCC‑CA‑TTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTA‑cc                                                                                                                               <  2:5004989/150‑1 (MQ=255)

                        aTGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTG‑CCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTAC‑GCC‑CA‑TTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTA‑cc                                                                                                                               <  2:544691/150‑1 (MQ=255)

                        aTGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTG‑CCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAAAACCCTCCAAAAGGTATTAA‑GCC‑CA‑ATTAAAAATTTGTGCAGTAACTGTTTGTGCGGGTTTaac                                                                                                                                >  1:6794608/1‑148 (MQ=255)

                         tGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTG‑CCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAAACCCCTCCTACAGGTATTAC‑GCC‑CA‑TTTCAAATGTGGTGCAGTACCGGTAGGAGCGGGTTTT‑ccc                                                                                                                              >  2:4515610/1‑150 (MQ=255)

                                   cGGCTATGGTGATTGGGCTTTTTGTG‑CCAACTCTGAGCGTTTGATAGCGCGTGTGCGGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTACTGC‑ACCACA‑TTTTAACTG‑GGCTTAATACCTGTAGGCGCGGGTTTA‑CCCGCGAAAGGGc                                                                                                                    <  2:4759867/150‑1 (MQ=14)

                                    ggCTATGGTGTTTTGGCATTTGGTG‑CCAACTCTGAGGGTTTGATAGCGCCTGTGCGGGCCCTTTCGGGGTTAAACCCGCTCCTACAGGTATTAC‑GCC‑CA‑TTTCAAATGTGGTGCAGTACCTGTAGGAGGGGGTTTA‑CCCGCGAAAGGGcc                                                                                                                   <  2:5773311/150‑1 (MQ=37)

                                     gCTATGGTTTTTTGGCTTTTTGTG‑CCAACTGTTTGCTTTTGATAGCGCCTGTGCCGGCCCTTTCGGGGGTAAACCCGCTCCTACTGTTTTTAC‑GCC‑AA‑TTTCAATTTTGGTGCAGTACCTGTAGGAGCGGGTTGA‑CCCGCGAAAGGGCCg                                                                                                                  <  2:4795799/150‑1 (MQ=255)

                                       tAAGGTGATTTGGCATTTTGTG‑CCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTAC‑GCC‑CA‑TTTCAAATGTGGTGCAGTACCGGTAGGAGCGGGTGTA‑ACCCCGAAAGGGCCGgc                                                                                                                >  2:4483738/1‑150 (MQ=255)

                                        aaGGTGATTTGGCATTTTGTG‑CCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTAC‑GCC‑CA‑TTTCAAAGGTGGTGCAGTACCTGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCGgcg                                                                                                               >  2:7602112/1‑150 (MQ=255)

                                               tttGGGATTTTTTGCCCATCTGTTAGCTTTTGTTCGCGCCTGTCCGGGCCCTTTCGGGTTTAAACCCGATCCAACCGTTATTAC‑GCA‑AT‑TTTCAAATGGGGTGCTCACCCTTTTGGAGCGGTTTTA‑CCCGCGAAAGGGCCGGCGCAGttt                                                                                                         <  2:7867758/150‑1 (MQ=255)

                                                     aTTTTGTG‑CCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAAACCGCTCCAAAAGGTATTAC‑GCC‑CC‑ATTCAAATGTGGTGCAGTACCGGTAGGAGCGGGTTTA‑CCCGCGAAAGGGGCGGCGCAGTTTGTAgaga                                                                                                  >  2:2907118/1‑150 (MQ=37)

                                                           tg‑CCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCAACAGGTATTAA‑GCC‑CA‑ATTCAAATGTGGTGCAGTACCCGGAGGAGCGGGTTTA‑CCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGcatc                                                                                            >  2:4400786/1‑150 (MQ=38)

                                                               cAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCACCTACAGGTATTAC‑GCC‑CA‑TTTAAAAAGTGGTGCAGTACCAGTAGGAGCCGGTTTA‑CCCCCGGAAGGGCCGGCGCAGGTTGTAGAGATGCATCCtt                                                                                         >  2:3690749/1‑150 (MQ=37)

                                                                 aCTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAAACCCCTACTAAAGGTATTAA‑GCC‑CA‑ATTAAAATGTTGTGCAGTACCAGTAGGAGCGGGTTTA‑CCCCCGCAAGGGCCCGCGCAGTTTGTATATATGCAAAAtttt                                                                                       >  2:4324281/1‑150 (MQ=255)

                                                                           tttGATCGCGCCTGTGCCGGCCCTTTCGCGGTTAAACCCGCTCCTCCCGTTTTTAC‑GCC‑CA‑TTTAAAATGTGGTGCAGTACCTGTAGGAGCGGTTTTA‑CCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGc                                                                             <  2:2887896/150‑1 (MQ=37)

                                                                               aTAGCGCCTGTGCCGGCCCTTTCGCGGGTAAAACCCCACCAAAACGGATTAA‑GCC‑CC‑AATCAAATTGGGGGCCGTAGCCGTCGGGGCGGGGTTA‑ACCCCCCAAGGGGCCGCGCCGGTTGTATAGATGCCACATTTTTTAGAAGGGCTTTa                                                                         >  2:4242223/1‑150 (MQ=255)

                                                                               aTAGCGCCTGTGCCGGCCCTTTCGCGGGGAAAACCCCGCCAACACGTATTAC‑GCC‑CC‑AATGAAATGTGGTGGCGTACCGGGAGGGGCGGGGTTA‑ACCCCGAAACGGCCCGCGCGGGTGGTTGGGATGGCTCATTTTTTGGACGGGCgttt                                                                         >  2:6216785/1‑149 (MQ=255)

                                                                                tggcgcGTGTGCGGGCCGTTGCG‑GGGTTAACCCGGCCCCGCCGGTTTTACGGCC‑CT‑TTTATATTGTGGTGCCTTACGTGTGGGTGGGGTTTTC‑CCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTAt                                                                        <  2:4167978/148‑1 (MQ=255)

                                                                                  gcggCTGTGCCGCCCCTTTCGCGGTTAACCCCGATCCCACCGGTCTCAC‑CCC‑CA‑TTTCACATGTGGTGCAGTACCTGTAGGTGGGGTTTTA‑CCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGa                                                                      <  2:8768650/150‑1 (MQ=37)

                                                                                     ccTGTGCCGGCCCTTTCGCGGGTAAACCCCCTCCTACAGGTATTAA‑GCC‑CC‑TTTGAAATGGGGGGCAGTAACGGTAGGAGCGGGGTTA‑ACCGCGAAAGGGCCCGCGCCGGTTGTAGAGAAGCCTCAATTTTTAGACGGGGGTTATGActc                                                                   >  2:8121291/1‑150 (MQ=18)

                                                                                         tgCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTAC‑GCC‑CA‑TTTCAAATGTGGTGCAGTACCTGTAGTATCGGGTTTA‑CCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGAt                                                               <  2:6029560/150‑1 (MQ=255)

                                                                                             ggCCCTTGCGCGGGTAACCCCGCTCCTACAGGTATTAC‑GCC‑CA‑TTTCAATTGTGGTGCCGTACCTGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTa                                                           <  2:8970281/150‑1 (MQ=255)

                                                                                             gccccTTCCGGGGTTAACCCCGCCCCTACCGTTA‑TAC‑CCC‑CATTTTCAATGGTGGTGCTGTACGTTTTGGCGGGGTTTTA‑CCCGCGAACGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTa                                                           <  2:3991863/148‑1 (MQ=255)

                                                                                              gCCCTTTCGGGGTTAAACCCGCCCCTTCAGTTATTAC‑GCC‑AT‑TTTCAACTGTGGTTCAGTACCTGTTGGAGCGGGTTTA‑CCCGCGAATGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAg                                                          <  2:5833729/150‑1 (MQ=38)

                                                                                              gCCCTTCGGCGGTAAACCCCGCTCCTACCGTTATTCC‑CCC‑CA‑TTTCAACGGTGGTGCTGAACCTGTTGGCGCGGGTTTA‑CCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAg                                                          <  1:8362381/150‑1 (MQ=255)

                                                                                                ccTTTAGCGGGTAAAACCGCACCAAAAGGTAATAA‑GCC‑CA‑TTTCAAATGTGGTGCAGTAACTGTAGGAGCCGGTTTA‑ACCCCGAAAGGGCCGGCGCAGTTTGTAGAGATTCATCATTTTTTCGACGGGCGTTAAGACAATGTTTCGAGGt                                                        >  2:3786890/1‑150 (MQ=255)

                                                                                                      gcgGGTACCCCCGCTCCCCCCGGTATTCC‑GCC‑CT‑GTTCAATTGTGGTGCTGTACCTGTTGGCGGGGTTTTA‑CCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGAc                                                  <  2:7737907/150‑1 (MQ=255)

                                                                                                      gcgGGTAAACCCGCTCCTACAGGTATTAC‑GCC‑CA‑TTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGAc                                                  >  2:6891310/1‑150 (MQ=255)

                                                                                                        gggTAAACCCGCTCCTACAGGTACTGC‑ACCACA‑TTTGAAATG‑GGCGTAATACCTGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGAccc                                                >  1:3357009/1‑150 (MQ=17)

                                                                                                           tAAACCCGCTCCTACAGGTATTAC‑GCC‑CA‑TTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCgg                                             >  1:6156676/1‑150 (MQ=255)

                                                                                                              aCCCGCTCCTACAGGTATTAC‑GCC‑CA‑TTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCGGGCGCAGTTTGTAGAGATGCATCATTTTTTCGAGGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTc                                          >  2:733015/1‑150 (MQ=255)

                                                                                                                ccGCTCCTACAGGTATTAC‑GCC‑CA‑TTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTA‑CCCGCGACAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTcgc                                        >  1:3020454/1‑150 (MQ=255)

                                                                                                                ccGCTCCTACAGGTATTAC‑GCC‑CA‑TTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTA‑CCCGCGACAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTcgc                                        >  2:7003516/1‑150 (MQ=255)

                                                                                                                ccGCTCCTACAGGTATTAC‑GCC‑CA‑TTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTcgc                                        >  1:2553293/1‑150 (MQ=255)

                                                                                                                ccGCTCCTACAGGTACTGC‑ACCACA‑TTTGAAATG‑GGCGTAATACCTGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTcgc                                        >  2:3534747/1‑150 (MQ=18)

                                                                                                                  gCTCCTACAGGTATTAC‑GCC‑CA‑TTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCgcgc                                      >  1:8871883/1‑150 (MQ=255)

                                                                                                                      cTACAGGTACTGC‑ACCACA‑TTTGAAATG‑GGCGTAATACCTGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTAt                                                                        >  1:1307246/1‑112 (MQ=11)

                                                                                                                      cTACAGGTACTGC‑ACCACA‑TTTGAAATG‑GGCGTAATACCTGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTAt                                                                        <  2:1307246/112‑1 (MQ=11)

                                                                                                                       tACAGGTATTAC‑GCC‑CA‑TTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCAc                                 >  1:4506199/1‑150 (MQ=255)

                                                                                                                       tACAGGTATTAC‑GCC‑CA‑TTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCAc                                 >  2:5127937/1‑150 (MQ=255)

                                                                                                                       tACAGGTATTAC‑GCC‑CA‑TTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCAGCAc                                 >  2:5203434/1‑150 (MQ=255)

                                                                                                                       tACAGGTATTAC‑GCC‑CA‑TTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATGCAAGGTATCGACCCCGGCTCGCGCATCAc                                 >  1:2230271/1‑150 (MQ=255)

                                                                                                                       tACAGGTACTGC‑ACCACA‑TTTGAAATG‑GGCGTAATACCTGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCAc                                 >  2:7232081/1‑150 (MQ=21)

                                                                                                                       tACAGGTACTGC‑ACCACA‑TTTGAAATG‑GGCGTAATACCTGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCAc                                 >  1:7886107/1‑150 (MQ=21)

                                                                                                                        aCAGGTATTAC‑GCC‑CA‑TTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCAcc                                >  2:3840291/1‑150 (MQ=255)

                                                                                                                        aCAGGTATTAC‑GCC‑CA‑TTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCAcc                                >  1:7359071/1‑150 (MQ=255)

                                                                                                                          aGGTATTAC‑GCC‑CA‑TTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCgg                              >  2:1069806/1‑150 (MQ=255)

                                                                                                                          aGGTATTAC‑GCC‑CA‑TTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCgg                              >  2:5618554/1‑150 (MQ=255)

                                                                                                                            gTATTAC‑GCC‑CA‑TTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCt                            >  1:1105/1‑150 (MQ=255)

                                                                                                                               tcaC‑GCC‑CT‑TTACACCTGGGGTGCAGTCCCTGTAGCCG‑GGGTTAACCCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATg                         <  1:5208359/148‑1 (MQ=255)

                                                                                                                                 aC‑GCC‑CT‑TTTCAATTGTGGTGATGTACCGGTTGGAGGGGGTTTA‑CCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGc                       <  2:670963/150‑1 (MQ=255)

                                                                                                                                    gCC‑CA‑TTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTa                                                           <  1:838999/112‑1 (MQ=255)

                                                                                                                                    gCC‑CA‑TTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTa                                                           >  2:838999/1‑112 (MQ=255)

                                                                                                                                    gCC‑CA‑TTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGt                     >  1:6900087/1‑150 (MQ=255)

                                                                                                                                     cc‑cA‑TTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTg                    <  1:1069806/150‑1 (MQ=255)

                                                                                                                                      c‑cA‑TTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTgg                   >  1:8555458/1‑150 (MQ=255)

                                                                                                                                      c‑cA‑TTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTgg                   >  2:2353434/1‑150 (MQ=255)

                                                                                                                                      c‑cA‑TTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTgg                   >  2:2009741/1‑150 (MQ=255)

                                                                                                                                      c‑cA‑TTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTgg                   >  2:8073981/1‑150 (MQ=255)

                                                                                                                                        cA‑TTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTA‑ACCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATAGACACCGGCTCGCGCATCACCGGCTATGGCGTGGt                  >  2:8152749/1‑150 (MQ=255)

                                                                                                                                         a‑TTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTa                 >  1:7425311/1‑150 (MQ=255)

                                                                                                                                           tttCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTAc                >  2:8035081/1‑150 (MQ=255)

                                                                                                                                           tttCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTAc                >  2:5277580/1‑150 (MQ=255)

                                                                                                                                            ttCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACg               >  2:456987/1‑150 (MQ=255)

                                                                                                                                            ttCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACg               >  2:5111542/1‑150 (MQ=255)

                                                                                                                                               aaaTGTGGTGCAGTACCTGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGCCa            >  2:1940900/1‑150 (MQ=255)

                                                                                                                                               aaaTGTGGTGCAGTACCTGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGCCa            >  1:6871272/1‑150 (MQ=255)

                                                                                                                                               aaaTGTGGTGCAGTACCTGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGCCa            >  2:8376434/1‑150 (MQ=255)

                                                                                                                                               aaaTGTGGTGCAGTACCTGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGCCa            >  1:349649/1‑150 (MQ=255)

                                                                                                                                               aaaTGTGGTGCAGTACCTGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGCCa            >  1:8685612/1‑150 (MQ=255)

                                                                                                                                               aaaTGTGGTGCAGTACCTGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGGATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGCCa            >  1:4921687/1‑150 (MQ=255)

                                                                                                                                                       tgCAGTACCTGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGCCAGACCGccc    >  1:440706/1‑150 (MQ=255)

                                                                                                                                                       tgCAGTACCTGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGCCAGACCGccc    >  1:8266781/1‑150 (MQ=255)

                                                                                                                                                        gCAGTACCTGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGCCAGACCGCCCg   >  1:5615240/1‑150 (MQ=255)

                                                                                                                                                        gCAGTACCTGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCAGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGAGTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGCCAGACCGCCCg   >  1:4739765/1‑150 (MQ=255)

AATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTG‑CCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTAC‑GCC‑CA‑TTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGCCAGACCGCCCGT  >  NC_002947/1392795‑1393093

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