BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,886,825	C→G	12.7%	A83P (GCA→CCA)	PP_5161 ←	lysophospholipase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,886,825	0	C	G	12.7%	334.6 / 34.0	117	A83P (GCA→CCA)	PP_5161	lysophospholipase
Reads supporting (aligned to +/- strand): ref base C (54/48); new base G (7/8); total (61/56)
Fisher's exact test for biased strand distribution p-value = 7.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 120 total aligned reads displayed.

GCTGGACAGGCCGTGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTGCGCAGGC  >  NC_002947/5886676‑5886971

    gACAGGCCGTGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGGTCGGGCAGGCAtggctgc                                                                                                                                                 >  1:176155/1‑150 (MQ=255)

    gACAGGCCGTGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAtgcctgc                                                                                                                                                 <  1:2262224/150‑1 (MQ=255)

       aGGCCGTGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGC‑TGGCTGCTCGGGCAGGCATGGCTGCCCGa                                                                                                                                             <  2:1828801/150‑1 (MQ=255)

        gggcttGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGGCATGGCTGACCGa                                                                                                                                             <  2:602672/145‑1 (MQ=255)

        ggCCGTGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGGCATGGCTGCCCGa                                                                                                                                             <  2:1270736/150‑1 (MQ=255)

         gCCGTGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGGCATGGCTGCCCGag                                                                                                                                            <  1:1578618/150‑2 (MQ=255)

          ccGTGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGACCCTGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGAcc                                                                                                                                           >  2:192227/1‑150 (MQ=255)

           cGTGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGGCATGGCTGCCCGAGCa                                                                                                                                          <  2:2253694/150‑1 (MQ=255)

           cGTGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCa                                                                                                                                          >  1:790667/1‑150 (MQ=255)

            gTGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAg                                                                                                                                         <  2:1422866/150‑1 (MQ=255)

             tGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAgg                                                                                                                                        >  2:1945108/1‑150 (MQ=255)

                ccGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGGTCGGGCAGGCATGGCTGCCCGAGCAgccaa                                                                                                                                     <  2:176155/150‑5 (MQ=255)

                   ggCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AAc                                                                                                                                 <  2:109751/150‑1 (MQ=255)

                        gTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGc                                                                                                                            >  2:2373659/1‑150 (MQ=255)

                         tCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCt                                                                                                                           >  2:2105238/1‑150 (MQ=255)

                           acaGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGc                                                                                                                         >  2:1029179/1‑150 (MQ=255)

                             aGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAgg                                                                                                                                        >  2:974227/1‑134 (MQ=255)

                                 gATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGaaa                                                                                                                   <  1:2149663/150‑1 (MQ=255)

                                     aCGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTa                                                                                                               >  1:1490662/1‑150 (MQ=255)

                                      cGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AAcccc                                                                                                                              <  2:1470192/134‑1 (MQ=255)

                                       ggCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTAcc                                                                                                             >  2:457680/1‑150 (MQ=255)

                                               ccTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGt                                                                                                                                       >  2:2249658/1‑117 (MQ=255)

                                               ccTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGt                                                                                                                                       <  1:2249658/117‑1 (MQ=255)

                                               ccTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCcagccagc                                                                                                     >  1:1701987/1‑150 (MQ=255)

                                                cTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCt                                                                                                    >  1:2191923/1‑150 (MQ=255)

                                                cTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCt                                                                                                    >  1:2191401/1‑150 (MQ=255)

                                                 tGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCtg                                                                                                   >  1:1111534/1‑150 (MQ=255)

                                                   cTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCtgtg                                                                                                 <  1:1396583/150‑1 (MQ=255)

                                                    ttGCGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCtgtgt                                                                                                <  1:821668/150‑1 (MQ=255)

                                                         tGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCgg                                                                                           >  1:2046599/1‑150 (MQ=255)

                                                         tGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCgg                                                                                           >  1:67139/1‑150 (MQ=255)

                                                         tGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCgg                                                                                           >  1:629364/1‑150 (MQ=255)

                                                         tGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCgg                                                                                           <  2:846862/150‑1 (MQ=255)

                                                           cccACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCCTGCCCCTTTGGGCAGCCAAACCTGCCCG‑CGGTGACGCCATCTGGTTGGCGGAAAGCTACCCAGCCAGCTGTGTGCCGGCa                                                                                         >  1:2378476/1‑150 (MQ=255)

                                                             cACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAgg                                                                                       >  1:1060757/1‑150 (MQ=255)

                                                              aCTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTTGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGt                                                                                      >  1:532620/1‑150 (MQ=255)

                                                              aCTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAg                                                                                        >  2:587446/1‑148 (MQ=255)

                                                              aCTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAg                                                                                        <  1:587446/148‑1 (MQ=255)

                                                                tCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTAACGCTTGGCTAATGAATCAGTCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTc                                                                                    >  1:885440/1‑150 (MQ=255)

                                                                tCGATCACATGTCGGTAGAGGCCAATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGGCATGGCTGCCCGag                                                                                                                                            <  1:994907/95‑2 (MQ=255)

                                                                  gATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGCCCAGGTCG‑AACCCCGCTGCCAGACACCTCCCCCGCAAGCGTGGTGCCGGCAGGGTCAg                                                                                  >  2:515379/1‑150 (MQ=255)

                                                                  gATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAg                                                                                  <  2:349295/150‑1 (MQ=255)

                                                                    tCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCAGTGAAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGc                                                                                <  2:1085245/150‑1 (MQ=255)

                                                                       caTGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCAt                                                                             <  1:2216439/150‑1 (MQ=255)

                                                                        aTGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATa                                                                            >  1:1598970/1‑150 (MQ=255)

                                                                           tCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGaa                                                                         <  2:629364/150‑1 (MQ=255)

                                                                           tCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGaa                                                                         <  1:2373659/150‑1 (MQ=255)

                                                                             ggTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAgc                                                                       >  1:120447/1‑150 (MQ=255)

                                                                                  agGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTgg                                                                  <  1:2360259/150‑1 (MQ=255)

                                                                                  agGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGCCGCGCTTGGCTGCTAGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTgg                                                                  <  1:1594657/150‑1 (MQ=255)

                                                                                    gCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTa                                                                <  2:2191401/150‑1 (MQ=255)

                                                                                    gCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTa                                                                <  2:2191923/150‑1 (MQ=255)

                                                                                      ccATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGt                                                              >  2:208682/1‑150 (MQ=255)

                                                                                           tgGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTgc                                                         <  2:120447/150‑1 (MQ=255)

                                                                                             gTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTa                                                                                                               <  2:1650840/94‑1 (MQ=255)

                                                                                              tCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGTTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCt                                                      >  2:250931/1‑150 (MQ=255)

                                                                                               cGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCtt                                                     <  2:1060757/150‑1 (MQ=255)

                                                                                               cGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCtt                                                     <  1:457680/150‑1 (MQ=255)

                                                                                                gTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTtg                                                    >  2:1435844/1‑150 (MQ=255)

                                                                                                gTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTtg                                                    >  1:394636/1‑150 (MQ=255)

                                                                                                  aaTAGCCGTGCAGCAGGAATTGTGTCGCGCGTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTtgtg                                                  <  1:1945108/150‑1 (MQ=255)

                                                                                                  aaTAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTtgtg                                                  <  2:885440/150‑1 (MQ=255)

                                                                                                         gTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCtgt                                                                                                  <  2:157233/95‑1 (MQ=255)

                                                                                                         gTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCtgt                                                                                                  >  1:157233/1‑95 (MQ=255)

                                                                                                          tGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCGGAAATATACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGt                                          <  2:790667/150‑1 (MQ=255)

                                                                                                           gcagcaGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGtg                                         <  1:1464742/150‑1 (MQ=255)

                                                                                                             agcagGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGtggc                                       >  1:2061100/1‑150 (MQ=255)

                                                                                                                 ggAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGtggctggc                                   <  1:2329944/150‑1 (MQ=255)

                                                                                                                  gAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTtgt                                                   <  1:357467/133‑1 (MQ=255)

                                                                                                                   aaCAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGa                                 >  2:782218/1‑150 (MQ=255)

                                                                                                                     cAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCgcgg                               >  2:1743304/1‑147 (MQ=255)

                                                                                                                      aGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAgtg                              <  2:2061100/150‑1 (MQ=255)

                                                                                                                        cGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAgtggt                            >  2:884961/1‑150 (MQ=255)

                                                                                                                         gTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTc                                                                                    >  2:2439415/1‑93 (MQ=255)

                                                                                                                         gTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAgtggtg                           >  1:571274/1‑150 (MQ=255)

                                                                                                                           cgcgCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCATGGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGaa                         <  1:1106583/150‑1 (MQ=255)

                                                                                                                              gcTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACgg                      >  1:2351723/1‑150 (MQ=255)

                                                                                                                              gcTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACgg                      >  2:744299/1‑150 (MQ=255)

                                                                                                                              gcTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACgg                      >  2:1440145/1‑150 (MQ=255)

                                                                                                                                ttGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTg                                                                   <  1:1694400/103‑1 (MQ=255)

                                                                                                                                ttGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTg                                                                   >  2:1694400/1‑103 (MQ=255)

                                                                                                                                  acctggtCGGGCAGGCATGGCTGCCCGAGCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGcc                  >  2:859035/7‑150 (MQ=255)

                                                                                                                                    ctggtCGGGCAGGCATGGCTGCCCGAGCAGGTCG‑ATCTCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAg                >  1:1378304/1‑150 (MQ=255)

                                                                                                                                    ctggtCGGGCAGGCATGGCTGCCCGAGCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAg                >  1:844413/1‑150 (MQ=255)

                                                                                                                                        tCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCgg            <  1:192227/150‑1 (MQ=255)

                                                                                                                                        tCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCgg            <  1:250931/150‑1 (MQ=255)

                                                                                                                                        tCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCgg            <  1:1435844/150‑1 (MQ=255)

                                                                                                                                         cGGGCAGGCATGGCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCggg           >  2:1828045/1‑150 (MQ=255)

                                                                                                                                          gggCAGGCATGGCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGc              >  2:2199100/1‑146 (MQ=255)

                                                                                                                                          gggCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCgggg          >  1:892679/1‑149 (MQ=255)

                                                                                                                                          gggCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGt          <  1:2271088/150‑1 (MQ=255)

                                                                                                                                           ggCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTg         <  2:1870014/150‑1 (MQ=255)

                                                                                                                                            gCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTgc        >  1:1434934/1‑150 (MQ=255)

                                                                                                                                            gCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTgc        >  2:1875609/1‑150 (MQ=255)

                                                                                                                                               gCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTGCGCa     <  2:67139/150‑1 (MQ=255)

                                                                                                                                                ccATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGTCTGGCGAGTGGTGAACGGCGCCAGGCGGGTGCGCAg    >  1:125220/1‑150 (MQ=255)

                                                                                                                                                ccATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTGCGCAg    >  2:701296/1‑150 (MQ=255)

                                                                                                                                                ccATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTGCGCAg    >  2:2418660/1‑150 (MQ=255)

                                                                                                                                                 cATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTGCGCAgg   <  1:271413/150‑1 (MQ=255)

                                                                                                                                                  aTGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTGCGCAGGc  >  1:2425776/1‑150 (MQ=255)

GCTGGACAGGCCGTGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCG‑AACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTGCGCAGGC  >  NC_002947/5886676‑5886971

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: `‑`