BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,392,944	C→T	18.8%	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	18.8%	304.7 / 51.6	133	intergenic (+108/‑78)	PP_1214/ruvC	transcriptional regulator/crossover junction endodeoxyribonuclease RuvC
Reads supporting (aligned to +/- strand): ref base C (62/42); new base T (8/14); total (77/56)
Fisher's exact test for biased strand distribution p-value = 5.94e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

Only 100 of 133 total aligned reads displayed.

GCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGC‑GGGTAAA‑CCCGCTCCTACAGGTATTACG‑CC‑CATTTCAAATGTGGTGCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGCCAGACCGCCCG  >  NC_002947/1392798‑1393092

gCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGC‑GGGTAAA‑CCCGCTCCTACAGGTATTACG‑CC‑CATTTCAAATGTGGTGCAGt                                                                                                                                                    <  1:3484319/150‑1 (MQ=255)

 cccAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGC‑GGGTAAA‑CCCGCTCCTACAGGTATTACG‑CC‑CATTTCAAATGTGGTGCAGTa                                                                                                                                                   >  2:3832984/1‑150 (MQ=255)

  ccAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGC‑GGGTAAA‑CCCGCTCCTACAGGTATTACG‑CC‑CATTTCAAATGTGGTGCAGTAc                                                                                                                                                  <  2:1774636/150‑1 (MQ=255)

  ccAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGC‑GGGTAAA‑CCCGCTCCTACAGGTATTACG‑CC‑CATTTCAAATGTGGTGCAGTAc                                                                                                                                                  <  2:2368626/150‑1 (MQ=255)

    aGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGC‑GGGTAAA‑CCCGCTCCTACAGGTATTACG‑CC‑CATTTCAAATGTGGTGCAGTACCt                                                                                                                                                <  2:4399851/150‑1 (MQ=255)

    aGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGC‑GGGTAAA‑CCCGCTCCTACAGGTATTACG‑CC‑CATTTCAAATGTGGTGCAGTACCt                                                                                                                                                <  2:2217255/150‑1 (MQ=255)

     gATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGC‑GGGTAAA‑CCCGCTCCTACAGGTATTACG‑CC‑CATTTCAAATGTGGTGCAGTACCTg                                                                                                                                               <  1:2998880/150‑1 (MQ=255)

      aTTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGC‑GGGTAAA‑CCCGCTCCTACAGGTACTGCA‑CCACATTTGAAATG‑GGCGTAATACCTGt                                                                                                                                              <  2:1016199/150‑1 (MQ=21)

       tttCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGC‑GGGTAAA‑CCCGCTCCTACAGGTATTACG‑CC‑CATTTCAAATGTGGTGCAGTACCTGTa                                                                                                                                             <  2:185507/150‑1 (MQ=255)

        ttCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGC‑GGGTAAA‑CCCGCTCCTACAGTTATTACG‑CC‑CATTTCAAATGTGGTGCAGTACCTGTAg                                                                                                                                            <  1:1103546/150‑1 (MQ=255)

        ttCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGC‑GGGTAAA‑CCCGCTCCTACAGGTATTACG‑CC‑CATTTCAAATGTGGTGCAGTACCTGTAg                                                                                                                                            <  1:3211052/150‑1 (MQ=255)

        ttCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGC‑GGGTAAA‑CCCGCTCCTACAGGTATTACG‑CC‑CATTTCAAATGTGGTGCAGTACCTGTAg                                                                                                                                            <  1:3281105/150‑1 (MQ=255)

        ttCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGC‑GGGTAAA‑CCCGCTCCTACAGGTATTACG‑CC‑CATTTCAAATGTGGTGCAGTACCTGTAg                                                                                                                                            <  2:3786896/150‑1 (MQ=255)

         tCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGC‑GGGTAAA‑CCCGCTCCTACAGGTACTGCA‑CCACATTTGAAATG‑GGCGTAATACCTGTAgg                                                                                                                                           <  1:739310/150‑1 (MQ=18)

         tCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGC‑GGGTAAA‑CCCGCTCCTACAGGTACTGCA‑CCACATTTGAAATG‑GGCGTAATACCTGTAgg                                                                                                                                           <  1:1916310/150‑1 (MQ=18)

         tCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGC‑GGGTAAA‑CCCGCTCCTACAGGTACTGCA‑CCACATTTGAAATG‑GGCGTAATACCTGTAgg                                                                                                                                           <  1:1916230/150‑1 (MQ=18)

          ccGACGAGATCATGGAAACCCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGCGCCGGCCCTTTCGC‑GGGTAAA‑CCCGCTCCTACAGGTACTGCA‑CCACATTTGAAATG‑GGCGTAATACCTGTAGGa                                                                                                                                          <  1:285505/150‑1 (MQ=17)

          ccGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGC‑GGGTAAA‑CCCGCTCCTACAGGTATTACG‑CC‑CATTTCAAATGTGGTGCAGTACCTGTAGGa                                                                                                                                          <  1:2618769/150‑1 (MQ=255)

           cgacgaGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGC‑GGGTAAA‑CCCGCTCCTACAGGTATTACG‑CC‑CATTTCAAATGTGGTGCAGTACCGGTAGGAg                                                                                                                                         >  2:2081676/1‑150 (MQ=255)

           cgacgaGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGC‑GGGTAAA‑CCCGCTCCTACAGGTACTGCA‑CCACATTTGAAATG‑GGCGTAATACCTGTAGGAg                                                                                                                                         <  2:2454211/150‑1 (MQ=18)

            gacgaGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGC‑GGGTAAA‑CCCGCTCCTACAGGTATTACG‑CC‑CATTTCAAATGTGGTGCAGTACCTGTAGGAGc                                                                                                                                        <  2:3642651/150‑1 (MQ=255)

            gacgaGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGC‑GGGTAAA‑CCCGCTCCTACAGGTATTACG‑CC‑CATTTCAAATGTGGTGCAGTACCTGTAGGAGc                                                                                                                                        <  1:4128711/150‑1 (MQ=255)

              cgaGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGC‑GGGTAAA‑CCCGCTCCTACAGGTATTACG‑CC‑CATTTCAAATGTGGTGCAGTACCTGTAGGAGC‑gg                                                                                                                                     <  1:2688420/150‑1 (MQ=255)

              cgaGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGC‑GGGTAAA‑CCCGCTCCTACAGGTATTACG‑CC‑CATTTCAAATGTGGTGCAGTACCTGTAGGAGC‑gg                                                                                                                                     <  1:372196/150‑1 (MQ=255)

              cgaGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGC‑GGGTAAA‑CCCGCTCCTACAGGTATTACG‑CC‑CATTTCAAATGTGGTGCAGTACCTGTAGGAGC‑gg                                                                                                                                     <  1:372107/150‑1 (MQ=255)

              cgaGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGC‑GGGTAAA‑CCCGCTCCTACAGGTATTACG‑CC‑CATTTCAAATGTGGTGCAGTACCTGTAGGAGC‑gg                                                                                                                                     <  2:4160555/150‑1 (MQ=255)

              cgaGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGCGCCGGCCCTTTCGC‑GGGTAAA‑CCCGCTCCTACAGGTACTGCA‑CCACATTTGAAATG‑GGCGTAATACCTGTAGGAGC‑gg                                                                                                                                     <  2:3769324/150‑1 (MQ=17)

               gagaTCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGCTAGCGCCTGTGCCGGCCCTTTCGC‑GGGTAAA‑CCCGCTCCTCCAGGTATTACG‑CC‑CATTTCAAATGTGGTGCAGTACCTGTAGGAGC‑ggg                                                                                                                                    <  2:4005132/150‑1 (MQ=255)

               gagaTCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGC‑GGGTAAA‑CCCGCTCCTACAGGTATTACG‑CC‑CATTTCAAATGTGGTGCAGTACCTGTAGGAGC‑ggg                                                                                                                                    <  2:4002427/150‑1 (MQ=255)

               gagaTCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGC‑GGGTAAA‑CCCGCTCCTACAGGTATTACG‑CC‑CATTTCAAATGTGGTGCAGTACCTGTAGGAGC‑ggg                                                                                                                                    <  2:3629026/150‑1 (MQ=255)

                 gaTCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGC‑GGGTAAA‑CCCCCGCCAACACGTATTAAG‑CC‑CCAATAAAAAGTGGTGCAGGACCGGGAGGGGC‑ggggt                                                                                                                                  >  2:1907814/1‑148 (MQ=255)

                  aTCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGATTGATAGCGCCTGTGCCGGCCCTTTCGC‑GGGTAAA‑CCCGCTCCTACAGGTATTACG‑CC‑CATTTCAAATGTGGTGCAGTACCTGTAGGAGC‑GGGttt                                                                                                                                 <  1:2024226/150‑1 (MQ=255)

                    cATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGC‑GGGTAAA‑CCCGCTCCTACAGGTATTACG‑CC‑CATTTCAAATGTGGTGCAGTACCTGTAGGAGC‑GGGTTTAc                                                                                                                               <  1:490568/150‑1 (MQ=255)

                     aTGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGC‑GGGTAAA‑CCCGCTCCTACAGGTATTACG‑CC‑CAATTGAAATGTTGTGGCGTACCTGTAGGAGC‑GGGTTTAcc                                                                                                                              >  2:2667504/1‑150 (MQ=38)

                     aTGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCGTTCGC‑GGGTAAA‑ACCCCACCGACAGGTATTACG‑CC‑CAATTAAAATGTGGTGGCGTACCGGTGGGAGG‑GGGTTTaac                                                                                                                              >  2:2667444/1‑148 (MQ=255)

                       ggAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGC‑GGGTTAA‑ACCCCGCCGACAGGGATTAAG‑CC‑CATTTGAAATGTGGGGGAGTAACGGTGTGAGCGGGGTTTAccg                                                                                                                             >  2:1554695/1‑149 (MQ=255)

                       ggAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGC‑GGGTAAA‑CCCGCTCCTACAGGTACTGCA‑CCACATTTGAAATG‑GGCGTAATACCTGTAGGAGC‑GGGTTTACCcg                                                                                                                            <  1:3052465/150‑1 (MQ=17)

                       ggAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGC‑GGGTAAA‑CACCCTCCAACAGGTATTAAGCCC‑CATTACAAA‑GGGGGGCCGTACCAGTAGGAGC‑GGGTTTAACcg                                                                                                                            >  2:1554687/1‑150 (MQ=37)

                       ggAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGC‑GGGGAAA‑CCCGCGCCGACAGGTATTACG‑CC‑CCTTTCCAATTTGGTGGAGTACCGGGTGGTGC‑GGGTTTagacg                                                                                                                            >  1:2218242/1‑146 (MQ=255)

                        gAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGC‑GGGTAAA‑CCCGCTCCTACAGGTATTACG‑CC‑CATTTCAAATGTGGTGCAGTACCTGTAGGAGC‑GGGTTTACCcgc                                                                                                                           <  1:2732012/150‑1 (MQ=255)

                            aCCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGC‑GGGTAAA‑CCCGCTCCTACAGGTATTACG‑CC‑CATTTCAAATGTGGTGCAGTACCTGTAGGAGC‑GGGTTTACCCGCGaaa                                                                                                                       <  2:1411519/150‑1 (MQ=255)

                             ccTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGC‑GGGTAAA‑ACCGCGCCAACAGGTATTACG‑CC‑CATTTCAAATGTGGGGCAGTACCGGTAGGAGC‑GGGTTTACCCGCGAAAg                                                                                                                      >  2:711280/1‑150 (MQ=38)

                               tCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGC‑GGGTAAA‑CCCCCAACAACAGGTATTACG‑CC‑CCATTCAAATGGGGTGCCGTACCGGTAGGGGC‑GGGTTGAACCGCGCAAggg                                                                                                                    >  2:172372/1‑150 (MQ=255)

                                cGGCTAAGGTGATTGGCCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTCCGC‑GGGTAAA‑CCCGCTCCTACTGGTATTCCG‑CC‑CATTTCAACTGTGGTGCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGc                                                                                                                   <  2:3850003/150‑1 (MQ=37)

                                  gCTAAGGTGATTTGGCAGTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGGTAAA‑CCC‑CACCGACCGGGATTAAG‑CC‑CCATTCAAATGTGGTGCGCTACCGGGTGGAGG‑GGGGTTGCCCCCGACAAGGCCg                                                                                                                 >  2:3779046/1‑150 (MQ=255)

                                            tttGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGC‑GGGGAAA‑CCCCCGCCTACAGGTATTAAG‑CC‑CAATTCAAATGTGGTGCAGTTCCTGTTGGGGC‑GGGGTTACCCGCGGAACGGGCGGGGCCGGTTg                                                                                                       >  2:1948853/1‑150 (MQ=255)

                                                        tgCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGC‑GGGTAAA‑ACCGCTACAACAGGTATTAAG‑CC‑CATTTAAAATGTTGTGGAGTACATGTAGGAGC‑GGGTTTAACCGCGAAAGGGCCGGCGCAGGTTGTAGAGATTcaaa                                                                                           >  2:238531/1‑148 (MQ=25)

                                                          ccAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGC‑GGGTAAA‑CCCGCTCCTACAGGTATTACG‑CC‑CATTTCAAATGTGGTGCAGTACCTGTAg                                                                                                                                            >  2:291253/1‑100 (MQ=255)

                                                          ccAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGC‑GGGTAAA‑CCCGCTCCTACAGGTATTACG‑CC‑CATTTCAAATGTGGTGCAGTACCTGTAg                                                                                                                                            <  1:291253/100‑1 (MQ=255)

                                                            aaCTCTGTGCGTTGGATAGCGCGTGGGCCGCCCTTTTCG‑‑GGTTAAACCCCCCCCCCCCTGGTTTGCCG‑CC‑ATTTTCAATTTTGGTGCTGTACCGTTAGGTGG‑GGTTTGCCCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCAttt                                                                                       <  2:3223686/150‑1 (MQ=255)

                                                                 tGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGC‑GGGTGAA‑AACCCGCCAACACGTAGTACG‑CC‑CCATTCAAAGGTGGTGGCGTACCGGTTGGGGC‑GGGGGTTCCCCCGCAAGGGCCCGCGCCGGGTGTTGAGGTGGAACCttttttgg                                                                                  >  2:843466/1‑148 (MQ=255)

                                                                      gTTTGATAGCGCCTGTGCCGGCCCTTTCGC‑GGGTAAA‑ACCGCGCCTAAACGGCTTAAG‑CC‑CAGTTAAAATGTGGTGCCGGACCGGTAGGAGC‑GGGTTTAACCGCGCAAGGGCCGGCGCGGGGTGGAGAGATGCCTGATTTTTTCGACggg                                                                             >  2:3095040/1‑150 (MQ=255)

                                                                           aTAGCGCCTGTGCCGGCCCATTCGC‑GGGTAAA‑ACCGCACCAAAAGGTGTTACC‑CC‑CCATTCAAAACTGGGGCCGTACCGGTGGGAGC‑GGGGTTACCCCCGCAAGGGCCGGCGCCGTTTGGAGCGGAGCCACATTTTTTAGAAGGGCtttt                                                                        >  2:3128043/1‑149 (MQ=255)

                                                                                   tgtgCCGGCCCTTTCGC‑GGGTAAA‑CCCGCTCCTACAGGTATTACG‑CC‑CATTTCAAATGTGGTGCAGTACCTGTAGGAGC‑gg                                                                                                                                     >  1:508792/1‑81 (MQ=255)

                                                                                        cGGCCCTTTCGG‑GGCAAAA‑CCCGCACCTCCTGTTATTACC‑CC‑CTTTTCACTTGTGGTGCCTAACCTTTCGGCGG‑GGTTTAACCCGCGAACGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCt                                                           <  2:1990483/150‑1 (MQ=255)

                                                                                             cTTTCGC‑GGGTAAA‑ACCGCTCCTACAGGTATTAAG‑CC‑CCTTTCAAATGTGGTGCAGTACCTGTATGAGC‑GGGTTTAACCGCGAAAGGGCCGGCGCAGTTTTTAGAGATGCAACAATTTTTATAAGGGCTTTATTAATCTGATTCTATGTa                                                      >  2:2073681/1‑150 (MQ=255)

                                                                                             cTTTCGC‑GGGTAAA‑ACCGCTCCTACAGGTATTAAG‑CC‑CATTTCAAATGTGGTGCAGTACCTGTAGGAGC‑GGGTTTAACCGCGAAAGGGCCGGCGCAGTTTGTATAGATGCATCAATTTTTCGACGGGCTTTATTACTATGATTCTAGGTa                                                      >  2:3124403/1‑150 (MQ=37)

                                                                                                tCGC‑GGGTAAA‑CCCGCTCCTACAGGTATTACG‑CC‑CATTTCAAATGTGGTGCCGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCg                                                   >  2:266095/1‑150 (MQ=255)

                                                                                                       taaaa‑CCCGCTCCTACTGTTATTACG‑CC‑AATTTCAATTGTGGTGCTGTACCTGTTGGTGG‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCg                                             <  2:2218242/148‑1 (MQ=255)

                                                                                                        tAAA‑CCCGCTCCTACAGGTATTACG‑CC‑CATTTCAAATGTGGTGCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCgg                                            >  1:4235469/1‑150 (MQ=255)

                                                                                                          aa‑CCCGCTCCTACAGGTATTACG‑CC‑CATTTCAAATGTGGTGCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCt                                          >  2:744205/1‑150 (MQ=255)

                                                                                                           a‑CCCGCTCCTACAGGTATTACG‑CC‑CATTTCAAATGTGGTGCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTc                                         >  1:176582/1‑150 (MQ=255)

                                                                                                             cccGCTCCTACAGGTATTACG‑CC‑CATTTCAAATGTGGTGCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTcg                                        >  1:3045565/1‑150 (MQ=255)

                                                                                                             cccGCTCCTACAGGTACTGCA‑CCACATTTGAAATG‑GGCGTAATACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTcg                                        >  1:561863/1‑150 (MQ=18)

                                                                                                              ccGCTCCTACAGGTATTACG‑CC‑CATTTCACATGTGGTGCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGGTTGTAGAGATGCATCATTTTTTCGACGGGCGTTATGACTCTGATTCTAGGTATCGACCCCGGCTcgc                                       >  2:3996063/1‑150 (MQ=255)

                                                                                                              ccGCTCCTACAGGTATTACG‑CC‑CATTTCAAATGTGGTGCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGCGCTGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTcgc                                       >  2:3032821/1‑150 (MQ=255)

                                                                                                              ccGCTCCTACAGGTATTACG‑CC‑CATTTCAAATGTGGTGCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTcgc                                       >  1:4409362/1‑150 (MQ=255)

                                                                                                                gCTCCTACAGGTACTGCA‑CCACATTTGAAATG‑GGCGTAATACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCgcgc                                     >  2:478922/1‑150 (MQ=18)

                                                                                                                 cTCCTACAGGTACTGCA‑CCACATTTGAAATG‑GGCGTAATACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCa                                    >  1:3821846/1‑150 (MQ=18)

                                                                                                                  tCCTACAGGTATTACG‑CC‑CATTTCAAATGTGGTGCAGTACCTGTAGGa                                                                                                                                          >  2:322888/1‑48 (MQ=255)

                                                                                                                  tCCTACAGGTATTACG‑CC‑CATTTCAAATGTGGTGCAGTACCTGTAGGa                                                                                                                                          <  1:322888/48‑1 (MQ=255)

                                                                                                                  tCCTACAGGTACTGCA‑CCACATTTGAAATG‑GGCGTAATACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTGTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGGATCGACCCCGGCTCGCGCAt                                   >  1:2026848/1‑150 (MQ=18)

                                                                                                                  tCCTACAGGTACTGCA‑CCACATTTGAAATG‑GGCGTAATACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCACAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCAt                                   >  2:4261762/1‑150 (MQ=18)

                                                                                                                   ccTACAGGTATTACG‑CC‑CATTTCAAATGTGGTGCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATc                                  >  1:3565190/1‑150 (MQ=255)

                                                                                                                    cTACAGGTACTGCA‑CCACATTTGAAATG‑GGCGTAATACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTCGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCa                                 >  1:1897715/1‑150 (MQ=18)

                                                                                                                     tACAGGTATTACG‑CC‑CATTTCAAATGTGGTGCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGGAGAGATGCATCATTTTTTCGACGGGCGTTATGACTCTGAGTCTAGGTATCGACCCCGGCTCGCGCATCAc                                >  2:4124303/1‑150 (MQ=255)

                                                                                                                         ggTATTACG‑CC‑CATTTCAAATGTGGTGCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGc                            >  1:1038516/1‑150 (MQ=255)

                                                                                                                          gTATTACG‑CC‑CATTTCAAATGTGGTGCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCt                           >  1:3694618/1‑150 (MQ=255)

                                                                                                                               aCG‑CC‑CATTTCAAATGTGGTGCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGc                      >  1:241693/1‑150 (MQ=255)

                                                                                                                                 g‑CC‑CATTTCAAATGTGGTGCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGt                    >  2:2821119/1‑150 (MQ=255)

                                                                                                                                 g‑CC‑CATTTCAAATGTGGTGCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGt                    >  1:2877975/1‑150 (MQ=255)

                                                                                                                                 g‑CC‑CATTTCAAATGTGGTGCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGt                    >  2:3719808/1‑150 (MQ=255)

                                                                                                                                 g‑CC‑CATTTCAAATGTGGTGCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGt                    >  1:1941954/1‑150 (MQ=255)

                                                                                                                                 g‑CC‑CATTTCAAATGTGGTGCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGt                    >  1:479090/1‑150 (MQ=255)

                                                                                                                                    c‑cATTTCAAATGTGGTGCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTgg                  >  2:1577085/1‑150 (MQ=255)

                                                                                                                                    c‑cATTTCAAAGGTGGTGCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTgg                  >  2:144029/1‑150 (MQ=255)

                                                                                                                                      cATTTCAAATGTGGTGCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGt                 >  1:1197523/1‑150 (MQ=255)

                                                                                                                                      cATTTCAAATGTGGTGCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGt                 >  2:1997444/1‑150 (MQ=255)

                                                                                                                                       aTTTCAAATGTGGTGCAGTACCTGTAAGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTa                >  2:325698/1‑150 (MQ=255)

                                                                                                                                        tttCAAATGTGGTGCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTAc               >  2:3643439/1‑150 (MQ=255)

                                                                                                                                         ttCAAATGTGGTGCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACg              >  1:2240641/1‑150 (MQ=255)

                                                                                                                                           cAAATGTGGTGCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGcc            >  1:3464222/1‑150 (MQ=255)

                                                                                                                                            aaaTGTGGTGCAGTACCTGTAGGAGC‑GGGTTTAGCCGCGAGAGGGCCGGCGCAGTTTGTCGAGATGCATCATTTTTTCGACGGGGTTTATGACTCTGATTCTAGGTATCGGCCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGCCa           >  2:171391/1‑150 (MQ=255)

                                                                                                                                            aaaTGTGGTGCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGCCa           >  2:2550106/1‑150 (MQ=255)

                                                                                                                                            aaaTGTGGTGCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGCCa           >  1:128268/1‑150 (MQ=255)

                                                                                                                                            aaaTGTGGTGCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACACCGGCTCGCGCATCACCGGCTATGGCGTGGTACGCCa           >  2:171393/1‑150 (MQ=255)

                                                                                                                                              aTGTGGTGCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGCCAGa         <  1:1577085/150‑1 (MQ=255)

                                                                                                                                              aTGTGGTGCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGCCAGa         <  2:4409362/150‑1 (MQ=255)

                                                                                                                                                 tggtgCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGCCAGACCg      >  1:1694382/1‑150 (MQ=255)

                                                                                                                                                   gtgCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCCGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGCCAGgccgcc    >  2:2019929/1‑150 (MQ=255)

GCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGC‑GGGTAAA‑CCCGCTCCTACAGGTATTACG‑CC‑CATTTCAAATGTGGTGCAGTACCTGTAGGAGC‑GGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGCCAGACCGCCCG  >  NC_002947/1392798‑1393092

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