BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,299,500	G→A	5.6%	L863F (CTC→TTC)	sbcC ←	exonuclease subunit SbcC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,299,500	0	G	A	5.6%	557.2 / 14.1	177	L863F (CTC→TTC)	sbcC	exonuclease subunit SbcC
Reads supporting (aligned to +/- strand): ref base G (78/89); new base A (3/7); total (81/96)
Fisher's exact test for biased strand distribution p-value = 3.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 178 total aligned reads displayed.

GGTATTGGCCTTGAGCTCGCTGGCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGT  >  NC_002947/2299351‑2299648

ggTATTGGCCTTGAGCTCGCTGGCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAg                                                                                                                                                      >  2:3240562/1‑150 (MQ=255)

       gCCTTGAGCTCGCTGGCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCg                                                                                                                                               >  1:3292576/1‑150 (MQ=255)

       gCCTTGAGCTCGCTGGCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCg                                                                                                                                               >  2:4141364/1‑150 (MQ=255)

        ccTTGAGCTCGCTGGCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAACTGCTCgg                                                                                                                                              <  1:5205262/150‑1 (MQ=255)

        ccTTGAGCTCGCTGGCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAACTGCTCgg                                                                                                                                              <  1:5206006/150‑1 (MQ=255)

         cTTGAGCTCGCTGGCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGc                                                                                                                                             <  1:1001075/150‑1 (MQ=255)

         cTTGAGCTCGCTGGCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGc                                                                                                                                             <  2:2607970/150‑1 (MQ=255)

         cTTGAGCTCGCTGGCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGc                                                                                                                                             >  1:965689/1‑150 (MQ=255)

         cTTGAGCTCGCTGGCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGc                                                                                                                                             <  2:1224329/150‑1 (MQ=255)

           tGAGCTCGCTGGCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGa                                                                                                                                           <  2:2752798/150‑1 (MQ=255)

           tGAGCTCGCTGGCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGa                                                                                                                                           <  1:411598/150‑1 (MQ=255)

            gAGCTCGCTGGCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAACTGCTCGGCGaa                                                                                                                                          <  2:3212145/150‑2 (MQ=255)

              gCTCGCTGGCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTg                                                                                                                                        >  2:5022380/1‑150 (MQ=255)

                  gCTGGCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGc                                                                                                                                    >  2:4446154/1‑150 (MQ=255)

                     ggCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAg                                                                                                                                 <  2:4149172/150‑1 (MQ=255)

                     ggCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAg                                                                                                                                 <  1:1460122/150‑1 (MQ=255)

                     ggCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAg                                                                                                                                 <  1:4088807/150‑1 (MQ=255)

                     ggCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAg                                                                                                                                 <  1:1660314/150‑1 (MQ=255)

                      gCCAGTTGCACGCCTTGGGTGTGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAgg                                                                                                                                <  2:5424452/150‑1 (MQ=255)

                      gCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAgg                                                                                                                                >  2:2191313/1‑150 (MQ=255)

                       ccAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAggg                                                                                                                               <  2:216498/150‑1 (MQ=255)

                       ccAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAggg                                                                                                                               <  1:1259919/150‑1 (MQ=255)

                        cAGTTGCACGCGTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAgggg                                                                                                                              >  2:155140/1‑149 (MQ=255)

                         aGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCt                                                                                                                             <  2:4515129/150‑1 (MQ=255)

                           ttGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTg                                                                                                                           <  1:2310575/150‑1 (MQ=255)

                            tGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGc                                                                                                                          <  1:12579/150‑1 (MQ=255)

                             gCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCt                                                                                                                         >  2:4191823/1‑150 (MQ=255)

                                  ccTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGt                                                                                                                    <  1:4048978/150‑1 (MQ=255)

                                     tGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGc                                                                                                                 <  2:191188/150‑1 (MQ=255)

                                      gggTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCt                                                                                                                >  1:916767/1‑150 (MQ=255)

                                        gTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCtcg                                                                                                              >  2:1648595/1‑150 (MQ=255)

                                              aGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAgg                                                                                                        <  1:2685790/150‑1 (MQ=255)

                                              aGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAgg                                                                                                        <  1:1026282/150‑1 (MQ=255)

                                              aGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTCGCTCGTCGAgg                                                                                                        <  2:1953590/150‑1 (MQ=255)

                                                 tCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGggcg                                                                                                     >  1:1325286/1‑150 (MQ=255)

                                                  cTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGggcgg                                                                                                    >  2:5444731/1‑150 (MQ=255)

                                                       aGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGt                                                                                               >  2:1736270/1‑150 (MQ=255)

                                                       aGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGt                                                                                               >  2:5364560/1‑150 (MQ=255)

                                                         gcgcTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGttt                                                                                             >  2:839477/1‑150 (MQ=255)

                                                          cgcTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTc                                                                                            <  2:2303212/150‑1 (MQ=255)

                                                           gcTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTct                                                                                           >  1:1709575/1‑150 (MQ=255)

                                                                        ggCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCa                                                                                                 >  1:3165777/1‑131 (MQ=255)

                                                                           gTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAActgctgctg                                                                           >  1:530596/1‑150 (MQ=255)

                                                                            tCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTctg                                                                                          <  1:4695392/134‑1 (MQ=255)

                                                                            tCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTctg                                                                                          >  2:4695392/1‑134 (MQ=255)

                                                                            tCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTgg                                                                          <  2:5181936/150‑1 (MQ=255)

                                                                               cgggTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCCTGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGGTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCt                                                                       >  2:3564396/2‑150 (MQ=255)

                                                                                  gTGCGGGCGTGCTCCAGCGCGGGGTCCATGTGCTGCTGCCAGGCTTCGGAGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGc                                                                    >  1:4453675/1‑150 (MQ=255)

                                                                                    gCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCAc                                                                  <  2:5375957/150‑1 (MQ=255)

                                                                                      gggCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCAcgc                                                                >  2:761182/1‑150 (MQ=255)

                                                                                       ggCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAggg                                                                                                       >  2:2405601/1‑110 (MQ=255)

                                                                                       ggCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAggg                                                                                                       <  1:2405601/110‑1 (MQ=255)

                                                                                        gCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCtt                                                                      <  1:5353415/142‑1 (MQ=255)

                                                                                        gCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCAcgcgc                                                              <  2:361595/150‑1 (MQ=255)

                                                                                        gCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCAcgcgc                                                              <  2:4975359/150‑1 (MQ=255)

                                                                                          gTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGt                                                                                                                                         <  2:14610/73‑1 (MQ=255)

                                                                                          gTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACgcgcgc                                                            <  1:5022380/150‑1 (MQ=255)

                                                                                          gTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACgcgcgc                                                            >  2:2029055/1‑150 (MQ=255)

                                                                                              tCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCctgc                                                        >  2:3828169/1‑150 (MQ=255)

                                                                                              tCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCctgc                                                        >  2:4893305/1‑150 (MQ=255)

                                                                                                cAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCctgctg                                                      >  2:2068214/1‑150 (MQ=255)

                                                                                                 aGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGt                                                     <  1:5071649/150‑1 (MQ=255)

                                                                                                  gcgcGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGTTGGCCTTGCACGCGCGCctgcaaga                                                    <  2:1325286/150‑5 (MQ=255)

                                                                                                  gcgcGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGtt                                                    <  2:22960/150‑1 (MQ=255)

                                                                                                  gcgcGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGtt                                                    <  1:4228702/150‑1 (MQ=255)

                                                                                                   cgcgGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCgg                                                                                                                                              >  2:2274192/1‑59 (MQ=255)

                                                                                                   cgcgGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCgg                                                                                                                                              <  1:2274192/59‑1 (MQ=255)

                                                                                                   cgcgGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTg                                                   <  1:1073516/150‑1 (MQ=255)

                                                                                                    gcgGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTgg                                                  <  2:1116453/150‑1 (MQ=255)

                                                                                                    gcgGTGGCCATGGGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTgg                                                  <  2:2978736/150‑1 (MQ=255)

                                                                                                     cgGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGt                                                 <  1:1972554/150‑1 (MQ=255)

                                                                                                            cATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAg                                          >  2:1038568/1‑150 (MQ=255)

                                                                                                              tgtgCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCt                                        >  1:4338352/1‑150 (MQ=255)

                                                                                                                  ctgctgcCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTc                                    <  1:4994090/150‑1 (MQ=255)

                                                                                                                    gctgcCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCa                                  >  1:884913/1‑150 (MQ=255)

                                                                                                                       gcCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCt                               <  1:928361/150‑1 (MQ=255)

                                                                                                                        ccAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTg                              >  1:5307991/1‑150 (MQ=255)

                                                                                                                          aGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTggg                            >  1:3673073/1‑150 (MQ=255)

                                                                                                                          aGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTggg                            >  2:3989962/1‑150 (MQ=255)

                                                                                                                           ggCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGt                           <  1:482786/150‑1 (MQ=255)

                                                                                                                           ggCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGt                           <  2:2719485/150‑1 (MQ=255)

                                                                                                                              ttCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTg                        >  1:1885490/1‑150 (MQ=255)

                                                                                                                                 ggCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTg                                                   >  1:3213682/1‑120 (MQ=255)

                                                                                                                                  gcgcTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCgg                    >  1:1975893/1‑150 (MQ=255)

                                                                                                                                  gcgcTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCgg                    >  2:2836513/1‑150 (MQ=255)

                                                                                                                                  gcgcTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCgg                    <  1:1648595/150‑1 (MQ=255)

                                                                                                                                  gcgcTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGGTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCgg                    >  2:3332236/1‑150 (MQ=255)

                                                                                                                                      tGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCtt                <  2:884913/150‑1 (MQ=255)

                                                                                                                                        gCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTtgc              <  1:3628450/150‑1 (MQ=255)

                                                                                                                                            gTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTgctgct          <  2:3177416/150‑1 (MQ=255)

                                                                                                                                            gTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTgctgct          >  1:4479351/1‑150 (MQ=255)

                                                                                                                                            gTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTgctgct          >  1:1098756/1‑150 (MQ=255)

                                                                                                                                            actCGCCGAACTGCTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTtgct             <  2:4698753/145‑1 (MQ=255)

                                                                                                                                            actCGCCGAACTGCTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTtgct             >  1:4698753/3‑147 (MQ=255)

                                                                                                                                            actCGCCGAACTGCTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTgctgct          >  1:5391542/3‑150 (MQ=255)

                                                                                                                                                gCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTc      >  2:1126926/1‑150 (MQ=255)

                                                                                                                                                gCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTc      <  2:5316210/150‑1 (MQ=255)

                                                                                                                                                 ccGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTcc     <  2:3292576/150‑1 (MQ=255)

                                                                                                                                                  cGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACACAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCa    <  1:761182/150‑1 (MQ=255)

                                                                                                                                                    aGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCTGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGt  <  1:1762469/150‑1 (MQ=255)

GGTATTGGCCTTGAGCTCGCTGGCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGT  >  NC_002947/2299351‑2299648

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