BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,289,835	A→T	5.7%	L180Q (CTG→CAG)	PP_4663 ←	transmembrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,289,835	0	A	T	5.7%	398.0 / 9.4	122	L180Q (CTG→CAG)	PP_4663	transmembrane protein
Reads supporting (aligned to +/- strand): ref base A (54/61); new base T (4/3); total (58/64)
Fisher's exact test for biased strand distribution p-value = 7.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 122 total aligned reads displayed.

CGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCAAGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGCTGGCATCACGAACACCCCGGTCGACGCAGTGACCACACCCGTGACCAGCCCGCAAAGCGGTCCCAGCCAGCGCTCGTGTCCCGGGGC  >  NC_002947/5289688‑5289977

cGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCCAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCAAGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAgg                                                                                                                                              <  1:738521/150‑1 (MQ=255)

cGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCAAGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAgg                                                                                                                                              <  2:84874/150‑1 (MQ=255)

    gACGCAGCCACGGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCAAGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCt                                                                                                                                          <  1:599819/150‑1 (MQ=255)

       gCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCAAGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGa                                                                                                                                       >  2:442071/1‑150 (MQ=255)

             aCTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCAAGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCt                                                                                                                                 >  2:1344832/1‑150 (MQ=255)

             aCTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCAAGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCt                                                                                                                                 >  1:1343347/1‑150 (MQ=255)

                   ccGCCAGCATGCCGAGCAAGGCCGGCGCCAGAATCAGCAGAGACGCACCCAGCGCCTGGCCACCCAAGGCATCCTGCCCGGCCAGGCCTAGGGACAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTcgc                                                                                                                           <  2:1427630/150‑1 (MQ=255)

                    cGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCAAGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTcgcg                                                                                                                          <  2:1249520/150‑1 (MQ=255)

                     gCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCAAGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCgg                                                                                                                         <  2:968953/150‑1 (MQ=255)

                       cAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCAAGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCGCGCGGCt                                                                                                                       >  1:36404/1‑150 (MQ=255)

                        aGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCAAGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTc                                                                                                                      >  1:2152100/1‑150 (MQ=255)

                        aGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCAAGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTc                                                                                                                      >  1:979287/1‑150 (MQ=255)

                               cGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCAAGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCa                                                                                                               <  2:1957445/150‑1 (MQ=255)

                               cGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCAAGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCa                                                                                                               >  2:1095787/1‑150 (MQ=255)

                                   cAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCAAGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGct                                                                                                           >  1:1333076/1‑150 (MQ=255)

                                       gCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCAAGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGc                                                                                                       <  2:812911/150‑1 (MQ=255)

                                       gCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCAAGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGc                                                                                                       <  2:1652195/150‑1 (MQ=255)

                                       gCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCAAGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGc                                                                                                       <  2:1333076/150‑1 (MQ=255)

                                       gCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCAAGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGc                                                                                                       <  1:1344832/150‑1 (MQ=255)

                                        ccGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCAAGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCTGGGc                                                                                                                                            <  2:1063849/112‑1 (MQ=255)

                                        ccGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCAAGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCTGGGc                                                                                                                                            >  1:1063849/1‑112 (MQ=255)

                                           gcgcCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCAAGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGt                                                                                                   <  1:2153017/150‑1 (MQ=255)

                                             gcCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCAAGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTaa                                                                                                 <  2:1214090/150‑1 (MQ=255)

                                               cAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCAAGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAgg                                                                                               >  2:1141772/1‑150 (MQ=255)

                                                 gCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCAAGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCa                                                                                             <  2:979287/150‑1 (MQ=255)

                                                  cATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCAAGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCAc                                                                                            >  1:109605/1‑150 (MQ=255)

                                                  cATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCAAGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCAc                                                                                            <  2:473951/150‑1 (MQ=255)

                                                     cagcagCGACGCCCCCAGCGCCTGGCCACCCAAGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGc                                                                                         >  2:369726/1‑150 (MQ=255)

                                                       gcagcGACGCCCCCAGCGCCTGGCCACCCAAGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGCTg                                                                                       <  2:696647/150‑1 (MQ=255)

                                                              cGCCCCCAGCGCCTGGCCACCCAAGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGCTGGCATCAc                                                                                <  2:1783385/150‑1 (MQ=255)

                                                                 ccccAGCGCCTGGCCACCCAAGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGCTGGCATCACGaa                                                                             <  1:2024470/150‑1 (MQ=255)

                                                                 ccccAGCGCCTGGCCACCCAAGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGCTGGCATCACGaa                                                                             <  1:820041/150‑1 (MQ=255)

                                                                 ccccAGCGCCGGGTCACCCAAGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGCTGGCATCACGaa                                                                             <  2:2243501/150‑1 (MQ=255)

                                                                  cccAGCGCCTGGCCACCCAAGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCGGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGCTGGCATCACGAac                                                                            >  2:542041/1‑150 (MQ=255)

                                                                   ccAGCGCCTGGCCACCCAAGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGCTGGCATCACGAaca                                                                           >  1:1628480/1‑150 (MQ=255)

                                                                          cTGGCCACCCAAGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGCTGGCATCACGAACACCCCGGt                                                                    >  2:954294/1‑150 (MQ=255)

                                                                          cTGGCCACCCAAGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCACGCTCTGCAGGTAAGGCACCGCTGGCATCACGAACACCCCGGt                                                                    >  1:112776/1‑150 (MQ=255)

                                                                           tGGCCACCCAAGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGCTGGCATCACGAACACCCCGGTc                                                                   <  2:2176159/150‑1 (MQ=255)

                                                                            ggCCACCCAAGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAg                                                                                                >  1:2322209/1‑120 (MQ=255)

                                                                            ggCCACCCAAGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAg                                                                                                <  2:2322209/120‑1 (MQ=255)

                                                                            ggCCACCCAAGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGCTGGCATCACGAACACCCCGGTca                                                                  >  1:1220586/1‑149 (MQ=255)

                                                                               cACCCAAGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCTGGGCCTgtc                                                                                                                                       >  1:2318138/1‑76 (MQ=255)

                                                                               cACCCAAGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCTGGGCCTgtc                                                                                                                                       <  2:2318138/78‑3 (MQ=255)

                                                                               cACCCAAGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAg                                                                                                     <  2:198789/112‑1 (MQ=255)

                                                                               cACCCAAGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAg                                                                                                     >  1:198789/1‑112 (MQ=255)

                                                                               cACCCAAGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGCTGGCATCACGAACACCCCGGTCGACg                                                               <  1:78582/150‑1 (MQ=255)

                                                                                aCCCAAGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGCTGGCATCACGAACACCCCGGTCGACGc                                                              >  2:2117383/1‑150 (MQ=255)

                                                                                 cccAAGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTaa                                                                                                 >  1:693998/1‑114 (MQ=255)

                                                                                   cAAGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGCTGGCATCACGAACACCCCGGTCGACGCAGt                                                           >  1:1470516/1‑150 (MQ=255)

                                                                                    aaGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGCTGGCAt                                                                                   <  2:2115599/125‑1 (MQ=255)

                                                                                    aaGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGCTGGCAt                                                                                   >  1:2115599/1‑125 (MQ=255)

                                                                                    aaGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGCTGGCATCACGAACACCCCGGTCGACGCAGTg                                                          >  1:1386314/1‑150 (MQ=255)

                                                                                     aGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCa                                                                                                      <  1:1429604/105‑1 (MQ=255)

                                                                                     aGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCa                                                                                                      >  2:1429604/1‑105 (MQ=255)

                                                                                      ggCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGCTGGCATCACGAACACCCCGGTCGACGCAGTGAc                                                        <  2:1626721/150‑1 (MQ=255)

                                                                                      ggCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGCTGGCATCACGAACACCCCGGTCGACGCAGTGAc                                                        <  1:369726/150‑1 (MQ=255)

                                                                                      ggCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGCTGGCATCACGAACACCCCGGTCGACGCAGTGAc                                                        <  2:2248379/150‑1 (MQ=255)

                                                                                            cTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCATGCTCTGCAGGTAAGGCACCGCTGGCATCACGAACACCCCGGTCGACGCAGTGAccggggt                                                  <  2:112776/150‑6 (MQ=255)

                                                                                               cccGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGCTGGCATCACGAACACCCCGGTCGACGCAGTGACCACACCCGt                                               <  1:1095787/150‑1 (MQ=255)

                                                                                                ccGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGCTGGCATCACGAACAccc                                                                        <  2:642577/124‑1 (MQ=255)

                                                                                                ccGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGCTGGCATCACGAACAccc                                                                        >  1:642577/1‑124 (MQ=255)

                                                                                                 cGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGGACCATTGCCACCTCGAACACCCCGGTCGACGCAGTGACCACACCCGTGa                                             >  2:635420/1‑150 (MQ=255)

                                                                                                  ggCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGCTGGCATCACGAACACCCCGGTCGACGCAGTGACCACACCCGTGAc                                            <  1:1475625/150‑1 (MQ=255)

                                                                                                  ggCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGCTGGCATCACGAACACCCCGGTCGACGCAGTGACCACACCCGTGAc                                            <  2:503684/150‑1 (MQ=255)

                                                                                                  ggCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACAACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGCTGGCATCACGAACACCCCGGTCGACGCAGTGACCACACCCGTGAc                                            <  1:635420/150‑1 (MQ=255)

                                                                                                          cTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGCTGGCATCACGAACACCCCGGTCGACGCAGTGACCACACCCGTGACCAGCCCGc                                    >  2:1986337/1‑150 (MQ=255)

                                                                                                          cTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGCTGGCATCACGAACACCCCGGTCGACGCAGTGACCACACCCGTGACCAGCCCGc                                    >  1:366367/1‑150 (MQ=255)

                                                                                                               gCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGATAACGCAGAGCTGGCATCACGAACACCCCGGTCGACGCAGTGACCACACCCGTGACCAGCCCGCAAAGc                               >  1:1488904/1‑150 (MQ=255)

                                                                                                                ccAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGCTGGCATCACGAACACCCCGGTCGACGCAGTGACCACACCCGTGa                                             <  1:426134/135‑1 (MQ=255)

                                                                                                                ccAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGCTGGCATCACGAACACCCCGGTCGACGCAGTGACCACACCCGTGACCAGCCCGCAAAGCg                              <  1:2117383/150‑1 (MQ=255)

                                                                                                                ccAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGCTGGCATCACGAACACCCCGGTCGACGCAGTGACCACACCCGTGACCAGCCCGCAAAGCg                              >  1:2313035/1‑150 (MQ=255)

                                                                                                                ccAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGGAAGGCACTGAAGGCAGCACGAACACCCCGGTCGACGCAGTCACCACACCCGTGACCAGCCCGCAAAGCg                              >  2:1531904/1‑150 (MQ=255)

                                                                                                                 cAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGCTGGCATCACGAACACCCCGGTCGACGCAGTGACCACACCCGTGACCAGCCCGCAAAGCgg                             >  2:648155/1‑150 (MQ=255)

                                                                                                                       cAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGCTGGCATCACGaa                                                                             <  2:90205/96‑1 (MQ=255)

                                                                                                                        aGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGCTGGCATCACGAACACCCCGGTCGACGCAGTGACCACACCCGTGACCAGCCCGCAAAGCGGTCccagc                      <  2:1470516/150‑1 (MQ=255)

                                                                                                                        aGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGCTGGCATCACGAACACCCCGGTCGACGCAGTGACCACACCCGTGACCAGCCCGCAAAGCGGTCccagc                      <  1:954294/150‑1 (MQ=255)

                                                                                                                        aGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGCTGGCATCACGAACACCCCGGTCGACGCAGTGACCACACCCGTGACCAGCCCGCAAAGCGGTCccagc                      <  2:718296/150‑1 (MQ=255)

                                                                                                                         gggTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGCTTGCATCACGAACACCCCGGTCGACGCAGTGACCACACCCGTGACCAGCCCGCAAAGCGGTCccagcc                     >  2:2373615/1‑150 (MQ=255)

                                                                                                                         gggTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGCTGGCATCACGAACACCCCGGTCGACGCAGTGACCACACCCGTGACCAGCCCGCAAAGCGGTCccagcc                     >  2:604409/1‑150 (MQ=255)

                                                                                                                           gTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGCTGGCATCACGAACACCCCGGTCGACGCAGTGACCACACCCGTGACCAGCCCGCaaa                                 <  2:256998/136‑1 (MQ=255)

                                                                                                                           gTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGCTGGCATCACGAACACCCCGGTCGACGCAGTGACCACACCCGTGACCAGCCCGCAAAGCGGTCccagccag                   >  1:1949533/1‑150 (MQ=255)

                                                                                                                                  ccGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGCTGGCATCAc                                                                                >  1:238434/1‑82 (MQ=255)

                                                                                                                                   cGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAg                                                                                                     <  1:859259/60‑1 (MQ=255)

                                                                                                                                   cGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAg                                                                                                     >  2:859259/1‑60 (MQ=255)

                                                                                                                                   cGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGCTGGCATCACGAACACCCCGGTCGACGCAGTGACCACACCCGTGACCAGCCCGCAAAGCGGTCCCAGCCAGCGCTCgtg           >  2:2052789/1‑150 (MQ=255)

                                                                                                                                   cGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGCCACCGCTGGCATCACGAACACCCCGGTCGACGCAGTGACCACACCCGTGACCAGCCCGCAAAGCGGTCCCAGCCAGCGCTCgtg           >  1:870380/1‑150 (MQ=255)

                                                                                                                                    gTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGCTGGCATCACGAACACCCCGGTCGACGCAGTGACCACACCCGTGACCAGCCCGCAAAGCGGTCCCAGCCAGCGCTCgtgt          >  1:1597254/1‑150 (MQ=255)

                                                                                                                                    gTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGCTGGCATCACGAACACCCCGGTCGACGCAGTGACCACACCCGTGACCAGCCAAGAATGTGGAGCCAGCCAGCGCTCgtgt          >  2:560986/1‑150 (MQ=255)

                                                                                                                                     tAAACGACAGGCCCAGGGCCTGGATCATCTGCTcgcg                                                                                                                          >  2:1865570/1‑37 (MQ=255)

                                                                                                                                     tAAACGACAGGCCCAGGGCCTGGATCATCTGCTcgcg                                                                                                                          <  1:1865570/37‑1 (MQ=255)

                                                                                                                                     tAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGCTGGCATCACGAACACCCCGGTCGACGCAGTGACCACACCCGTGACCAGCCCGCAAAGCGGTCCCAGCCAGCGCTc              >  1:1593036/1‑145 (MQ=255)

                                                                                                                                     tAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGCTGGCATCACGAACACCCCGGTCGACGCAGTGACCACACCCGTGACCAGCCCGCAAAGCGGTCCCAGCCAGCGCTc              <  2:1593036/145‑1 (MQ=255)

                                                                                                                                        gatccAGGCCCTGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGCTGGCATCACGAACACCCCGGTCGACGCAGTGa                                                         <  1:1342377/95‑1 (MQ=255)

                                                                                                                                        aCGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGCTGGCATCACGAACACCCCGGTCGACGCAGTGACCACACCCGTGACCAGCCCGCAAAGCGGTCCCAGCCAGCGCTCGTGTCCCg      <  2:393321/150‑1 (MQ=255)

                                                                                                                                          gACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGCTGGCATCACGAACACCCCGGTCGACGCAGTGACCACACCCGTGACCAGCCCGCAAAGCGGTCCCAGCCAGCGCTCGTGTCCCggg    <  1:1589804/150‑1 (MQ=255)

                                                                                                                                          gACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGCTGGCATCACGAACACCCCGGTCGACGCAGTGACCACACCCGTGACCAGCCCGCAAAGCGGTCCCAGCCAGCGCTCGTGTCCCggg    >  1:2177079/1‑150 (MQ=255)

                                                                                                                                           aCAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGCTGGCATCACGAACACCCCGGTCGACGCAGTGACCACACCCGTGACCAGCCCGCAAAGCGGTCCCAGCCAGCGCTCGTGTCCCgggg   >  2:1991379/1‑150 (MQ=255)

                                                                                                                                            cAGGCCCTGGGCCTGGATCATCTGCTCGCGGCTAAAGCCCAGGCTCTGCAGGTAAGGCACCGCTGGCATCACGAACACCCCGGTCGACGCAGTGACCACACCCGTGACCAGCCCGCAAAGCGGTCCCAGCCAGCGCTCGCATCCCGGGGc  >  1:1575457/1‑150 (MQ=255)

                                                                                                                                             aGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGCTGGCATCACGAACACCCCGGTCGACGCAGTGACCACACCCGTGACCAGCCCGCaaa                                 >  2:426115/1‑118 (MQ=255)

                                                                                                                                             aGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGCTGGCATCACGAACACCCCGGTCGACGCAGTGACCACACCCGTGACCAGCCCGCaaa                                 <  1:426115/118‑1 (MQ=255)

CGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCAAGGCATCCTGCCCGGCCAGGCCTAGGGCCAAGGCCAGGGTCGATACCGTAAACGACAGGCCCAGGGCCTGGATCATCTGCTCGCGGCTCAAGCCCAGGCTCTGCAGGTAAGGCACCGCTGGCATCACGAACACCCCGGTCGACGCAGTGACCACACCCGTGACCAGCCCGCAAAGCGGTCCCAGCCAGCGCTCGTGTCCCGGGGC  >  NC_002947/5289688‑5289977

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