BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,297,315	T→G	12.4%	L536R (CTG→CGG)	PP_2893 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,297,315	0	T	G	12.4%	833.4 / 25.9	316	L536R (CTG→CGG)	PP_2893	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (75/179); new base G (42/4); total (132/184)
Fisher's exact test for biased strand distribution p-value = 1.27e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 334 total aligned reads displayed.

GCTCAAGGCCGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTGGATGCCCGCCTGGACA  >  NC_002947/3297166‑3297464

  tCAAGGCCGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAGTCGACCTTGACTTGCCGCGCCTCGGCCAGCtgt                                                                                                                                                     >  2:1671262/1‑150 (MQ=255)

  tCAAGGCCGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCCGCCAGCtgt                                                                                                                                                     >  1:1108147/1‑150 (MQ=255)

  tCAAGGCCGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCACGGCCCGCtgt                                                                                                                                                     >  1:4834463/1‑150 (MQ=255)

    aaGGCCGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTgg                                                                                                                                                   >  2:7041245/1‑150 (MQ=255)

    aaGGCCGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTgg                                                                                                                                                   <  1:3030187/150‑1 (MQ=255)

     aGGCCGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGc                                                                                                                                                  <  2:4169275/150‑1 (MQ=255)

     aGGCCGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGc                                                                                                                                                  <  2:2843192/150‑1 (MQ=255)

     aGGCCGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGc                                                                                                                                                  <  1:6405522/150‑1 (MQ=255)

      ggCCGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACCACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGACTCGACCTTGACTTGCCGCGCCTCGGCCAGCCGTGGcc                                                                                                                                                 <  1:2886369/150‑1 (MQ=255)

       gCCGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCg                                                                                                                                                >  2:4281991/1‑150 (MQ=255)

        ccGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGa                                                                                                                                               <  2:5699155/150‑1 (MQ=255)

        ccGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGa                                                                                                                                               <  2:2999223/150‑1 (MQ=255)

        ccGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGa                                                                                                                                               <  1:4736826/150‑1 (MQ=255)

        cagagaCCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCATCTGTGGCCGa                                                                                                                                               <  2:1632245/148‑1 (MQ=255)

         cGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAg                                                                                                                                              >  2:812608/1‑150 (MQ=255)

         cGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCGCGGCCAGCGGTGGCCGAg                                                                                                                                              >  1:7074320/1‑150 (MQ=255)

          gagaCCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAgg                                                                                                                                             <  2:5443500/150‑1 (MQ=255)

          gagaCCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCACAGCTGGCCGAgg                                                                                                                                             <  2:5178126/150‑1 (MQ=255)

          gagaCCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCAATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAgg                                                                                                                                             <  2:3890594/150‑1 (MQ=255)

             aCCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTa                                                                                                                                          <  1:5381591/150‑1 (MQ=255)

                      gcagggcaAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCgg                                                                                                                                 <  1:3653207/150‑1 (MQ=255)

                      gcagggcaAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCgg                                                                                                                                 <  1:974339/150‑1 (MQ=255)

                       cagggcaAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGGCAGCTGTGGCCGAGGCTACAAGGCCggg                                                                                                                                >  1:7349147/1‑150 (MQ=255)

                          ggcaAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGGCGCGGCCCCCGGTGGCCGCGGCTACAAGGCCGGGTCa                                                                                                                             >  1:3331517/1‑150 (MQ=255)

                          ggcaAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCCGCGGTGGCCGAGGCTACAAGGCCGGGTCa                                                                                                                             >  1:3509611/1‑150 (MQ=255)

                            caAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCACGGCCCCCTGTGGCCGAGGCTACAAGGCCGGGTCAAg                                                                                                                           >  1:6558197/1‑150 (MQ=255)

                             aaaaCTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAgg                                                                                                                          <  1:101242/150‑1 (MQ=255)

                             aaaaCTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAgg                                                                                                                          >  2:5555337/1‑150 (MQ=255)

                             aaaaCTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCGCGGCCCGCGGGGGCCGAGGCGACCAGGCCGGGGCAAgg                                                                                                                          >  1:1684621/1‑150 (MQ=255)

                              aaaCTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCGCGGCCAGCAGTGGCCGAGGGGCCAAGTCCCGGGCAAggg                                                                                                                         >  1:5631812/1‑150 (MQ=255)

                                aCTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGcc                                                                                                                       <  2:3732313/150‑1 (MQ=255)

                                aCTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGcc                                                                                                                       <  2:6913718/150‑1 (MQ=255)

                                   ggACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGc                                                                                                                    <  2:1282667/150‑1 (MQ=255)

                                    gACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGTAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCCGCGGTGGCCGCGGCTACAAGGCCGGGGGCAGGGCCGGcg                                                                                                                   >  1:3156513/1‑149 (MQ=255)

                                    gACCCTGGGCACACTGGACATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCt                                                                                                                   <  2:7324497/150‑1 (MQ=255)

                                     aCCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGTTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTg                                                                                                                  <  1:4452734/150‑1 (MQ=255)

                                       ccTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCCGCGGTGGCCGCGGCTACAAGGCCGGGTCAAGGGCCGGCTGGa                                                                                                                >  2:2145087/1‑150 (MQ=255)

                                        cTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCGACAAGGCCGGGTCAAGGGCCGGCGGGAt                                                                                                               >  1:493934/1‑150 (MQ=255)

                                           ggCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTc                                                                                                            <  1:812608/150‑1 (MQ=255)

                                             cacaCTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCCGCGGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGc                                                                                                          >  2:5835577/1‑150 (MQ=255)

                                             cacaCTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGc                                                                                                          <  2:2009947/150‑1 (MQ=255)

                                                  tGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGGCGCGGCCCCCAGTGGCCGAGGCGACAAGGCCCGGTCAAAGGCCGGCGGGAAGTCGCCGGCa                                                                                                     >  1:3823963/1‑150 (MQ=255)

                                                    gCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCAcc                                                                                                   <  1:4527871/150‑1 (MQ=255)

                                                      cATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCtt                                                                                                 >  2:6643841/1‑150 (MQ=255)

                                                      cATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCtt                                                                                                 <  2:3453466/150‑1 (MQ=255)

                                                      cATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCCGTGGCCGAGGCTACAAGTCCGGGTCAAGGGCCGGCGGGATGTCGCCGGCaacat                                                                                                 >  1:5534007/1‑146 (MQ=255)

                                                      cATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCGCGGCCCCCGGGGGCCGCGGCGACAACGCCGGGGCCACGGCCCGCGGGATGTCGCCCGccccat                                                                                                 >  1:1259872/1‑148 (MQ=255)

                                                         ccAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCa                                                                                              <  2:3565084/150‑1 (MQ=255)

                                                         ccAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCa                                                                                              <  1:1194872/150‑1 (MQ=255)

                                                          cAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAg                                                                                             <  2:12019/150‑1 (MQ=255)

                                                               gggTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCGCGGCCAGCGGTGGCCGCCGCGACAAAGCCCCGGCAAAGGCCGGCAGGATGTACCCCGCAACATGCACGcccc                                                                                        >  1:5837926/1‑150 (MQ=255)

                                                                  tGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGGCACCGCCCGCTGTGGGCGGGGCGACAAGGCCGGGGCAAGGGGCGGCGGGGTGTGGCCCGCCCCGTTCAGGccccccc                                                                                     >  1:2517349/1‑147 (MQ=255)

                                                                   gACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCaa                                                                                    <  2:6791815/150‑1 (MQ=255)

                                                                   gACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCaa                                                                                    <  2:1847422/150‑1 (MQ=255)

                                                                   gACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCaa                                                                                    <  2:1733980/150‑1 (MQ=255)

                                                                   gACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCACCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGGTGTCGCCGGCACCTTGCAGGCCCCGCag                                                                                    >  1:4914900/1‑149 (MQ=255)

                                                                      aaCCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGt                                                                                 <  1:1396261/150‑1 (MQ=255)

                                                                       aCCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTa                                                                                <  1:7039687/150‑1 (MQ=255)

                                                                       aCCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTa                                                                                <  1:2046906/150‑1 (MQ=255)

                                                                        ccGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCCCGGCCCGCGGTGGCCCAGGCGACAATGCCCGGTCAAGGGGCGGCGGGCTGTCGCCCGCACCAAGCACGCCCCCCCAGGTAc                                                                               >  1:5939216/1‑150 (MQ=255)

                                                                         cGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCGCGGCCCCCTGTGGGCGCGGCGCCAAAGCCCGGGGAAAGGCCGGCGGGATGTCGCCGGCACCATGCAGGCCCCCCaacggacc                                                                              >  1:2606434/1‑144 (MQ=255)

                                                                            aTCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCtgtg                                                                                                                                                    >  1:4556508/1‑77 (MQ=255)

                                                                                aCGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCCAGGTACGCTCAccc                                                                       >  2:6635901/1‑150 (MQ=255)

                                                                                 cGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCCGGTCAAGGGCCGGCTGGATGTCGCCGGCACCGTGCAGGCCCCGCAAGGTACGCTCAccca                                                                      >  1:6810832/1‑149 (MQ=255)

                                                                                  ggCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCACAGCTGGCCGAg                                                                                                                                              >  1:905682/1‑77 (MQ=38)

                                                                                    cAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCa                                                                   <  2:1048315/150‑1 (MQ=255)

                                                                                       tGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCGCGGCCAGCGGTGGCCGAGGGTACAAGGCCGGGTCAAGGGCCGGCGGGGTGTCGCCGGCACCGTGCAGGCCCCGCAAGGTACGCTCACCCGGCCggg                                                                >  2:1491479/1‑150 (MQ=255)

                                                                                        ggCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGc                                                               <  1:3226635/150‑1 (MQ=255)

                                                                                         gCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGGCGCGGCCCCCTGTGGGCGCGGCGACAACGCCGGGGCAAGGGCCGGCGGGGTGTGGCCGGCCCCATGCCGGCCCCCCCAGGTACCCGCACCCCGCACGGcc                                                              >  1:6242719/1‑150 (MQ=255)

                                                                                          cAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCTGCAGGGTGAGCGTACCTTGCGGGGcct                                                             <  2:5005639/150‑2 (MQ=255)

                                                                                          cAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCa                                                             <  2:2275737/150‑1 (MQ=255)

                                                                                           aGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGGCTCGGGCCGCTGTGGGCGAGGCTACAACGCCGGGTCCAGGGCCCGCGGGATGTCGCCGGCCCCATGCAGGCCCCGCCAGGGACGCGCACCCAGCAggagatc                                                            >  2:3342063/1‑145 (MQ=255)

                                                                                            gCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAgc                                                           <  1:4044819/150‑1 (MQ=255)

                                                                                               tGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCt                                                        <  2:6040276/150‑1 (MQ=255)

                                                                                               tGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCGGTGGCCGCGGCTACAAGGCCGGGTCAAGGGCCGGCGGGATGTCGCCGGCCCCTTGCCGGCCCCGCAAGGTACGCTCACCCTGCCGGGCCAGCGCt                                                        >  2:4465900/1‑150 (MQ=255)

                                                                                                gCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCtt                                                       <  1:1343938/150‑1 (MQ=255)

                                                                                                 cAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTGGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTg                                                      <  2:6232031/150‑1 (MQ=255)

                                                                                                  aaCAGCGCCTGGCTGGGCGAATCGCCCTTGACGTGCCGCGCCTCGGCCACCGGTGGCCGAGGCTACAAGGCCGGGGCAAGGGCCCGCGGGGTGTCGCCGGCACCGTGCAGGCCCCCCAAGGTACCCTCACCCGGCAGGGCCAgcgcgggg                                                     >  2:4519825/1‑146 (MQ=255)

                                                                                                  aaCAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTgg                                                     <  1:4213159/150‑1 (MQ=255)

                                                                                                   aCAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCCGCTGTGGCCGAGGCGACAAGGCCGGGTCAAGGGCCGGCGGGATGTCGCCGGCACCCTGCAGGCCCCGCCAGGTACGCTCACCCGGCAGGGCCACCGCGTGGc                                                    >  2:763593/1‑150 (MQ=255)

                                                                                                      gcgcCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCGCGGCCAGCGGTGGCCGAGGCGACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCGTGCAGGCCCCCCCAGGTACGCGCACCCGGCAGGGCCCGCGCGTGGCCCa                                                 >  1:1138421/1‑150 (MQ=255)

                                                                                                        gcCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGCTGTCGCCGGCACCTTGCAGGCCCCGCCAGGTACGCTCACCCGGCAGGGCCAGCGCGTGGCCCAgg                                               >  2:5762530/1‑150 (MQ=255)

                                                                                                         ccTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCGCGGCCACCTGTGGCCGAGGCTACAAGGCCCGGGCAAGGGCCGGCGGGAAGTCGCCCGCACCATGCAGGCCCCCCAACGTACGCGCACCCAGCAGGGCCCCCGCGTGGCCCACGc                                              >  1:6070628/1‑150 (MQ=255)

                                                                                                                gggCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAAc                                       <  1:5322124/150‑1 (MQ=255)

                                                                                                                  gCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGGTGTCGCCGGCACCGTGCAGGCCCCGCAAGGTACGCGCACCCAGCAGGGCCAGCGCTTTGCCCAGGCCGAAAACCg                                     >  2:2700402/1‑150 (MQ=255)

                                                                                                                     aaTCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCt                                  <  2:2019450/150‑1 (MQ=255)

                                                                                                                     aaTCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCt                                  <  1:3671727/150‑1 (MQ=255)

                                                                                                                       tCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTgc                                >  2:6187840/1‑150 (MQ=255)

                                                                                                                               gACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTg                        <  1:660464/150‑1 (MQ=255)

                                                                                                                                    gCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCt                   <  1:2779675/150‑1 (MQ=255)

                                                                                                                                      cgcgCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTgg                 <  1:2095317/150‑1 (MQ=255)

                                                                                                                                        cgcCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTGGAt               <  1:237595/150‑1 (MQ=255)

                                                                                                                                            tCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTGGATGccc           >  2:4707222/1‑150 (MQ=255)

                                                                                                                                             cGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTGGATGCCCg          >  2:4053405/1‑150 (MQ=255)

                                                                                                                                                ccAGCTGTGGCCGAGGCTACACGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTGGATGCCCGCCt       <  1:1921781/150‑1 (MQ=255)

                                                                                                                                                ccAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAACACCGCCTGCAGCAGTTGGACCTGGATGCCCGCCt       >  2:4207780/1‑150 (MQ=255)

                                                                                                                                                ccAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTGGATGCCCGCCt       <  1:1671262/150‑1 (MQ=255)

                                                                                                                                                ccAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTGGATGCCCGCCt       <  2:253948/150‑1 (MQ=255)

                                                                                                                                                   gCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTGGATGCCCGCCTGGa    <  2:5419946/150‑1 (MQ=255)

                                                                                                                                                   gCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTGGATGCCCGCCTGGa    <  2:4610271/150‑1 (MQ=255)

GCTCAAGGCCGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTGGATGCCCGCCTGGACA  >  NC_002947/3297166‑3297464

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