BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,549,929	C→G	9.9%	A96G (GCT→GGT)	groL →	chaperonin GroEL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,549,929	0	C	G	9.9%	611.6 / 45.0	192	A96G (GCT→GGT)	groL	chaperonin GroEL
Reads supporting (aligned to +/- strand): ref base C (87/86); new base G (13/6); total (100/92)
Fisher's exact test for biased strand distribution p-value = 1.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

Only 100 of 195 total aligned reads displayed.

GCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGCCC  >  NC_002947/1549780‑1550078

gcCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGc                                                                                                                                                       >  1:2908653/1‑150 (MQ=255)

   gACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCa                                                                                                                                                    <  1:1867578/150‑1 (MQ=255)

    aCCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAg                                                                                                                                                   <  2:4189310/150‑1 (MQ=255)

     gcATCACCAAGGACGGTGTATCGGTTGCCAAGGAAATCGAACTGGAAGACAAGTTCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGACCGTTCTCGGTCAgg                                                                                                                                                  >  1:4828461/2‑150 (MQ=255)

     gcATCACCAAGGACGGTGTATCGGTTGCCAAGGAAATCGAACTGGAAGACAAGTTCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGACCGTTCTCGGTCAgg                                                                                                                                                  >  1:2093184/2‑150 (MQ=255)

     ccATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAgg                                                                                                                                                  <  1:184726/150‑1 (MQ=255)

          aCCAAGGACGGTGTATCGGTTGCCAAGGAAATCGAACTGGAAGACAAGTTCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGACCGTTCTCGGTCAGGCTATc                                                                                                                                             >  1:1434979/1‑150 (MQ=255)

          aCCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATc                                                                                                                                             <  1:1967786/150‑1 (MQ=255)

          aCCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGGTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATc                                                                                                                                             <  2:2783581/150‑1 (MQ=255)

               ggACGGTGTATCGGTTGCCAAGGAAATCGAACTGGAAGACAAGTTCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGACCGTTCTCGGTCAGGCTATCGttca                                                                                                                                        <  2:1493699/150‑4 (MQ=255)

                gACGGTGTATCGGTTGCCAAGGAAATCGAACTGGAAGACAAGTTCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGACCGTTCTCGGTCAGGCTATCGttcag                                                                                                                                       >  2:4901322/1‑146 (MQ=255)

                gACGGTGTATCGGTTGCCAAGGAAATCGAACTGGAAGACAAGTTCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGACCGTTCTCGGTCAGGCTATCGttcag                                                                                                                                       >  1:1087694/1‑146 (MQ=255)

                   ggCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGaa                                                                                                                                    >  1:1465819/1‑150 (MQ=255)

                    gCGTTTCCGTCGCCAAAGAAATCGAGCTGAACGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAg                                                                                                                                   <  2:1310292/150‑1 (MQ=255)

                         tCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTg                                                                                                                              >  1:4458275/1‑150 (MQ=255)

                          ccGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGGTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCTTCGTTGACGATGGCCTGa                                                                                                                             <  1:2738230/150‑1 (MQ=255)

                            gTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGaaa                                                                                                                           <  1:4020978/150‑1 (MQ=255)

                            gTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGaaa                                                                                                                           <  1:2074508/150‑1 (MQ=255)

                                cctaaGAATTCGAGCTGAAAGTCGCTTTCGAATACCTGGGCGCCCAGCTGGTCTAGGATGTTGCTTCCAAGGCCAACGACGGTGCCGGTGACGGCAGCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCg                                                                                                                       <  2:2002704/147‑1 (MQ=255)

                                 cAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCGCCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGt                                                                                                                      >  2:4522324/1‑150 (MQ=255)

                                 cAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGt                                                                                                                      >  1:4284876/1‑150 (MQ=255)

                                  aaaGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTc                                                                                                                     <  1:4070381/150‑1 (MQ=255)

                                           gAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGc                                                                                                            >  2:4107260/1‑150 (MQ=255)

                                            aGCTGAAAGACGCTTTCGAAAACATGGGCGCCCATCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGTTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCAATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCa                                                                                                           <  1:4672526/150‑1 (MQ=255)

                                             gCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCAt                                                                                                          <  2:2931253/150‑1 (MQ=255)

                                               tGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGa                                                                                                        >  2:4323395/1‑150 (MQ=255)

                                                gAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGaa                                                                                                       >  1:3952303/1‑150 (MQ=255)

                                                     acaagttCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGACCGTTCTCGGTCAGGCTATCGTTCAGGAAGGCGCCAAGGCTGTTGCCGCTGGCATGAAcccaa                                                                                                  >  1:4172471/6‑148 (MQ=255)

                                                     acaagttCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGACCGTTCTCGGTCAGGCTATCGTTCAGGAAGGCGCCAAGGCTGTTGCCGCTGGCATGAAcccaa                                                                                                  >  1:54209/6‑148 (MQ=255)

                                                     aCGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGa                                                                                                  >  2:782035/1‑150 (MQ=255)

                                                     aCGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTCCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGa                                                                                                  >  1:3587997/1‑150 (MQ=255)

                                                       gCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATg                                                                                                >  2:4961161/1‑150 (MQ=255)

                                                       gCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATg                                                                                                <  1:55358/150‑1 (MQ=255)

                                                       gCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATg                                                                                                <  1:782035/150‑1 (MQ=255)

                                                       gCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATg                                                                                                >  2:1853486/1‑150 (MQ=255)

                                                       gCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATg                                                                                                <  1:2744422/150‑1 (MQ=255)

                                                       gCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATg                                                                                                <  2:2341306/150‑1 (MQ=255)

                                                       gCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATg                                                                                                >  2:3197772/1‑150 (MQ=255)

                                                          ttCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGAc                                                                                             >  2:3452918/1‑150 (MQ=255)

                                                          ttCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGACCGTTCTCGGTCAGGCTATCGTTCAGGAAGGCGCCAAGGCTGTTGCCGCTGGCATGaa                                                                                                       <  2:1174108/140‑1 (MQ=255)

                                                          ttCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGACCGTTCTCGGTCAGGCTATCGTTCAGGAAGGCGCCAAGGCTGTTGCCGCTGGCATGaa                                                                                                       >  1:1174108/1‑140 (MQ=255)

                                                          ttCGAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGACCGTTCTCGGTCAGGCTATCGTTCAGGAAGGCGCCAAGGCTGTTGCCGCTGGCATGAACCCAATGGAc                                                                                             >  1:2847167/1‑150 (MQ=255)

                                                           tCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGAcc                                                                                            >  2:2087280/1‑150 (MQ=255)

                                                             gAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTg                                                                                          >  1:1246359/1‑150 (MQ=255)

                                                             gAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAGGGACCTg                                                                                          >  1:932950/1‑150 (MQ=255)

                                                             gAAAACATGGGCGCACAGATGCTGCGCGAAGTGGCTTCGAAGACCAACGACACTGCCGGTGACGGCACCACCACTGCGACCGTTCTCGGTCAGGCTATCGTTCAGGAAGGCGCCAAGGCTGTTGCCGCTGGCATGAACCCAATGGACCTg                                                                                          >  2:1624040/1‑150 (MQ=255)

                                                               aaaCATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGaa                                                                                        <  1:3451200/150‑1 (MQ=255)

                                                                    tGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGcgcg                                                                                   <  2:3971086/150‑1 (MQ=255)

                                                                    tGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGcgcg                                                                                   <  1:3414326/150‑1 (MQ=255)

                                                                     gggCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGc                                                                                                                <  2:4456802/120‑1 (MQ=255)

                                                                         gcCCAGCTGGCCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATc                                                                              <  2:1678546/150‑1 (MQ=255)

                                                                               cTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAg                                                                        <  1:2214659/150‑1 (MQ=255)

                                                                                tGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAgg                                                                       <  2:4078260/150‑1 (MQ=255)

                                                                                tGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAgg                                                                       <  1:2863731/150‑1 (MQ=255)

                                                                                 ggTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGc                                                                      <  2:549646/150‑1 (MQ=255)

                                                                                   tCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCa                                                                                                           <  2:4493767/111‑1 (MQ=255)

                                                                                   tCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCa                                                                                                           >  1:4493767/1‑111 (MQ=255)

                                                                                    cAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCt                                                                                                                               <  2:214223/90‑1 (MQ=255)

                                                                                     aaGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCAct                                                                  <  2:5036373/150‑1 (MQ=255)

                                                                                      aGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCActg                                                                 >  1:4341707/1‑150 (MQ=255)

                                                                                       ggAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCActgc                                                                >  1:3286703/1‑150 (MQ=255)

                                                                                               cTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGt                                                                                                                                           <  1:1729751/67‑1 (MQ=255)

                                                                                               cTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGt                                                                                                                                           >  2:1729751/1‑67 (MQ=255)

                                                                                                ttCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCcgtcgt                                                       >  1:3213694/1‑150 (MQ=255)

                                                                                                ttCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCcgtcgt                                                       >  1:1693006/1‑150 (MQ=255)

                                                                                                ttCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCcgtcgt                                                       >  1:2592348/1‑150 (MQ=255)

                                                                                                ttCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCcgtcgt                                                       >  1:356346/1‑150 (MQ=255)

                                                                                                 tCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGtt                                                      >  2:1955200/1‑150 (MQ=255)

                                                                                                 tCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGtt                                                      <  1:226352/150‑1 (MQ=255)

                                                                                                  ccAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCAc                                                                   <  2:4857979/136‑1 (MQ=255)

                                                                                                   cAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGc                                                    <  1:4107260/150‑1 (MQ=255)

                                                                                                     aGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCg                                                  <  2:4685471/150‑1 (MQ=255)

                                                                                                      ggCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGa                                                 <  1:2199462/150‑1 (MQ=255)

                                                                                                      ggCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGa                                                 <  1:3621565/150‑1 (MQ=255)

                                                                                                       gCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGaa                                                >  2:361273/1‑150 (MQ=255)

                                                                                                        ccAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAAc                                               <  2:2720040/150‑1 (MQ=255)

                                                                                                        ccAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAAc                                               <  1:1232408/150‑1 (MQ=255)

                                                                                                           acgacgCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAATGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGATGCCGAACTga                                            >  1:38475/1‑150 (MQ=255)

                                                                                                            cgacgCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTgaa                                           <  1:5373765/150‑1 (MQ=255)

                                                                                                              acgCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAATGCCACTGCTGCCGTCGTTGCCGAACTgaaga                                         >  1:2459749/1‑150 (MQ=255)

                                                                                                                    ccGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGc                                                                      >  2:4971936/1‑115 (MQ=255)

                                                                                                                    ccGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGc                                                                      <  1:4971735/115‑1 (MQ=255)

                                                                                                                       gTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCa                                <  2:2908653/150‑1 (MQ=255)

                                                                                                                        tGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCaa                               <  1:860358/150‑1 (MQ=255)

                                                                                                                        tGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCaa                               <  1:1492818/150‑1 (MQ=255)

                                                                                                                            ggCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCa                           >  1:4101378/1‑150 (MQ=255)

                                                                                                                               accaccaccGCTTCCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATgc                        <  2:4018257/150‑1 (MQ=255)

                                                                                                                               accaccaccGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATgc                        >  1:3342615/1‑150 (MQ=255)

                                                                                                                                  accaccGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGcc                     <  1:2958613/150‑1 (MQ=255)

                                                                                                                                   ccaccGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCTGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCg                    >  2:1789742/1‑150 (MQ=255)

                                                                                                                                   ccaccGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCg                    >  1:1527855/1‑150 (MQ=255)

                                                                                                                                        gCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTACCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTcc               >  2:4610996/1‑150 (MQ=255)

                                                                                                                                          tACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCaa             <  2:932950/150‑1 (MQ=255)

                                                                                                                                          tACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCaa             <  2:3383552/150‑1 (MQ=255)

                                                                                                                                           aCCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGCATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAg            <  2:4073065/150‑1 (MQ=255)

                                                                                                                                            ccGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAgg           <  1:3614990/150‑1 (MQ=255)

                                                                                                                                              gTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGcc         >  1:3897575/1‑150 (MQ=255)

                                                                                                                                               ttCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTg                                                     >  2:5128357/1‑105 (MQ=255)

                                                                                                                                               ttCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTg                                                     <  1:5128357/105‑1 (MQ=255)

                                                                                                                                                     cTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGccc  <  2:4567916/150‑1 (MQ=255)

GCCGACCATCACCAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGCCC  >  NC_002947/1549780‑1550078

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