BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	428,861:1	+C	6.5%	coding (149/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,861	1	.	C	6.5%	438.4 / 22.8	108	coding (149/957 nt)	PP_0351	transmembrane sensor
Reads supporting (aligned to +/- strand): ref base . (44/57); new base C (4/3); total (48/60)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.29e-01

Only 100 of 108 total aligned reads displayed.

ATGCCAGCAGCGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCAGGCAATCGCGGCTTCCAGGACCCGGGACGAAACCGGGAGGTTGCT  >  NC_002947/428715‑429006

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:1582046/1‑150 (MQ=255)

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

     agcagcGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAg                                                                                                                                           <  2:2105548/150‑1 (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)

                gAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCAcgtgcgtg                                                                                                                                <  1:1917069/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                                                                                                                          <  2:1685631/150‑1 (MQ=255)

                      aaaCGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTc                                                                                                                          <  1:307098/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                                                                                                                                         <  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)

                                 ccAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGACACGTGCGTGCTCTTCGCTGGCGGCAt                                                                                                               <  2:2211532/150‑1 (MQ=255)

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

                                     cccGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCa                                                                                                           <  1:342039/150‑1 (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                                                                                            <  2:542122/150‑1 (MQ=255)

                                                    gCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGtt                                                                                            <  1:468095/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                                                                                                     <  1:841333/133‑1 (MQ=255)

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

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

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

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

                                                                ccgcAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCAt                                                                                                     <  1:844274/129‑1 (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)

                                                                                         agcGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTtcat                                                                                    >  2:818070/1‑121 (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:1509785/150‑1 (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)

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

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

                                                                                                                                 gTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCAGGCAATCGCGGCTTCCAGGACCCGGGACg                  >  1:807340/1‑150 (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‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGTCGCTGCTCTCGCCAAGGCTGAGTTTCCAGGCAATCGCGGCTTCCAGGACCCGGGa                    >  2:1885836/1‑135 (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)

ATGCCAGCAGCGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGT‑‑G‑CCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCAGGCAATCGCGGCTTCCAGGACCCGGGACGAAACCGGGAGGTTGCT  >  NC_002947/428715‑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: `‑`