BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,729,085:1	+T	7.2%	coding (333/645 nt)	PP_3292 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,729,085	1	.	T	7.2%	630.9 / 42.5	153	coding (333/645 nt)	PP_3292	membrane protein
Reads supporting (aligned to +/- strand): ref base . (59/83); new base T (2/9); total (61/92)
Fisher's exact test for biased strand distribution p-value = 2.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 157 total aligned reads displayed.

GCGCTAGGCTGGGCGGTGCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGCACGTGCCT  >  NC_002947/3728936‑3729234

  gcTAGGCTGGGCGGTGCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAATGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑Gc                                                                                                                                                     <  2:2219876/150‑1 (MQ=255)

   cTAGGCTGGGCGGTGCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑Gcc                                                                                                                                                    <  2:2114956/150‑1 (MQ=255)

   cTAGGCTGGGCGGTGCAGCGGGGCGTGAAGCGCACCGTCGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGTCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATGCTGCTGGTGGTGCTGTTGTGCCTG‑Gcc                                                                                                                                                    <  1:5040385/150‑1 (MQ=255)

    tAGGCTGGGCGGTGCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑Gccc                                                                                                                                                   <  2:5390619/150‑1 (MQ=255)

    tAGGCTGGGCGGTGCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑Gccc                                                                                                                                                   >  2:1210562/1‑150 (MQ=255)

       gCTGGGCGGTGCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTg‑g                                                                                                                                                      <  1:3637312/144‑1 (MQ=255)

       gCTGGGCGGTGCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCctgc                                                                                                                                                <  2:2122881/150‑1 (MQ=255)

         tGGGCGGTGCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGACTG‑GCCctgctg                                                                                                                                              >  2:330467/1‑150 (MQ=255)

          gggCGGTGCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCTCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTgg                                                                                                                                             <  1:3820187/150‑1 (MQ=255)

          gggCGGTGCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTgg                                                                                                                                             <  2:4858505/150‑1 (MQ=255)

          gggCGGTGCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTgg                                                                                                                                             <  2:4858559/150‑1 (MQ=255)

           ggCGGTGCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTggg                                                                                                                                            <  2:414650/150‑1 (MQ=255)

             cGGTGCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGca                                                                                                                                          <  2:107740/150‑1 (MQ=255)

               gTGCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGTCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGcaca                                                                                                                                        <  2:1009558/150‑1 (MQ=255)

               gTGCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGcaca                                                                                                                                        <  1:2131948/150‑1 (MQ=255)

               gTGCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGcaca                                                                                                                                        >  2:1676624/1‑150 (MQ=255)

                tGCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGTGCCCAGCAGGGccag                                                                                                                                        <  1:5318035/150‑4 (MQ=255)

                tGCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCAGCAGGGccagg                                                                                                                                       <  2:246254/150‑5 (MQ=255)

                 gccgCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAgg                                                                                                                                      >  2:342808/4‑150 (MQ=255)

                 gCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGTGCCCAGCAGGGccagg                                                                                                                                       <  2:3507726/150‑5 (MQ=255)

                 gCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGTGCCCAGCAGGGccagg                                                                                                                                       <  1:2377092/150‑5 (MQ=255)

                  cAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGc                                                                                                                                     >  1:5201546/1‑150 (MQ=255)

                   aGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCa                                                                                                                                    >  1:2428453/1‑150 (MQ=255)

                    gCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGTGCCCAGCAGGGC‑CAGGCAc                                                                                                                                   <  2:69487/150‑1 (MQ=255)

                    gCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGTGCCCAGCAGGGC‑CAGGCAc                                                                                                                                   <  2:69500/150‑1 (MQ=255)

                        ggCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGTGCCCAGCAGGGC‑CAGGCAcaaca                                                                                                                               <  2:4164075/150‑5 (MQ=255)

                        ggCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGTGCCCAGCAGGGC‑CAGGCAcaaca                                                                                                                               <  2:4869124/150‑5 (MQ=255)

                          cGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTtg                                                                                                                             <  1:3598609/150‑1 (MQ=255)

                           gTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTtgg                                                                                                                            <  1:324179/150‑1 (MQ=255)

                           gTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTtgg                                                                                                                            <  2:1137040/150‑1 (MQ=255)

                             gAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTtggct                                                                                                                          <  1:5021909/150‑1 (MQ=255)

                                gcgcACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAgg                                                                                                                                      <  2:1261860/135‑1 (MQ=255)

                                gcgcACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAgg                                                                                                                                      >  1:1261860/1‑135 (MQ=255)

                                 cgcACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCa                                                                                                                      <  2:93025/150‑1 (MQ=255)

                                 cgcACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCa                                                                                                                      >  1:5161272/1‑150 (MQ=255)

                                   cACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAgg                                                                                                                    >  1:4259114/1‑150 (MQ=255)

                                    aCCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAggg                                                                                                                   >  2:3795671/1‑150 (MQ=255)

                                      cGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGgcg                                                                                                                 >  1:2881103/1‑150 (MQ=255)

                                           gCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAg                                                                                                            >  2:3082398/1‑150 (MQ=255)

                                             aaaGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTg                                                                                                          >  1:2634973/1‑150 (MQ=255)

                                                gCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACg                                                                                                       <  1:2353177/150‑1 (MQ=255)

                                                    ttttCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTtg                                                                                                                             <  1:3818170/124‑1 (MQ=255)

                                                     tttCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGAt                                                                                                  <  2:3910855/150‑1 (MQ=255)

                                                      ttCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATc                                                                                                 <  2:4259114/150‑1 (MQ=255)

                                                       tCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCg                                                                                                >  2:631367/1‑150 (MQ=255)

                                                          aCGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTgg                                                                                             <  1:956997/150‑1 (MQ=255)

                                                           cGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTATTGGCTGGCAGGGCGCTTAGTGACGCAGATCGTGGc                                                                                            <  2:3716722/150‑1 (MQ=255)

                                                               ggaacgccAGTCTGGCCTGCCAGCCTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAAc                                                                                        >  1:4097973/8‑150 (MQ=255)

                                                               cAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAAc                                                                                        <  2:3924011/150‑1 (MQ=255)

                                                               cAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAAc                                                                                        >  1:2218357/1‑150 (MQ=255)

                                                                aGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACt                                                                                       >  2:1450879/1‑150 (MQ=255)

                                                                   agaCAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCg                                                                                    <  1:3691667/150‑1 (MQ=255)

                                                                           tGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAAc                                                                            >  2:5233685/1‑150 (MQ=255)

                                                                                       cTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCaa                                                                <  1:2101747/150‑1 (MQ=255)

                                                                                        tGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGTGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAg                                                               <  1:330467/150‑1 (MQ=255)

                                                                                        tGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAg                                                               <  2:2130249/150‑1 (MQ=255)

                                                                                         ggCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGt                                                              <  1:2547633/150‑1 (MQ=255)

                                                                                          gCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTTGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTg                                                             >  1:5198647/1‑150 (MQ=255)

                                                                                          gCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGa                                                                                          <  2:1642444/121‑1 (MQ=255)

                                                                                          gCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTg                                                             <  2:338118/150‑1 (MQ=255)

                                                                                          gCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTg                                                             >  1:5198611/1‑150 (MQ=255)

                                                                                          gCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTg                                                             <  1:1333330/150‑1 (MQ=255)

                                                                                          gCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCATTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGa                                                                                          >  1:1642444/1‑121 (MQ=255)

                                                                                           cGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTgg                                                            <  1:1119440/150‑1 (MQ=255)

                                                                                            gTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGa                                                           <  2:3154105/150‑1 (MQ=255)

                                                                                             ttCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACAGCGAGGCCAACCTGCAGGCTCAAGTGGat                                                          >  1:1410356/1‑150 (MQ=255)

                                                                                              tCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGata                                                         <  1:5422397/150‑1 (MQ=255)

                                                                                                cTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGAtata                                                       >  1:4379166/1‑150 (MQ=255)

                                                                                                cTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCGTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACAGCGAGGCCAACCTGCAGGCTCAAGTGGAtata                                                       >  1:1923585/1‑150 (MQ=255)

                                                                                                   cGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATAccc                                                    <  1:794382/150‑1 (MQ=255)

                                                                                                      tCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGaa                                                 <  1:342808/150‑1 (MQ=255)

                                                                                                        aGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGa                                               >  2:4040624/1‑150 (MQ=255)

                                                                                                             cTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGaaa                                          <  2:564189/150‑1 (MQ=255)

                                                                                                             cTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGaaa                                          <  1:1807775/150‑1 (MQ=255)

                                                                                                             cTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCATTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGaaa                                          <  1:1046462/150‑1 (MQ=255)

                                                                                                              tCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGa                                                                                   <  2:1937808/108‑1 (MQ=255)

                                                                                                              tCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCCACCTGCAGGCTCAAGTGGATATACCCGAAGATGaaaa                                         >  1:2382540/1‑150 (MQ=255)

                                                                                                              tCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCGTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGa                                                                                   >  1:1937808/1‑108 (MQ=255)

                                                                                                               cGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGaaaaa                                        <  1:4049804/150‑1 (MQ=255)

                                                                                                                     tCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAACAACTCAcc                                  >  1:4811235/1‑150 (MQ=255)

                                                                                                                       gAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACcgc                                >  2:2706197/1‑150 (MQ=255)

                                                                                                                        aaCATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATAGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATAACCGAAGATGAAAAACTCACcgcg                               >  1:4250565/1‑150 (MQ=255)

                                                                                                                         aCATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCgcgc                              <  1:86639/150‑1 (MQ=255)

                                                                                                                              ctgctgGTGGTGCTGTTGTGCCTG‑GCCCGGCTGGGCACAGGCACCGTTTGGCTGGGAGGGCGCTCAGTGACGCAGATCGTGGCGCACTGCGAGGCCAACCAGCAGGCTCAAGTGGATATACCCGCAGATGAAAAACTCACCGCGCGtctc                         >  1:2172/1‑150 (MQ=255)

                                                                                                                                gctgGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCCACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCgg                       >  1:4159457/1‑150 (MQ=255)

                                                                                                                                 ctgGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGata                                                         <  2:591797/115‑1 (MQ=255)

                                                                                                                                 ctgGTGGTGCTGTTGTGCCTG‑GCCCAGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGata                                                         >  1:591797/1‑115 (MQ=255)

                                                                                                                                     tggtgCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAAc                  <  1:464903/150‑1 (MQ=255)

                                                                                                                                        tgCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCTGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCTAGTGGATATACCCGAAGTTGAAAAACTCACCGCGCTTCTCGGCAAACCgg               >  2:3627078/1‑150 (MQ=255)

                                                                                                                                        tgCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCgg               <  2:3751543/150‑1 (MQ=255)

                                                                                                                                        ccaaaaggtgCCTGTGCCCAGCAGGGC‑CAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCgg               >  2:4058839/8‑150 (MQ=255)

                                                                                                                                         gCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGgt              <  2:1403450/150‑1 (MQ=255)

                                                                                                                                           tgttgtGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGgtgt            <  2:2145051/150‑1 (MQ=255)

                                                                                                                                           tgttgtGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGgtgt            >  1:1534787/1‑150 (MQ=255)

                                                                                                                                               gtgCCTGTGCCCAGCAGGGC‑CAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGCAc        <  2:501630/150‑1 (MQ=255)

                                                                                                                                                tgCCTGTGCCCAGCAGGGC‑CAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGCACg       >  1:2410738/1‑150 (MQ=255)

                                                                                                                                                tgCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGCACg       >  1:3218583/1‑150 (MQ=255)

                                                                                                                                                tgCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGCACg       >  1:1621604/1‑150 (MQ=255)

                                                                                                                                                  ccTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGCACGTg     >  1:4217684/1‑150 (MQ=255)

                                                                                                                                                     g‑gCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGCACGTGCCt  >  2:2038696/1‑150 (MQ=255)

GCGCTAGGCTGGGCGGTGCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGCACGTGCCT  >  NC_002947/3728936‑3729234

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