BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,908,430	A→T	5.2%	L112Q (CTG→CAG)	PP_4315 ←	epimerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,908,430	0	A	T	5.2%	1060.7 / 29.3	308	L112Q (CTG→CAG)	PP_4315	epimerase
Reads supporting (aligned to +/- strand): ref base A (149/143); new base T (6/10); total (155/153)
Fisher's exact test for biased strand distribution p-value = 3.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 330 total aligned reads displayed.

TCAACGTTTCGAGCAGCATGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGTGGGGGT  >  NC_002947/4908281‑4908578

  aaCGTTTCGAGCAGCATGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCcagc                                                                                                                                                    >  2:6783891/1‑150 (MQ=255)

  aaCGTTTCGAGCAGCATGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCcagc                                                                                                                                                    >  2:3093290/1‑150 (MQ=255)

  aaCGTTTCGAGCAGCATGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCAAGGCcagc                                                                                                                                                    >  2:5987781/1‑150 (MQ=255)

  aaCGTTTCGAGCAGCATGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGGTGTTTGGCGGGGAAGTCCATGGCcagc                                                                                                                                                    >  2:8413044/1‑150 (MQ=255)

           agcagcATGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCTGCTGGCCAtgg                                                                                                                                           <  1:2009673/150‑3 (MQ=255)

           agcagcATGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCa                                                                                                                                           <  2:4081849/150‑1 (MQ=255)

             cagcaTGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCTGCTGGCCAtggac                                                                                                                                         <  2:6188863/150‑5 (MQ=255)

             cagcaTGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCTGCTGGCCAtggac                                                                                                                                         <  2:3276538/150‑5 (MQ=255)

              agcaTGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCTGCTGGCCAtggact                                                                                                                                        <  1:5766085/150‑6 (MQ=255)

               gcaTGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAgcgc                                                                                                                                       >  2:1101904/1‑150 (MQ=255)

                  tGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCAc                                                                                                                                    <  1:2308508/150‑1 (MQ=255)

                  tGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCAc                                                                                                                                    >  1:7415847/1‑150 (MQ=255)

                   ggCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACg                                                                                                                                   >  1:6378225/1‑150 (MQ=255)

                    gCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACgg                                                                                                                                  >  2:2568009/1‑150 (MQ=255)

                      gTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCt                                                                                                                                <  1:6477928/150‑1 (MQ=255)

                       tCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTg                                                                                                                               >  1:3079994/1‑150 (MQ=255)

                         aTCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGAc                                                                                                                             <  1:5137232/150‑1 (MQ=255)

                          tCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGAcg                                                                                                                            >  1:1553919/1‑150 (MQ=255)

                          tCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGCGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGAcg                                                                                                                            >  1:8418448/1‑150 (MQ=255)

                           cGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGAcgc                                                                                                                           <  2:4337937/150‑1 (MQ=255)

                            gATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGAcgcg                                                                                                                          <  1:8938550/150‑1 (MQ=255)

                            gATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGAcgcg                                                                                                                          <  2:3437929/150‑1 (MQ=255)

                              tCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCa                                                                                                                        >  2:6275361/1‑150 (MQ=255)

                                 ccacGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACt                                                                                                                     >  2:2741083/1‑150 (MQ=255)

                                 ccacGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGGTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACt                                                                                                                     >  2:3149001/1‑150 (MQ=255)

                                  cacGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTc                                                                                                                    >  2:4988019/1‑150 (MQ=255)

                                  cacGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTc                                                                                                                    >  2:5561155/1‑150 (MQ=255)

                                  cacGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTc                                                                                                                    >  1:2994609/1‑150 (MQ=255)

                                  cacGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTc                                                                                                                    <  2:8852714/150‑1 (MQ=255)

                                  cacGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTc                                                                                                                    >  1:4919444/1‑150 (MQ=255)

                                  cacGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTc                                                                                                                    >  2:4848651/1‑150 (MQ=255)

                                        aTAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACt                                                                                                              >  1:3417093/1‑150 (MQ=255)

                                        aTAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACt                                                                                                              >  2:4928502/1‑150 (MQ=255)

                                         tAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTg                                                                                                             <  1:2515000/150‑1 (MQ=255)

                                          aaTCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGc                                                                                                            >  1:5572964/1‑150 (MQ=255)

                                          aaTCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGc                                                                                                            >  2:2532883/1‑150 (MQ=255)

                                            tcgtcgCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCgg                                                                                                          <  2:6749126/150‑1 (MQ=255)

                                            tcgtcgCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCgg                                                                                                          <  1:6276351/150‑1 (MQ=255)

                                            tcgtcgCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCgg                                                                                                          <  1:3136374/150‑1 (MQ=255)

                                            tcgtcgATATGGTACACCGCCTTCGCCTGGCTCAAGCACAGCGCCTTCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTTTTTGGCGGTGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCgg                                                                                                          <  1:6848684/150‑1 (MQ=255)

                                             cgtcgCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCggg                                                                                                         >  1:108283/1‑150 (MQ=255)

                                              gtcgCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGt                                                                                                        >  1:7923160/1‑150 (MQ=255)

                                                cgCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTAt                                                                                                      <  1:1601673/150‑1 (MQ=255)

                                                cgCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTAt                                                                                                      <  2:7491632/150‑1 (MQ=255)

                                                cgCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTAt                                                                                                      <  1:9026325/150‑1 (MQ=255)

                                                 gCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTAtt                                                                                                     <  1:3885162/150‑1 (MQ=255)

                                                 gCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTAtt                                                                                                     >  1:8179337/1‑150 (MQ=255)

                                                  cTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTg                                                                                                    >  1:5816404/1‑150 (MQ=255)

                                                  cTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTg                                                                                                    >  1:567894/1‑150 (MQ=255)

                                                   tAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTga                                                                                                   >  1:2539203/1‑150 (MQ=255)

                                                       gTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTgaagcg                                                                                               <  1:4061596/150‑1 (MQ=255)

                                                       gTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTgaagcg                                                                                               >  1:5331948/1‑150 (MQ=255)

                                                        tACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGGTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTgaagcga                                                                                              >  1:4807065/1‑150 (MQ=255)

                                                             cGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCAc                                                                                         >  2:239948/1‑150 (MQ=255)

                                                             cGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCAc                                                                                         >  1:6105563/1‑150 (MQ=255)

                                                               ccTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACtt                                                                                       <  1:6571047/150‑1 (MQ=255)

                                                                     ccTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCt                                                                                 >  2:1277507/1‑150 (MQ=255)

                                                                      cTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTg                                                                                <  1:3323326/150‑1 (MQ=255)

                                                                      cTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTGGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTg                                                                                >  2:8073401/1‑150 (MQ=255)

                                                                       tGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGc                                                                               <  2:2360729/150‑1 (MQ=255)

                                                                         gCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGc                                                                                                                  <  1:5456573/113‑1 (MQ=255)

                                                                          cTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCa                                                                            <  1:3475904/150‑1 (MQ=255)

                                                                           tCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCCCTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCAc                                                                           >  2:7075113/1‑150 (MQ=255)

                                                                            cAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCAcc                                                                          <  2:5227181/150‑1 (MQ=255)

                                                                                 ccAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGtt                                                                     >  2:4471483/1‑150 (MQ=255)

                                                                                 ccAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGtt                                                                     <  1:4464601/150‑1 (MQ=255)

                                                                                  cAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTg                                                                    >  1:917146/1‑150 (MQ=255)

                                                                                   aGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGc                                                                   <  2:5709226/150‑1 (MQ=255)

                                                                                    gcgcCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCt                                                                  <  1:3729959/150‑1 (MQ=255)

                                                                                     cgcCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTc                                                                 <  1:8734358/150‑1 (MQ=255)

                                                                                     cgcCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTc                                                                 <  1:6673308/150‑1 (MQ=255)

                                                                                       ccTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGa                                                               <  2:823875/150‑1 (MQ=255)

                                                                                        cTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGaa                                                              <  2:2055877/150‑1 (MQ=255)

                                                                                         tGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAAc                                                             <  1:6783891/150‑1 (MQ=255)

                                                                                              aaGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCAc                                                        <  2:3615329/150‑1 (MQ=255)

                                                                                                    cGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGc                                                  <  1:2162581/150‑1 (MQ=255)

                                                                                                      cTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCg                                                <  1:5948095/150‑1 (MQ=255)

                                                                                                        ggCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAg                                              <  2:7859691/150‑1 (MQ=255)

                                                                                                                 aCCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTa                                     <  2:773782/150‑1 (MQ=255)

                                                                                                                  ccGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAg                                    >  2:1231587/1‑150 (MQ=255)

                                                                                                                        cGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGACGAGGCCAAGGTAGCATGGc                              >  1:8875526/1‑150 (MQ=255)

                                                                                                                         gggTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTAt                                                                                                      <  2:169933/77‑1 (MQ=255)

                                                                                                                           gttgttTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGc                           >  2:5407929/1‑150 (MQ=255)

                                                                                                                           gttgttTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGc                           >  1:32306/1‑150 (MQ=255)

                                                                                                                               tttGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAg                       <  1:6328020/150‑1 (MQ=255)

                                                                                                                                  ggCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTc                    <  2:7354740/150‑1 (MQ=255)

                                                                                                                                  ggCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTc                    <  2:8824037/150‑1 (MQ=255)

                                                                                                                                  ggCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTc                    <  2:8979740/150‑1 (MQ=255)

                                                                                                                                  ggCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTc                    <  2:6367600/150‑1 (MQ=255)

                                                                                                                                  ggCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTc                    >  2:381863/1‑150 (MQ=255)

                                                                                                                                  ggCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTc                    >  1:7330953/1‑150 (MQ=255)

                                                                                                                                  ggCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTc                    >  1:7035567/1‑150 (MQ=255)

                                                                                                                                     ggggAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGAc                 <  1:7762/150‑1 (MQ=255)

                                                                                                                                           gctgaTGGCCTGCTGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCAc           >  2:2226913/5‑150 (MQ=255)

                                                                                                                                           gctgaTGGCCTGCTGGCCATCAGAGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCTAAGCACTTGCTCCTGCTCACCCAGTAGATCGAACAGCACGAAAGCCGAGGCCAAGGTAGCATGGCAGCACAGGTCGACTTCCAc           >  2:6589316/5‑150 (MQ=255)

                                                                                                                                                tGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGTgg      >  1:1370762/1‑150 (MQ=255)

                                                                                                                                                 ggCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGTggg     >  2:5118464/1‑150 (MQ=255)

                                                                                                                                                  gCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGTgggg    <  1:5077744/150‑1 (MQ=255)

                                                                                                                                                   ccAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGTggggg   <  2:1418182/150‑1 (MQ=255)

                                                                                                                                                    cagcagGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGTGGGGGt  >  1:6262633/1‑150 (MQ=255)

TCAACGTTTCGAGCAGCATGGCGTCATCGATCACCACGACATAATCGTCGCTAAGGTACACCGCCTTCGCCTGGCTCAAGCCCAGCGCCTGCAAAAGGCCCGCTGGCATGTCGACCGCTACGGGTTGTTTGGCGGGGAAGTCCATGGCCAGCAGGCCATCAGCGCCACGGCTGACGCGCAACTCGCCACTGCGGGTATTGAAGCGAAGCACTTGCGCCTGCTCACCCAGTTGCTCGAACAGCACGTAAGCCGAGGCCAAGGTAGCATGGCCGCACAGGTCGACTTCCACGGTGGGGGT  >  NC_002947/4908281‑4908578

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