BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	428,857	Δ1 bp	6.4%	coding (153/957 nt)	PP_0351 ←	transmembrane sensor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	428,857	0	G	.	6.4%	355.8 / 17.8	110	coding (153/957 nt)	PP_0351	transmembrane sensor
Reads supporting (aligned to +/- strand): ref base G (45/58); new base . (4/3); total (49/61)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.00e-01

Only 100 of 112 total aligned reads displayed.

CCGACCCATGCCAGCAGCGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCAGGCAATCGCGGCTTCCAGGACCCGGGACGAAACCGGGAGGTTGCT  >  NC_002947/428708‑429006

ccGACCCATGCCAGCAGCGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCg                                                                                                                                                        >  1:501966/1‑150 (MQ=255)

       aTGCCAGCAGCGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑Ggc                                                                                                                                                >  1:2105548/1‑150 (MQ=255)

       aTGCCAGCAGCGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑Ggc                                                                                                                                                <  2:417546/150‑1 (MQ=255)

       aTGCCAGCAGCGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑Ggc                                                                                                                                                >  2:1974178/1‑150 (MQ=255)

          ccAGCAGCGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGcc                                                                                                                                             >  2:1390592/1‑150 (MQ=255)

            agcagcGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAg                                                                                                                                           >  1:1849231/1‑150 (MQ=255)

            agcagcGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAg                                                                                                                                           <  2:2105548/150‑1 (MQ=255)

            agcagcGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAg                                                                                                                                           >  1:1582046/1‑150 (MQ=255)

             gcagcGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAgg                                                                                                                                          >  2:204947/1‑150 (MQ=255)

                 cgcgcCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCAcg                                                                                                                                      <  2:623759/150‑1 (MQ=255)

                  gcgcCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCAcgt                                                                                                                                     <  1:1070832/150‑1 (MQ=255)

                    gcCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCAcgtgc                                                                                                                                   <  1:1594972/150‑1 (MQ=255)

                      cGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGACAcgtgcgt                                                                                                                                 >  1:2211532/1‑150 (MQ=255)

                      cGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCCTGGGTGCATGGCGCCGAGCT‑GGCGCCAGGCAcgtgcgt                                                                                                                                 >  2:1998831/1‑150 (MQ=255)

                      cGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCCCTTACAGCAGCCCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCAcgtgcgt                                                                                                                                 <  2:1399103/150‑1 (MQ=255)

                       gAGCACAAACGTGCCGGCCAGCCCGCCGCTAACCTTGCCAATGCCTCGCCGCAGGCTTGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCAcgtgcgtg                                                                                                                                <  2:910549/150‑1 (MQ=255)

                            gAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAg                                                                                                                                           >  2:1174066/1‑134 (MQ=255)

                            gAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAg                                                                                                                                           <  1:1174066/134‑1 (MQ=255)

                             aaaCGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTc                                                                                                                          <  1:307098/150‑1 (MQ=255)

                             aaaCGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTc                                                                                                                          <  2:1685631/150‑1 (MQ=255)

                             aaaCGTGCCGGCCAGCCCGCCGCAAACCTTCCCAAAGCATTGCCGCAGGTGGGCGCGGGGTTGCAGCAGCGCCCCAGGTAAGATCGGAAGAGC‑‑ACACG‑TATCGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCAGTTGCGGGCTCTTc                                                                                                                          >  1:2351288/1‑150 (MQ=255)

                                 gTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCAAGTCACGTGCGTGCTCTTCGCTg                                                                                                                      >  2:656611/1‑150 (MQ=255)

                                  tGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCTCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGAGTGCTCTTCGCTgg                                                                                                                     <  1:1109329/150‑1 (MQ=255)

                                  tGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGc                                                                                                                                         >  2:2313786/1‑130 (MQ=255)

                                  tGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGc                                                                                                                                         <  1:2313786/130‑1 (MQ=255)

                                      ggCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTggcggc                                                                                                                 <  1:1638044/150‑1 (MQ=255)

                                        ccAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCAt                                                                                                               <  2:863950/150‑1 (MQ=255)

                                            cccGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCa                                                                                                           >  2:2384984/1‑150 (MQ=255)

                                              cgccgcCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAgc                                                                                                         >  2:1761738/1‑150 (MQ=255)

                                                     aaCCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGaa                                                                                                  <  2:1849231/150‑1 (MQ=255)

                                                      aCCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCAT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAAt                                                                                                 <  2:65402/150‑1 (MQ=255)

                                                          tGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGt                                                                                             <  1:204947/150‑1 (MQ=255)

                                                           gCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGtt                                                                                            <  2:1634474/150‑1 (MQ=255)

                                                           gCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGtt                                                                                            <  1:468095/150‑1 (MQ=255)

                                                           gCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGtt                                                                                            <  2:680157/150‑1 (MQ=255)

                                                             cAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGc                                                                                          >  1:2007854/1‑150 (MQ=255)

                                                              aaTGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCg                                                                                         >  1:2154120/1‑150 (MQ=255)

                                                                   gTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCAt                                                                                                     >  2:841333/1‑133 (MQ=255)

                                                                   gTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCAt                                                                                                     <  1:841333/133‑1 (MQ=255)

                                                                   gTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCAATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTtcat                                                                                    <  1:147654/150‑1 (MQ=255)

                                                                     cgccgcAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑A‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCCACCGGCCAGTGCGTGCTCTTCGCTGTCGGCGTGCCAGCTCATGTATTCGTTGCGTtcatca                                                                                  <  2:637212/150‑1 (MQ=255)

                                                                       ccgcAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGa                                                                                                                                                       <  2:1693765/80‑1 (MQ=255)

                                                                       ccgcAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGa                                                                                                                                                       >  1:1693765/1‑80 (MQ=255)

                                                                       ccgcAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCAt                                                                                                     <  1:844274/129‑1 (MQ=255)

                                                                       ccgcAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCAt                                                                                                     >  2:844274/1‑129 (MQ=255)

                                                                        cgcAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCCTCAGGc                                                                               >  1:2213949/1‑150 (MQ=255)

                                                                        cgcAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGc                                                                               <  2:1566081/150‑1 (MQ=255)

                                                                               tGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCgc                                                                        >  1:842583/1‑150 (MQ=255)

                                                                                  cgcgcgATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGAGCATGAATTCGTTGCGTTCATCAGGCGTGCCgctgc                                                                     <  1:1998831/150‑1 (MQ=255)

                                                                                      cgATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCtctc                                                                 <  1:2384984/150‑1 (MQ=255)

                                                                                        aTTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCGTCAGGCGTGCCGCTGCTCTCGc                                                               <  2:698614/150‑1 (MQ=255)

                                                                                             agcagcGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAg                                                           <  2:1185249/149‑1 (MQ=255)

                                                                                             agcagcGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAg                                                           >  1:1185249/1‑149 (MQ=255)

                                                                                              gcagcGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGc                                                         >  1:917167/1‑150 (MQ=255)

                                                                                               cagcGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCt                                                        >  1:1339883/1‑150 (MQ=255)

                                                                                                agcGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTtcat                                                                                    <  1:818070/121‑1 (MQ=255)

                                                                                                   gcCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCaa                                                            <  1:1709220/142‑1 (MQ=255)

                                                                                                   gcCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCaa                                                            >  2:1709220/1‑142 (MQ=255)

                                                                                                      tGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCATGCCGCTGCTCTCGCCAAGGCTGAGTTTc                                                 >  1:1980/1‑150 (MQ=255)

                                                                                                        gCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCa                                               <  2:2149782/150‑1 (MQ=255)

                                                                                                        gCGT‑‑G‑CCGGGCCAGCTGCCCCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCa                                               <  2:1582046/150‑1 (MQ=255)

                                                                                                              g‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCAGGCa                                           >  1:382354/1‑150 (MQ=255)

                                                                                                                     ccAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCAGGCAATcgcg                                     >  2:101911/1‑150 (MQ=255)

                                                                                                                      cAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCAGGCAATCGCgg                                    <  1:101911/150‑1 (MQ=255)

                                                                                                                       aGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCAGGCAATCGCGGc                                   <  1:1761738/150‑1 (MQ=255)

                                                                                                                       aGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCAGGCAATCGCGGc                                   <  2:382354/150‑1 (MQ=255)

                                                                                                                           gCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCAGGCAATCGCGGCTTcc                               <  1:2055780/150‑1 (MQ=255)

                                                                                                                                  tAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCAGGCAATCGCGGCTTCCAGGAccc                        >  2:1922666/1‑150 (MQ=255)

                                                                                                                                   aaCGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCAGGCAATCGCGGCTTCCAGGACCCg                       <  2:831853/150‑1 (MQ=255)

                                                                                                                                   aaCGCGGTGTTCCATTGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCAGGCAATCGCGGCTTCCAGGACCCg                       <  2:1339883/150‑1 (MQ=255)

                                                                                                                                    aCGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCAGGCAATCGCGGCTTCCAGGACCCgg                      <  2:842583/150‑1 (MQ=255)

                                                                                                                                      gcgTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGctc                                                                   <  1:2097312/103‑1 (MQ=255)

                                                                                                                                      gcgTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGctc                                                                   >  2:2097312/1‑103 (MQ=255)

                                                                                                                                      gcgTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGACATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCAGGCAATCGCGGCTTCCAGGACCCGGGa                    >  1:575730/1‑150 (MQ=255)

                                                                                                                                       cgTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCAGGCAATCGCGGCTTCCAGGACCCGGGAc                   <  1:51638/150‑1 (MQ=255)

                                                                                                                                        gTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCAGGCAATCGCGGCTTCCa                              <  2:926088/138‑1 (MQ=255)

                                                                                                                                        gTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCAGGCAATCGCGGCTTCCa                              >  1:926088/1‑138 (MQ=255)

                                                                                                                                        gTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCAGGCAATCGCGGCTTCCAGGACCCGGGACg                  >  1:2042695/1‑150 (MQ=255)

                                                                                                                                           ggTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCAGGCAATCGCGGCTTCCAGGACCCGGGACGaaa               <  2:2042695/150‑1 (MQ=255)

                                                                                                                                            gTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCAGGCAATCGCGGCTTCCAGGACCCGGGACGAAAc              >  1:800873/1‑150 (MQ=255)

                                                                                                                                             tCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGGTTCCAGGCAATCGCGGCTTCCAGGACCCGGGACGAAAcc             >  2:315117/1‑150 (MQ=255)

                                                                                                                                              gcctGGCGCC‑AGCTCGGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCAGGCAATCGCGGCTTCCAGGACCCGGGACGAAAccc            >  2:1384807/4‑149 (MQ=255)

                                                                                                                                              ccATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCAGGCAATCGCGGCTTCCAGGACCCGGGACGAAACCg            <  2:807340/150‑1 (MQ=255)

                                                                                                                                              ccATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCAGGCAATCGCGGCTTCCAGGACCCGGGACGAAACCg            >  2:346094/1‑150 (MQ=255)

                                                                                                                                               cATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCAGGCAATCGCGGCTTCCAGGACCCGGGACGAAACCgg           <  2:496458/150‑1 (MQ=255)

                                                                                                                                               cATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCAGGCAATCGCGGCTTCCAGGACCCGGGACGAAACCgg           >  1:2326702/1‑150 (MQ=255)

                                                                                                                                                ctGGCGCC‑AGCTCGGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCaa                                                            >  2:968882/2‑100 (MQ=255)

                                                                                                                                                ctGGCGCC‑AGCTCGGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCaa                                                            <  1:968882/99‑1 (MQ=255)

                                                                                                                                                ctGGCGCC‑AGCTCGGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCa                                               <  2:376055/112‑1 (MQ=255)

                                                                                                                                                ctGGCGCC‑AGCTCGGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCa                                               >  1:376055/2‑113 (MQ=255)

                                                                                                                                                ctGGCGCC‑AGCTCGGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCa                                               >  2:1111122/2‑113 (MQ=255)

                                                                                                                                                ctGGCGCC‑AGCTCGGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCa                                               <  1:1111122/112‑1 (MQ=255)

                                                                                                                                                 tGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCAGGCAATCGCGGCTTCCAGGACCCGGGACGAAACCggg          >  2:1211724/1‑149 (MQ=255)

                                                                                                                                                 tGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCAGGCAATCGCGGCTTCCAGGACCCGGGACGAAACCggg          <  1:1211724/149‑1 (MQ=255)

                                                                                                                                                 tGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCAGGCAATCGCGGCTTCCAGGACCCGGGACGAAACCGGGa         <  1:315117/150‑1 (MQ=255)

                                                                                                                                                    cgcCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCAGGCAATCGCGGCTTCCAGGACCCGGGACGAAACCGGGAGGt      >  1:784360/1‑150 (MQ=255)

                                                                                                                                                     gcCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGTCGCTGCTCTCGCCAAGGCTGAGTTTCCAGGCAATCGCGGCTTCCAGGACCCGGGa                    <  1:1885836/135‑1 (MQ=255)

                                                                                                                                                     gcCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCAGGCAATCGCGGCTTCCAGGACCCGGGACGAAACCGGGAGGtt     <  2:501966/150‑1 (MQ=255)

                                                                                                                                                     gcCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGACAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGTTGAGTTTCCAGGCAATCGCGGCTTCCAGGACCCGGGACGAAACCGGGAGGtt     >  1:826290/1‑150 (MQ=255)

                                                                                                                                                        gAGCT‑GGCGCCAGGCACGTTCGTGCTCGTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCAGGCAATCGCGGCTTCCAGGACCCGGGACGAAACCGGGAGGTTGCt  <  2:749565/150‑1 (MQ=255)

CCGACCCATGCCAGCAGCGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCAGGCAATCGCGGCTTCCAGGACCCGGGACGAAACCGGGAGGTTGCT  >  NC_002947/428708‑429006

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