BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,952,787	C→G	11.5%	intergenic (‑56/+114)	PP_1748 ← / ← PP_1749	hypothetical protein/acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,952,787	0	C	G	11.5%	231.3 / 3.4	101	intergenic (‑56/+114)	PP_1748/PP_1749	hypothetical protein/acetyltransferase
Reads supporting (aligned to +/- strand): ref base C (21/65); new base G (8/5); total (31/70)
Fisher's exact test for biased strand distribution p-value = 1.77e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 107 total aligned reads displayed.

TGTCGGTGAAGCCGGTGGGGCTGGGGATGGCGAGCATTTCCAGCAGCACGCGTTTGAGATAGTCGAGATCGGGTTCGGGG‑AGTCGTTCGGACATGTCCAGGCCTCCGGTCGGG‑TTAAAATTTTGCAGTCTTGTGGGAGCGGCCTTGTGT‑CGCGAAAGGGCCGCAAAGCGGCCCC‑GGCAATATCGGCAGCGAAGCTGAGACCCTAGGGGCGCTTCGCACCCCATTCGCGACACAAGGCCGCTCCCACAAAAAGCCGAGCAGTTGTCAGGCCAACGGCCGGCTATGGGGAAACAGCAGG  >  NC_002947/1952638‑1952934

tgtCGGTGAAGCCGGTGGGGCTGGGGATGGCGAGCATTTCCAGCAGCACGCGTTTGAGATAGTCGAGATCGGGTTCGGGG‑AGTCGTTCGGACATGTCCAGGCCTCCGGTCGGG‑TTAAAATTTTGCAGTCTTGTGGGAGCGGCCTTGTGT‑c                                                                                                                                                      <  2:3482355/150‑1 (MQ=255)

  tCGGTGAAGCCGGTGGGGCTGGGGATGGCGAGCATTTCCAGCAGCACGCGTTTGAGATAGTCGAGATCGGGTTCGGGG‑AGTCGTTCGGACATGTCCAGGCCTCCGGTCGGG‑TTAAAATTTTGCAGTCTTGTGGGAGCGGCCTTGTGT‑cgc                                                                                                                                                    <  1:2304133/150‑1 (MQ=255)

  tCGGTGAAGCCGGTGGGGCTGGGGATGGCGAGCATTTCCAGCAGCACGCGTTTGAGATAGTCGAGATCGGGTTCGGGG‑AGTCGTTCGGACATGTCCAGGCCTCCGGTCGGG‑TTAAAATTTTGCAGTCTTGTGGGAGCGGCCTTGTGT‑cgc                                                                                                                                                    <  1:2831849/150‑1 (MQ=255)

  tCGGTGAAGCCGGTGGGGCTGGGGATGGCGAGCATTTCCAGCAGCACGCGTTTGAGATAGTCGAGATCGGGTTCGGGG‑AGTCGTTCGGACATGTCCAGGCCTCCGGTCGGG‑TTAAAATTTTGCAGTCTTGTGGGAGCGGCCTTGTGT‑cgc                                                                                                                                                    <  1:504013/150‑1 (MQ=255)

     gTGAAGCCGGTGGGGCTGGGGATGGCGAGCATTTCCAGCAGCACGCGTTTGAGATAGTCGAGATCGGGTTCGGGG‑AGTCGTTCGGACATGTCCAGGCCTCCGGTCGGG‑TTAAAATTTTGCAGTCTTGTGGGAGCGGCCTTGTGT‑CGCGaa                                                                                                                                                 <  2:813352/150‑1 (MQ=255)

     gTGAAGCCGGTGGGGCTGGGGATGGCGAGCATGTCCAGCAGCACGCGTTTGAGATAGTCGAGATCGGGTTCGGGG‑AGTCGTTCGGACATGTCCAGGCCTCCGGTCGGG‑TTAAAATTTTGCAGTCTTGTGGGAGCGGCCTTGTGT‑CGCGaa                                                                                                                                                 <  2:4006028/150‑1 (MQ=255)

      tGAAGCCGGTGGGGCTGGGGATGGCGAGCATTTCCAGCAGCACGCGTTTGAGATAGTCGAGATCGGGTTCGGGG‑AGTCGTTCGGACATGTCCAGGCCTCCGGTCGGG‑TTAAAATTTTGCAGTCTTGTGGGAGCGGCCTTGTGT‑CGCGaaa                                                                                                                                                <  2:4742029/150‑1 (MQ=255)

       gAAGCCGGTGGGGCTGGGGATGGCGAGCATTTCCAGCAGCACGCGTTTGAGATAGTCGAGATCGGGTTCGGGG‑AGTCGTTCGGACATGTCCAGGCCTCCGGTCGGG‑TTAAAATTTTGCAGTCTTGTGGGAGCGGCCTTGTGT‑CGCGAAAg                                                                                                                                               <  1:3896276/150‑1 (MQ=255)

       gAAGCCGGTGGGGCTGGGGATGGCGAGCATTTCCAGCAGCACGCGTTTGAGATAGTCGAGATCGGGTTCGGGG‑AGTCGTTCGGACATGTCCAGGCCTCCGGTCGGG‑TTAAAATTTTGCAGTCTTGTGGGAGCGGCCTTGTGT‑CGCGAAAg                                                                                                                                               <  2:5351639/150‑1 (MQ=255)

       gAAGCCGGTGGGGCTGGGGATGGCGAGCATTTCCAGCAGCACGCGTTTGAGATAGTCGAGATCGGGTTCGGGG‑AGTCGTTCGGACATGTCCAGGCCTCCGGTCGGG‑TTAAAATTTTGCAGTCTTGTGGGAGCGGCCTTGTGT‑CGCGAAAg                                                                                                                                               <  2:4639237/150‑1 (MQ=255)

         aGCCGGTGGGGCTGGGGATGGCGAGCATTTCCAGCAGCACGCGTTTGAGATAGTCGAGATCGGGTTCGGGG‑AGTCGTTCGGACATGTCCAGGCCTCCGGTCGGG‑TTAAAATTTTGCAGTCTTGTGGGAGCGGCCTTGTGT‑CGCGAAAggg                                                                                                                                             <  2:5155668/150‑1 (MQ=255)

          gCCGGTGGGGCTGGGGATGGCGAGCATTTCCAGCAGCACGCGTTTGAGATAGTCGAGATCGGGTTCGGGG‑AGTCGTTCGGACATGTCCAGGCCTCCGGTCGGG‑TTAAAATTTTGCAGTCTTGTGGGAGCGGCCTTGTGT‑CGCGAAAGGGc                                                                                                                                            <  1:197485/150‑1 (MQ=255)

          gCCGGTGGGGCTGGGGATGGCGAGCATTTCCAGCAGCACGCGTTTGAGATAGTCGAGATCGGGTTCGGGG‑AGTCGTTCGGACATGTCCAGGCCTCCGGTCGGG‑TTAAAATTTTGCAGTCTTGTGGGAGCGGCCTTGTGT‑CGCGAAAGGGc                                                                                                                                            <  2:871734/150‑1 (MQ=255)

           ccGGTGGGGCTGGGGAGGGCGAGCATTTCCAGCAGCACGCGTTTGAGATAGTCGAGATCGGGTTCGGGG‑AGTCGTTCGGACATGTCCAGGCCTCCGGTCGGG‑TTAAAATTTTGCAGTCTTGTGGGAGCGGCCTTGTGT‑CGCGAAAGGGcc                                                                                                                                           <  2:2236038/150‑1 (MQ=255)

            cGGTGGGGCTGGGGATGGCGAGCATTTCCAGCAGCACGCGTTTGAGATAGTCGAGATCGGGTTCGGGG‑AGTCGTTCGGACATGTCCAGGCCTCCGGTCGGG‑TTAAAATTTTGCAGTCTTGTGGGAGCGGCCTTGTGT‑CGCGAAAGGGCCg                                                                                                                                          <  2:4283673/150‑1 (MQ=255)

             ggTGGGGCTGGGGATGGCGAGCATTTCCAGCAGCACGCGTTTGAGATAGTCGAGATCGGGTTCGGGG‑AGTCGTTCGGACATGTCCAGGCCTCCGGTCGGG‑TTAAAATTTTGCCGTCTTGTGGGAGCGGCCTTGTGT‑CGCGAAAGGGCCGc                                                                                                                                         <  2:2955601/150‑1 (MQ=255)

             ggTGGGGCTGGGGATGGCGAGCATTTCCAGCAGCACGCGTTTGAGATAGTCGAGATCGGGTTCGGGG‑AGTCGTTCGGACATGTCCAGGCCTCCGGTCGGG‑TTAAAATTTTGCAGTCTTGTGGGAGCGGCCTTGTGT‑CGCGAAAGGGCCGc                                                                                                                                         <  1:3136988/150‑1 (MQ=255)

             ggTGGGGCTGGGGATGGCGAGCATTTCCAGCAGCACGCGTTTGAGATAGTCGAGATCGGGTTCGGGG‑AGTCGTTCGGACATGTCCAGGCCTCCGGTCGGG‑TTAAAATTTTGCAGTCTTGTGGGAGCGGCCTTGTGT‑CGCGAAAGGGCCGc                                                                                                                                         <  2:2955599/150‑1 (MQ=255)

             ggTGGGGCTGGGGATGGCGAGCATTTCCAGCAGCACGCGTTTGAGATAGTCGAGATCGGGTTCGGGG‑AGTCGTTCGGACATGTCCAGGCCTCCGGTCGGG‑TTAAAATTTTGCAGTCTTGTGGGAGCGGCCTTGTGT‑CGCGAAAGGGCCGc                                                                                                                                         <  2:1089314/150‑1 (MQ=255)

              gTGGGGCTGGGGATGGCGAGCATTTCCAGCAGCACGCGTTTGAGATAGTCGAGATCGGGTTCGGGG‑AGTCGTTCGGACATGTCCAGGCCTCCGGTCGGG‑TGAAAATTTTGCAGTCTTGTGGGAGCGGCCCGGTGT‑CGCGCAACGGCCGCa                                                                                                                                        >  2:589416/1‑150 (MQ=255)

                   gCTGGGGATGTCGAGCATTTCCAGCAGCACGCGTTTGAGATAGTCGAGATCGGGTTCGGGG‑AGTCGTTCGGACATGTCCAGGCCTCCGGTCGGG‑TTAAAATTTTGCAGTCTTGTGGGAGCGGCCTTGTGT‑CGCGAAAGGGCCGCAAAGCg                                                                                                                                   <  1:1152337/150‑1 (MQ=255)

                   gCTGGGGATGGCGAGCATTTCCAGCAGCACGCGTTTGAGATAGTCGAGATCGGGTTCGGGG‑AGTCGTTCGGACATGTCCAGGCCTCCGGTCGGG‑TTAAAATTTTGCAGTCTTGTGGGAGCGGCCTTGTGT‑CGCGAAAGGGCCGCAAAGCg                                                                                                                                   <  1:380773/150‑1 (MQ=255)

                   gCTGGGGATGGCGAGCATTTCCAGCAGCACGCGTTTGAGATAGTCGAGATCGGGTTCGGGG‑AGTCGTTCGGACATGTCCAGGCCTCCGGTCGGG‑TTAAAATTTTGCAGTCTTGTGGGAGCGGCCTTGTGT‑CGCGAAAGGGCCGCAAAGCg                                                                                                                                   <  1:92738/150‑1 (MQ=255)

                     tGGGGATGGCGAGCATTTCCAGCAGCACGCGTTTGAGATAGTCGAGATCGGGTTCGGGG‑AGTCGTTCGGACATGTCCAGGCCTCCGGTCGGG‑TTAAAATTTTGCAGTCTTGTGGGAGCGGCCTTGTGT‑CGCGAAAGGGCCGCACAGCGGc                                                                                                                                 >  1:3285715/1‑150 (MQ=255)

                        ggATGGCGAGCATTTCCAGCAGCACGCGTTTGAGATAGTCGAGATCGGGTTCGGGG‑AGTCGTTCGGACATGTCCAGGCCTCCGGTCGGG‑TTAAAATTTTGCAGTCTTGTGGGAGCGGCCTTGTGT‑CGCGAAAGGGCCGCAAAGCGGcccc                                                                                                                              <  1:5557241/150‑1 (MQ=255)

                          aTGGCGAGCTTTTCCAGCAGCACGCGTTTGAGATAGTCGAGATCGGGTTCGGGG‑AGTCGTTCGGACATGTCCAGGCCTCCGGTCGGG‑TTAAAATTTGGCAGTCTTGTGGGAGCGGCCTTGTGT‑CGCGAAAGGGCCGCAAAGCGGCCCC‑gg                                                                                                                           <  1:4843380/150‑1 (MQ=255)

                             gcgcgcTTTTCCCGCAGCACGGTTTTGAGTTCGTCGTGATGGGGTTGGGGG‑TGTGGTTGGGCCTTGTCCCGGCCTCGGGTGGGG‑TTAAATTTTTGCGGTCTTGTGGGAGCGGCCTTGTGT‑GGCTAATGGGCCGCAACGCGGCCCC‑GGCaa                                                                                                                        <  2:2386398/150‑1 (MQ=255)

                               gAGCATTTCCAGCAGCACGCGTTTGAGATAGTCGAGATCGGGTTCGGGG‑AGTCGTTCGGACATGTCCAGGCCTCCGGTCGGT‑TTAAATTTTTGCAGTCTTGTGTGAGCGGCCTTGTGT‑CGCGAAAGGGCCGCAAAGCGGCCCC‑GGCAata                                                                                                                      <  2:413557/150‑1 (MQ=255)

                               gAGCATTTCCAGCAGCACGCGTTTGAGATAGTCGAGATCGGGTTCGGGG‑AGTCGTTCGGACATGTCCAGGCCTCCGGTCGGG‑TTAAAATTTTGCAGTCTTGTGGGAGCGGGCGTGTGG‑GGGGGAAGGGGCGCCAAGCGGCCCC‑GGcccga                                                                                                                      >  1:1834666/1‑146 (MQ=255)

                                       ccAGCAGCACGCGTTTGAGATAGTCGAGATCGGGTTCGGGG‑AGTCGTTCGGACATGTCCAGGCCTCCGGTCGGG‑TTAAAATTTTGCAGTCTTGTTGGAGCGGGCGTGGTT‑GGCGGAATGGGCGCGAAACGGCCCC‑GGCAAAATCGGcccc                                                                                                              >  1:4557536/1‑147 (MQ=255)

                                       ccAGCAGCACGCGTTTGAGATAGTCGAGATCGGGTTCGGGG‑AGTCGTTCGGACATGTCCAGGCCTCCGGTCGGG‑TTAAAATTTTGCAGTCTTGTGGGAGCGGCCGTGTGT‑CGCGCAAGGGCCCCAAAACCGCCCC‑CGCAAAATAGGCAGc                                                                                                              >  2:1300335/1‑150 (MQ=255)

                                          gcagcaCGCGTTTGCGATTGTCGAGATGGGGTTCGGGG‑AGTCGTTCGGACCTGTCCCGGCCTCGGGTCGGG‑TTAAATTTTTGCAGTCTTGTGGGAGCGGCCTTGTGT‑CGCGAACGGGCCGCAACGCGGCCCC‑GGCAATATCGGCAGCGaa                                                                                                           <  2:1857318/150‑1 (MQ=38)

                                          gcagcaCGCGTTTGAGATTGTCGAGATCGGGTTCGGGG‑AGTCGTTCGGACATGTCCAGGCCTCCGGTCGGG‑TTAAAATTTTGCAGTCTTGTGGGAGCGGCCTTGTGT‑CGCGAAAGGGCCGCAAAGCGGCCCC‑GGCAATATCGGCAGCGaa                                                                                                           <  1:5112917/150‑1 (MQ=255)

                                          gcagcaCGCGTTTGAGATAGTCGAGATCGGGTTCGGGG‑AGTCGTTCGGACATGTCCAGGCCTCCGGTCGGG‑TTAAATTTTTGCAGTCTTGTGGGAGCGGCCTTGTGT‑CGCGAAAGGGCCGCAAAGCGGCCCC‑GGCAATATCGGCAGCGaa                                                                                                           <  2:2659379/150‑1 (MQ=255)

                                          gcagcaCGCGTTTGAGATAGTCGAGATCGGGTTCGGGG‑AGTCGTTCGGACATGTCCAGGCCTCCGGTCGGG‑TTAAAATTTTGCAGTCTTGTGGGAGCGGCCTTGTGT‑CGCGAACGGGCCGCAAAGCGGCCCC‑GGCAATATCGGCAGCGaa                                                                                                           <  1:5112902/150‑1 (MQ=255)

                                            agcaCGCGTTTGAGATAGTCGAGATCGGGTTCGGGG‑AGTCGTTCGGACATGTCCAGGCCTCCGGTCGGG‑TTAAAATTTTGCAGTCTTGTGGGAGCGGCCTTGTGT‑CGCGAAAGGGCCGCAAAGCGGCCCC‑GGCAATATCGGCAGCGAAGc                                                                                                         <  2:2264739/150‑1 (MQ=255)

                                             gaaggcttttGAGTTCGTGTAGATGGGGTTGGGGG‑AGTCGTTCGTACATGTCCTGTCCTCGGTTCGGT‑TTAAAATTTGGCCGTGTTGTGGGAGCGGCCTGTTGT‑CGCGAAGGTGCCTCAAAGCGGCCCC‑GTCAATATCGGCAGCGAAGCt                                                                                                        <  2:962458/143‑1 (MQ=255)

                                              caggcgTTTGAGATAGTCGAGAGGGGGTGCGGGG‑AGTCGTTCGGAAATGTCCAGGCCTCCGGTCGGT‑TTAAAATTTTGCTGTGTTGTGGGTGCGGCGTTGTGT‑CGCGAAAGGGCCGAAAAGCGGCCCC‑GGCAATATCGGCAGCGAAGCTg                                                                                                       <  1:1235668/147‑1 (MQ=255)

                                              cacggttttGAGATCGTCTTGAGCGGGTTGGGGG‑AGTCTTTCGGACTTGTCCAGGCCTCGGGTGGGT‑TTAAAATTTGGCCGTCTTGTGGGAGCGGCCTTGGGT‑CGCGAACGGGCCCCAACGCGGCCCC‑GGCAATATCGGCAGCGAAGCTg                                                                                                       <  1:2926243/144‑1 (MQ=255)

                                               aCGCGTTTGAGATAGTCGAGATCGGGTTCGGGG‑AGTCGTTCGGACATGTCCAGGCCTCCGGTCGGG‑TTAAAATTTTGCAGTCTTGTGGGAGCGGCCTTGTGT‑CGCGAAAGGGCCGCAAAGCGGc                                                                                                                                 >  1:5561147/1‑124 (MQ=255)

                                               aCGCGTTTGAGATAGTCGAGATCGGGTTCGGGG‑AGTCGTTCGGACATGTCCAGGCCTCCGGTCGGG‑TTAAAATTTTGCAGTCTTGTGGGAGCGGCCTTGTGT‑CGCGAAAGGGCCGCAAAGCGGc                                                                                                                                 <  2:5561147/124‑1 (MQ=255)

                                               aCGCGTTTGAGATAGTCGAGATCGGGTTCGGGG‑AGTCGTTCGGACATGTCCAGGCCTCCGGTCGGG‑TTAAAATTTTGCAGTCTTGTGGGAGCGGCCATTTTT‑AGCGAAATGGCCCCCAAACCGCCCC‑CGCAAAAAAGGCACCGaaccaga                                                                                                      >  2:663286/1‑145 (MQ=255)

                                               aCGCGTTTGAGATAGTCGAGATCGGGTTCGGGG‑AGTCGTTCGGACATGACCAGGCCTCCGGTCGGG‑TTAAAATTTTGCAGTCTTGTGGGAGCGGGCGTTTGT‑GGGGGAAAGGGCGCCAAGCGGCCCC‑CGCAAAAAAAGCAGCGaaccaga                                                                                                      >  1:1327361/1‑145 (MQ=255)

                                                 gcgTTTGTGATAGTCGAGTTCGGGTTCGGGG‑ATTCTTTCGGACTTGTCCAGGCCTCCGTTCGGG‑TTAAAATTTTGCTGTCTTGTGGTAGCGGCCTTGTTT‑CGCGAAAGGGCCGCAAAGCGGCCCC‑GGCAATATCGGCAGCGAAGCTgaga                                                                                                    <  2:5300478/150‑1 (MQ=38)

                                                 gcgTTTGAGATAGTCGAGATGGGGTTCGGGG‑AGTCTTTCGGACTTGTCCTGGCCTCCGTTCGGG‑TTAAAATTTTGCAGTCTTGTGGGAGCGGCCTTGTGT‑CGCGAAAGGGCCGCAAAGCGGCCCC‑GGCAATATCGGCAGCGAAGCTgaga                                                                                                    <  2:5300151/150‑1 (MQ=38)

                                                     ttGAGATCGTCGAGATCGGGTTCGGGG‑AGTCGTTCGGACATGTCCAGGCCTCCGGTCGGT‑TTAAAATTTTGCCGTCTTGTGGGAGCGGCCTTGTGT‑CGCGAAAGGGCCGCAAAGCGGCCAC‑GGCAATATCGGCAGCGAAGCTGAGACCCt                                                                                                <  2:3214043/150‑1 (MQ=255)

                                                     ttGAGATAGTCGAGATCGGGTTCGGGG‑AGTCGTGCGGACATGTCCCGGCCTCCGGTGGGG‑TTAAAATTTTGCAGTCTTGTGGGAGCGGCCTTGTGT‑CGCGAAAGGGCCGCAAAGCGGCCAC‑GGCAATATCGGCAGCGAAGCTGAGACCCt                                                                                                <  2:3213771/150‑1 (MQ=255)

                                                        agaTAGTCGAGCTGGGGTTCGGGG‑AGTCGTTCGGACATGTCCTGGCCTCCGGTCGGT‑TTAAAATTTTGCAGTCTTGTGGGAGCGGCCTTGGGT‑GGCGAAAGGGCCGCAAAGCGGCCCC‑GGCAATATCGGCAGCGAAGCTGAGACCCTAgg                                                                                             <  2:3775839/150‑1 (MQ=38)

                                                         tatcgTCGTGATGGGGTTGGGGG‑AGTGGGTCGGACTTGTCCAGGCCTCGGGTGGGT‑TTAAAATTTTGCTGTGTTGTGGGAGCGGCGTTGTGT‑GGCGAACGGGCCGCAAAGCGGCCCC‑GGCAATATCGGCAGCGAAGCTGAGACCCTAggg                                                                                            <  1:4401848/146‑1 (MQ=25)

                                                         gaTCGTCGCGAGCGGGTTCGGGG‑AGTCGTTCGGACAGGGCCAGGCCTCCGGTGGGG‑TTAAAATTTTGCCGTCTTGGGGGAGCGGCCTGGTGT‑CGCGAAAGGGCCGCAAAGCGGCCCC‑GGCAATATCGGCAGCGAAGCTGAGACCCTAggg                                                                                            <  1:4157081/150‑1 (MQ=38)

                                                           tAGTCGAGTTGGGGTTCGGGG‑AGTCGTTCGTACTTGTCCTGGCCTCCGGTGGGG‑TTAAATTTTTGCTGTGTTGTGGGAGCGGCTTTGTTT‑CGCGAAAGGGCCGAAATGCGGCCCC‑GGCAATATCGGCAGCGAAGCTGAGACCCTAGGGgc                                                                                          <  1:5481969/150‑1 (MQ=255)

                                                             gTCGTGATGGGGTTGGGGG‑AGTCGTTGGGAATGGTCCCGGCCTCGGGTGGGT‑TTAAAATTTGGCTGTCTTGTGGGTGCGGCCTTGTGT‑CGCGAAAGGGCCGAAAAGCGGCCCC‑GGCAATATCGGCAGCGAAGCTGAGACCCTAGGGgcgc                                                                                        <  1:3962996/150‑1 (MQ=21)

                                                             gTCGAGATCGGGTTCGGGG‑AGTCGTTCGGACATGGCCAGGCCGCCGGTCGGG‑TTAAAATTTGGCCGTCTTGTGGGAGCGGCCTTGTGT‑CGCGAACGGGCCGCAAAGCGGCCCC‑GGCAATATCGGCAGCGAAGCTGAGACCCTAGGGgcgc                                                                                        <  1:4578294/150‑1 (MQ=38)

                                                                  gctggggTGCGGGG‑TGGGGTTCGGACTGGTCCCGGCCTCGGGTGGGT‑TTAAATTTTTGCCGTGGTGGGGGCGCGGCGTTGTGG‑GGCGACCGGGCCGCAACGCGGCCCC‑GGCAATATCGGCAGCGAAGCTGAGACCCTAGGGGCGCTTCGc                                                                                   <  1:3814584/146‑1 (MQ=255)

                                                                  gaTCGGGTTCGGGG‑AGTCGTTCGGACATGTCCAGGCCTCCGGTCGGG‑TTAAAATTTTGCAGTCTTGTGGGAGCGGGCGTGTGT‑GGCGGAAGGGGCGCAAAACCGCCCC‑CCCCAAAAACGCAACCAACCAGAGACCCCAGGGgcgcggagc                                                                                   >  1:5591420/1‑145 (MQ=255)

                                                                   ttggggTTCGGGG‑AGTCGTTCGGACTTGTCCAGGCCTCCGGTGGGG‑TTAAAATTTTGCAGTCTTGTGGGAGCGGCCTTGTGT‑CGCGAAAGGGCCGCAAAGCGGCCCC‑GGCAATATCGGCAGCGAAGCTGAGACCCTAGGGGCGCTTCGCa                                                                                  <  1:3404459/147‑1 (MQ=255)

                                                                   aTCGGGTTCGGGG‑AGTCGTTCGGACATGTCCAGGCCTCCGGTCGGG‑TTAAAATTTTGCAGTCTTGTGGGAGCGGGCGTGTGT‑GGCGAAATGGCCCCAAACCCGCCCC‑CGCAAAAAAGGCCGCCCACCTGAGTCCCAAGGGGGGCGTAGCa                                                                                  >  2:2839098/1‑150 (MQ=255)

                                                                     cGGGTTCGGGG‑AGTCGTTCGGACATGTCCAGGCCTCCGGTCGGG‑TTAAAATTTTGCAGTCTTGTGGGAGCGGCCTTGTGT‑CGCGAAAGGGCCGCAAAGCGGCCCC‑GGCAATATCGGCAGCGAAGCTGAGACCCTAGGGGCGCTTCGCAcc                                                                                <  2:3550791/150‑1 (MQ=255)

                                                                     cGGGTTCGGGG‑AGTCGTTCGGACATGTCCAGGCCTCCGGTCGGG‑TTAAAATTTTGCAGTCTTGTGGGAGCGGCCTTGTGT‑CGCGAAAGGGCCGCAAAGCGGCCCC‑GGCAATATCGGCAGCGAAGCTGAGACCCTAGGGGCGCTTCGCAcc                                                                                <  2:5111685/150‑1 (MQ=255)

                                                                     cGGGTTCGGGG‑AGTCGTTCGGACATGTCCAGGCCTCCGGTCGGG‑TTAAAATTTTGCAGTCTTGTGGGAGCGGCCTTGTGG‑CGCGACAGGGCCGCAAAGCGGCCCC‑GGCAATATCGGCAGCGAAGCTGAGACCCTAGGGGCGCTTCGCAcc                                                                                <  2:3550792/150‑1 (MQ=255)

                                                                      gggTTTGGGGTAGTCGTTGGTAATTGTCCCGGCC‑CCGGTTGGGTTTAAATTTTTGCAGTCTTGTGGGCGCGGCCTTGTGT‑CGCTAACGGGCCCCAAAGCGCCCCC‑GGCAATATCGGCAGCGAAGCTGAGACCCTAGGGGCGCTTCGCAcc                                                                                <  1:540189/150‑1 (MQ=255)

                                                                        gTTCGGGG‑AGTCGTTCGGACATGTCCAGGCCTCCGGTCGGG‑TTAAAATTTTGCAGTCTTGTGGGAGCGGGCCAGTTTGCCGGAATGGGCCCCAAACCGCCCCC‑GCCAA‑ATACGCCCCCCAGCGGTGAACCAACGGGCGCGTCGCAccccc                                                                             >  2:1018588/1‑149 (MQ=255)

                                                                           ggggg‑AGTGGTTCGGACTTGTCCAGGCCTCGGGTGGGT‑TTAAAATTTTGCTGTCTTGTGGGAGCGGCCTTGTGT‑CGCGAAAGGGCCGCAAAGCGGCCCC‑GGCAATATCGGCAGCGAAGCTGAGACCCTAGGGGCGCTTCGCACCCCATTc                                                                          <  1:1256601/149‑1 (MQ=37)

                                                                             ggg‑gGTCGTTCGGACTGGGCCCGGCCTCGGGGCGGG‑TTAAATTGTTGCAGTCTTGTGGGCGGGGCGGGGTGG‑CGCGAACGGGCCGCACCGCGGCCCC‑GGCAATATCGGCAGCGAAGCTGAGACCCTAGGGGCGCTTCGCACCCCATTcgc                                                                        <  1:5331908/150‑1 (MQ=255)

                                                                               g‑AGTCGTTGGGACATGTCCAGGCCTCCGGTCGGG‑TAAAAATTTTGCCGTCTTGTGGGCGCGGCCTTGTGT‑CGCGAACGGGCCGCAAAGCGGCCCC‑GGCAATATCGGCAGCGAAGCTGAGACCCTAGGGGCGCTTCGCACCCCATTCGCGa                                                                      <  2:2034609/150‑1 (MQ=38)

                                                                               g‑AGTCGTTCGGACATGTCCAGGCCTCCGGTCGGG‑TTAAAATTTTGCAGTCTTGTGGGAGCGGCCTTGTGT‑CGCGAGAGGGCCGCAAAGCGGCCCC‑GGCAATATCGGCAGCGAAGCTGAGACCCTAGGGGCGCTTCGCACCCCATTCGCGa                                                                      <  2:2034539/150‑1 (MQ=255)

                                                                                 cgTCGTTCGGACAGGGCCCGGCCTCCGGTCGGG‑TTAAAATTTGGCCGCCTGGTGGGAGCGGCCGGGGGG‑CGCGAACGGGCCGCAACGCGGCCCC‑GGCAATATCGGCAGCGAAGCTGAGACCCTAGGGGCGCTTCGCACCCCATTCGCGac                                                                     <  1:4078524/149‑1 (MQ=255)

                                                                                   tCGTTCGGACATGTCCAGGCCTCCGGTCGGG‑TTAAAATTTTGCAGTCTTGTGGGAGCGGCCGTTTGT‑CGCGAAAGGGCCGCAAAACCGCCCC‑CGCAAAATAGGCAGCCAAACAGAGAACCAACGGGCGCGGAGCACCCCCATAGCGacac                                                                   >  2:1616845/1‑150 (MQ=255)

                                                                                     gTTCGGACATGTCCAGGCCTCCGGTCGGG‑TTAAAATTTTGCAGTCTTGTGGGAGCGGCCATGTGT‑GGCGCAAGGGCCGCAAAGCCGCCCC‑CCCCAAAAAGGCACCGCACCGGAGACCCGAGGGGCGCGGAGCCCCCCCATCGCCACACaa                                                                 >  2:1660214/1‑150 (MQ=255)

                                                                                      ttCGTACTTGTCCTGGCCTCGGGTGGGT‑TTAAAATTTTGCAGTCTTGTGGGCGCGGCCTGGGGT‑CGCGACCGGGCCCCAACGCGGCCCC‑GGCAATATCGGCAGCGAAGCTGAGACCCTAGGGGCGCTTCGCACCCCATTCGCGACACAAg                                                                <  1:1041347/150‑1 (MQ=255)

                                                                                      gtCGGACCGGTCCCGGCCGCCGGTCGGG‑TTAACATTTGGCAGTCTGGTGGGAGCGGCGTTGTGT‑CGCGAACGGGCCGCAACGCGGCCCC‑GGCAATATCGGCAGCGAAGCTGAGACCCTAGGGGCGCTTCGCACCCCATTCGCGACACAAg                                                                <  2:3285715/149‑1 (MQ=255)

                                                                                          gACATGTCCAGGCCTCCGGTCGGG‑TTAAAATATTGCAGTCTTGTGGGAGCGGCCGTTTGT‑GGCGAAAGGGCCGCCAATCGGCCCC‑CGCAAAAAGCGCACCCCAGCAGAGAACCAACGGGCCATGCGCACCCCAATCGCCACACCACGCCg                                                            >  2:1840342/1‑150 (MQ=255)

                                                                                                tCCAGGCCTCCGGTCGGG‑TTAAAATTTTGCAGTCTTGTGGGAGCGGCCTTGTGT‑CGCGAAAGGGCCGCAAAGCGGCCCC‑GGCCATATCGGCAGCGACGCGGAGACCCAAGGGGCGCGTCGCACCCCATTCGCGACACAAGGCCCCTCCca                                                      >  2:3177807/1‑150 (MQ=37)

                                                                                                    ggCCTCCGGTCGGG‑TTAAAATTTTGCAGTCTTGTGGGAGCGGCCATGTGT‑GGCGAAAGGGCCCCAAACCGGCCCC‑CGCCAAAAAGGCAACGGATCGGAGAACCCACGGGCGCGGAGCCCCCCCATCCCGACACAAAGCCCCGccccccaa                                                  >  2:773360/1‑145 (MQ=255)

                                                                                                      ccgcgggTGGGT‑TTAAATTTTTGCCTTCTTGGGGGTGCGGCCTTGGGC‑CGCGAACGGCCCGAAACGCGGCCCC‑GGCAATATCGGCAGCGATGCTGAGACCCTCGGGGCGCTCCGCACCCAATTCGCGACACAAGGCCGCTCCCACaaaaa                                                <  2:1327361/145‑1 (MQ=21)

                                                                                                            ggcGGG‑TTAAATTTTTGCGGTCTTGTGGGAGGGGCCTGGGGT‑CGCGAATGGGCCGCATCGCGGCCCC‑GGCAATATCGGCAGCGAAGCTGAGACCCTTGGGGCGCTACGCACCCCATTCGCGACACAAGGCCGCTCCCACAAAAAGCCGAg                                          <  2:4604319/148‑1 (MQ=25)

                                                                                                             ggggttttAAATTTTTGCTTTCTTGTGGTCGCGGCCTTGTGT‑CGCTAACGGGCCGAAACGCGCCCCC‑GGCAATTTCGGCAGCTACGGTGAGACCCGCGGGGCGCTTCGCACCCCTTGCGCGACACACGGCCGCTCCCACAAAAAGCCGAg                                          <  2:1970421/144‑1 (MQ=255)

                                                                                                                 g‑TTAAAATTTTGCAGTCTTGTGGGAGCGGCCATGTGT‑CGCGAAAGGGCCGCAAAGCGGCCCC‑CGCAAAAACCGCCGCGAAGCAGAGACCCGAGGGGCGCTTCGCAACCCAATCGCCACACAAGGCCCCACCCACAAAAAGCCGAGCAGtt                                     >  2:4616300/1‑150 (MQ=25)

                                                                                                                     aaatTTTTGCTGTGTTGGGGGCGCGGCGTTGTGG‑CGCCACGGGGCCGAACCGCGGCCCG‑GGCAATATGGGCCGCGACGCTGACACCCTCGGGGCGCTCCCCACCCCATCCGCGACACACGGCCGCTCCCACAAAAAGCCGAGCAGTTGTc                                  <  2:1575456/150‑1 (MQ=255)

                                                                                                                      aatttttGCCTTCTTGTGGTCGCGGCCTGTTGT‑GGCGAACGGGCCGAAACGCGGCCCC‑GGCAATATCGGCAGCGAAGCTGAGACCCTAGGGGCGCTTCGCACCCCATTCGCGACACAAGGCCGCTCCCACAAAAAGCCGAGCAGTTGTCa                                 <  1:3672632/147‑1 (MQ=25)

                                                                                                                        aTTGTGCCGTCTGGTGGGAGCGGCCTTGTGT‑CGCGACCGGGCCGCAACGCGGCCCC‑GGCAATATCGGCAGCGACGCTGAGACCCTAGGGGCGCTTCGCACCCCATTCGCGACACAAGGCCGCTCCCACAAAAAGCCGAGCAGTTGTCAgg                               <  2:5602297/150‑1 (MQ=25)

                                                                                                                         ttttGCAGTCTTGTGGGAGCGGCCTTGTGT‑CGCGAAAGGGCCGCAAAGCGGCCCC‑GGCAATATCGGCAGCGAAGCTGAGACCCTAGGGGCGCTTc                                                                                     >  1:342322/1‑95 (MQ=255)

                                                                                                                         ttttGCAGTCTTGTGGGAGCGGCCTTGTGT‑CGCGAAAGGGCCGCAAAGCGGCCCC‑GGCAATATCGGCAGCGAAGCTGAGACCCTAGGGGCGCTTc                                                                                     <  2:342322/95‑1 (MQ=255)

                                                                                                                                 tCTTGTGGGAGCGGCCTTGTTT‑AGCGAAAGGGCCGCAAATCCGCCCC‑CGCAAAATAGGCCGCGAAGCTGAGACCCAAGGGGCGCTTCGCACCCCAATAGCCACACAAGGCCGCTCCCACAAAAAACCCAGCAGTTGTCAGGCCAAcgccc                      >  2:4860675/1‑150 (MQ=21)

                                                                                                                                  gttGGGGGCGCGGCCTGGGGG‑CGCGAACGGGCCCCAACGGGGCCCC‑GGCAATATCGGCCGCGAAGCTGAGACCCTCGGGGCGCTTCGCACCCCAGTCGCGACAAAAGGCCGCTCCCACAAAAAGCCGAGCAGTTGTCAGGCCAAcggccg                     <  2:4402045/149‑1 (MQ=255)

                                                                                                                                  ctggtgGGAGCGGCCTTGGGT‑CGGTACTGGGCCGCAACGGGCCCCG‑GGCAATATGGGCAGCGACGCTGAGACCCTCGGGGCGCTTCGCACCCCTTCCGCGACACAAGGCCGCTCCCACAAAAAGCCGAGCAGTTGTCAGGCCAAcggccg                     <  2:4402067/147‑1 (MQ=25)

                                                                                                                                     gtgGGAGCGGCCTTGTTT‑CGCCAAATGGCCGCAAAGCGGCCCC‑GGCAATATAGGCAGCGCAGCTGATACCCGAGGGGCGGTTAGCACCCCAATCGCGACACAAAGCCCCAAACAAAAAAAGCCGAGAAGTAGTCAGGCCAAcggccggca                  >  2:465609/1‑149 (MQ=21)

                                                                                                                                         ggggggCGTTGTTG‑CGCGACGGGGCCCAAATGGGGCCGC‑GGCTTTTGCGGCAGCTCCGGTGACACCGTGGGGGGGCTTCGAACCCCATCGGCCACAATTGGCCGCTCCCACAAAAAGCCGAGCAGTTGTCAGGCCAACGGCCGGCTATgg              <  2:232858/146‑1 (MQ=255)

                                                                                                                                         gcggggCGTTGTGG‑CGCTAAGGGCCCGCAACGGGCCCCG‑GGCAATTGCGGCAGCGATGCGGCGACCGTTGGGGCGCTTCGCACCCCTTTCGCTACACACGGCCGCTCCCACAAAAAGCCGAGCAGTTGTCAGGCCAACGGCCGGCTATgg              <  2:232788/146‑1 (MQ=255)

                                                                                                                                         gAGCGGCCTTGTGT‑CGCGAAAGTGCCGCAAATCGGCCCC‑GGCAATATCGGCAGCGAAGCTGAGACCCTAGGGGGGCTTCGCACCCCATTCGCGTCACAAGGCCGCTCCCACAAAAAGCCGAGCAGTTGTCAGGCCAACGGCCGGCTAttg              >  2:4244581/1‑148 (MQ=255)

                                                                                                                                            cGGCCTTGTGT‑CGCGAAAGGGCCGCAAAGCGGCCCC‑GGCAATATCGGCAGCGAAGCTGAGACCCTAGGGGCGCTTCGCACCCCATTCGCGACACAAGGCCGCTCCCACAAAAAGCCGAGCAGTTGt                                   <  2:1152960/126‑1 (MQ=255)

                                                                                                                                            cGGCCTTGTGT‑CGCGAAAGGGCCGCAAAGCGGCCCC‑GGCAATATCGGCAGCGAAGCTGAGACCCTAGGGGCGCTTCGCACCCCATTCGCGACACAAGGCCGCTCCCACAAAAAGCCGAGCAGTTGt                                   >  1:1152960/1‑126 (MQ=255)

                                                                                                                                               ccTTGTGT‑CGCGAAAGGGCCGCAAAGCGGCCCC‑GGCAATATCGGCAGCGAAGCTGAGACCCTAGGGGCGCTTCGCACCCCATTCGCGACACAAGGCCGCTCCCACAAAAAGCCGAGCAGTTGTCAGGCCAACGGCCGGTTATGGGGAAAc        >  1:5348108/1‑150 (MQ=255)

                                                                                                                                                 ttGTGT‑CGCGAAAGGGCCGCAAAGCGGCCCC‑GGCAATATCGGCAGCGAAGCTGAGACCCTAGGGGCGCTTCGCACCCCATTCGCGACACAAGGCCGCTCCCACAAAACGCCCAGCAGTTGTCAGGCCAACGGCCGGCTATGGGGAAAcag      >  2:3755175/1‑150 (MQ=255)

                                                                                                                                                 ttGTGT‑CGCGAAAGGGCCGCAAAGCGGCCCC‑GGCAATATCGGCAGCGAAGCTGAGACCCTAGGGGCGCTTCGCACCCCATTCGCGACACAAGGCCGCTCCCACAAAAAGCCGAGCAGTTGTCAGGCCAACGGCCGGCTATGGGGAAAcag      >  1:5386125/1‑150 (MQ=255)

                                                                                                                                                 ttGTGT‑CGCGAAAGGGCCGCAAAGCGGCCCC‑GGCAATATCGGCAGCGAAGCTGAGACCCTAGGGGCGCTTCGCACCCCATTCGCGACACAAGGCCGCTCCCACAAAAAGCCGAGCAGTTGTCAGGCCAACGGCCGGCTATGGGGAAAcag      >  1:5386572/1‑150 (MQ=255)

                                                                                                                                                   gtgt‑CGCTAACGGGCCGAAACGGGGCCCC‑GGAAATTTCGGCAGGGAAGCTTAGACCCTCGGGGCGCTTCGCACCCCATTCGCGACACAAGGCCGCTCCCACAAAAAGCCGAGCAGTTGTCAGGCCAACGGCCGGCTATGGGGAAAcagca    <  2:4851593/150‑1 (MQ=37)

                                                                                                                                                    ggggcgctAACGGGCCCCAACGCGCCCCC‑GGCAATATCGGCAGCGACGCTGACACCCTCGGGGCGCTTCGCACCCCATCCGCCAAACACGGCCGCTCCCACAAAAAGCCGAGCAGTTGTCAGGCCAACGGCCGGCTATGGGGAAAcagca    <  2:4851574/146‑1 (MQ=255)

                                                                                                                                                     gt‑CGCGAAAGGGCCGCAAAGCGGCCCC‑GGCAATATCGGCAGCGAAGCTGAGACCCTAGGGGCGCTTCGCACCCCATTCGCGACACAAGGCCGCTCCCACAAAAAGCCGAGCAGTTGTCAGGCCAACGGCCGGCTATGGGGAAACAGCAgg  >  2:4119391/1‑150 (MQ=255)

                                                                                                                                                     ggccgcgACCGGGCC‑CAAACGGGCCCCGGGCAATAGCGGCCCCGACGCTGAGCCCCTTGGGGCGCTCCGCACCCCTTCCGCGACACACGGCCGCTCCCACAAAAAGCCGAGCAGTTGTCAGGCCAACGGCCGGCTATGGGGAAAcagcag   <  2:544058/147‑1 (MQ=255)

TGTCGGTGAAGCCGGTGGGGCTGGGGATGGCGAGCATTTCCAGCAGCACGCGTTTGAGATAGTCGAGATCGGGTTCGGGG‑AGTCGTTCGGACATGTCCAGGCCTCCGGTCGGG‑TTAAAATTTTGCAGTCTTGTGGGAGCGGCCTTGTGT‑CGCGAAAGGGCCGCAAAGCGGCCCC‑GGCAATATCGGCAGCGAAGCTGAGACCCTAGGGGCGCTTCGCACCCCATTCGCGACACAAGGCCGCTCCCACAAAAAGCCGAGCAGTTGTCAGGCCAACGGCCGGCTATGGGGAAACAGCAGG  >  NC_002947/1952638‑1952934

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