BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,297,313	G→C	7.9%	Q535H (CAG→CAC)	PP_2893 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,297,313	0	G	C	7.9%	998.6 / 21.5	332	Q535H (CAG→CAC)	PP_2893	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (118/181); new base C (27/4); total (146/186)
Fisher's exact test for biased strand distribution p-value = 3.88e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 338 total aligned reads displayed.

ATGCTCAAGGCCGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTGGATGCCCGCCTGGA  >  NC_002947/3297164‑3297462

aTGCTCAAGGCCGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAg                                                                                                                                                       <  1:94100/150‑1 (MQ=255)

aTGCTCAAGGCCGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAg                                                                                                                                                       <  1:4726514/150‑1 (MQ=255)

    tCAAGGCCGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCCCACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCGCGGCCAGCggt                                                                                                                                                   >  2:4617114/1‑150 (MQ=255)

    tCAAGGCCGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCtgt                                                                                                                                                   >  2:7333816/1‑150 (MQ=255)

    tCAAGGCCGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCtgt                                                                                                                                                   <  2:1252679/150‑1 (MQ=255)

    tCAAGGCCGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCagt                                                                                                                                                   >  1:6797239/1‑150 (MQ=255)

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

    tCAAGGCCGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCGCGGCCAGCggt                                                                                                                                                   >  1:253948/1‑150 (MQ=255)

     cAAGGCCGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCACCtgtg                                                                                                                                                  >  2:6346136/1‑150 (MQ=255)

      aaGGCCGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCCGCCAGCTGTgg                                                                                                                                                 >  1:6385921/1‑150 (MQ=255)

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

       aGGCCGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGc                                                                                                                                                <  2:2843192/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)

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

           cGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCAGTGGGCGAg                                                                                                                                            >  1:3305551/1‑150 (MQ=255)

           cGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCGCCGCCACCGGTGGCCGAg                                                                                                                                            >  1:7074319/1‑150 (MQ=255)

            gagaCCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAgg                                                                                                                                           >  1:2405751/1‑150 (MQ=255)

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

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

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

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

                  cAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACaa                                                                                                                                     <  1:5184453/150‑1 (MQ=255)

                          agggcaAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGt                                                                                                                             <  2:1547304/150‑1 (MQ=255)

                            ggcaAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCGCGGCCACCGGTGGCCCAGGCTACAAGGCCCGGTCa                                                                                                                           >  1:1928320/1‑150 (MQ=255)

                              caAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACCAGGCCGGGTCAAg                                                                                                                         >  2:5777054/1‑150 (MQ=255)

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

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

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

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

                                    tGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCgg                                                                                                                   <  2:3323884/150‑1 (MQ=255)

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

                                       aCCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTg                                                                                                                <  1:221493/150‑1 (MQ=255)

                                         ccTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCGGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGa                                                                                                              >  2:1030125/1‑150 (MQ=255)

                                          cTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCGCGGCCAGCTGTGGCCGAGGCTACAACGCCGGGTCAAGGGCCGGCGGGAt                                                                                                             >  1:4126319/1‑150 (MQ=255)

                                             ggCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCGCGGCCCGCAGTGGCCGCGGCGACAAGGCCCGGGGAAGGGCCGGCGGGGTGTc                                                                                                          >  1:268495/1‑150 (MQ=255)

                                               cacaCTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCCGCGGTGGCCGGGGCGACCAGGCCGGGTCAAGGGCCGGCGGGCTGTCGc                                                                                                        >  2:5835606/1‑150 (MQ=255)

                                               cacaCTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCCGCGGTGGCCGCGGCTACAAGGCCGGGTCAACGGCCGGCGGGATGTCGc                                                                                                        >  2:2084630/1‑150 (MQ=255)

                                                acacTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGcc                                                                                                       <  1:3721561/150‑1 (MQ=255)

                                                 cacTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCg                                                                                                      <  2:5556694/150‑1 (MQ=255)

                                                  acTGGCCATCCAGCTGGGTGACAACCGCATCGACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCCGCGGTGGCCGAGGCTAAAAAGCCGGGTCAAGGGCCGGCTGGCTGTCGcccg                                                                                                     >  1:3892405/1‑148 (MQ=255)

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

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

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

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

                                                          tCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGc                                                                                             <  2:1831814/150‑1 (MQ=255)

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

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

                                                               cTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGcc                                                                                        <  1:3391156/150‑1 (MQ=255)

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

                                                                   gTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTGGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGc                                                                                    <  2:371696/150‑1 (MQ=255)

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

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

                                                                      aCAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCACGGCCAGCTGTGGGCGAGGCTCCCAGGCCGGGTCAAGGGCCGGCGGGGTGTCGCCGGCACCATGCAGGCCCCGCAAg                                                                                 >  2:222746/1‑150 (MQ=255)

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

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

                                                                          ccGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTAc                                                                             >  2:1977357/1‑150 (MQ=255)

                                                                           cGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCGGGGCCTGCAAGGtgcc                                                                            <  1:3671552/150‑4 (MQ=255)

                                                                               tCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCa                                                                        <  1:6857559/150‑1 (MQ=255)

                                                                                cAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCGCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCAc                                                                       <  2:1662154/150‑1 (MQ=255)

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

                                                                                  aCGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCAccc                                                                     >  2:1178191/1‑150 (MQ=255)

                                                                                       aGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAg                                                                <  1:88560/150‑1 (MQ=255)

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

                                                                                          ggCAGCGTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGc                                                             <  2:6226095/150‑1 (MQ=255)

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

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

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

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

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

                                                                                                   cAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTg                                                    <  1:5938998/150‑1 (MQ=255)

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

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

                                                                                                    aaCAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCGGTGGCCGAGGCGACAAGGCCGGGTCAAGGGCCGGCTGGAGGTCGCCGGCCCCGTGCCGGCCCCGCCCGGTACGCTCACCCGGCAGGGCCAGCGCGTgg                                                   >  2:5241709/1‑150 (MQ=255)

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

                                                                                                      cAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGGGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCAGCAGGGCCAGCGCTTGGcc                                                 >  2:147810/1‑150 (MQ=255)

                                                                                                        gcgcCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCa                                               <  1:4937339/150‑1 (MQ=255)

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

                                                                                                          gcCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCCAGGTACGCTCACCCTGCAGGGCCAGCGCGTGGCCCAgg                                             >  2:6649676/1‑150 (MQ=255)

                                                                                                          gcCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAgg                                             <  1:809726/150‑1 (MQ=255)

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

                                                                                                                ctggGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGaaa                                       <  2:2640899/150‑1 (MQ=255)

                                                                                                                   ggCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAAcc                                    <  1:3201265/150‑1 (MQ=255)

                                                                                                                   ggCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAAcc                                    <  1:313404/150‑1 (MQ=255)

                                                                                                                      gAATCGACCTTGACGTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGcc                                 <  1:4617114/150‑1 (MQ=255)

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

                                                                                                                             ccTTGACGTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTgcagca                          <  2:6797239/150‑1 (MQ=255)

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

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

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

                                                                                                                                         tagcCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTGGa              >  1:1702346/3‑150 (MQ=255)

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

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

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

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

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

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

                                                                                                                                                    aGCTGGGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTGGATGCCCGCCTgg   <  1:2856136/150‑1 (MQ=255)

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

ATGCTCAAGGCCGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTGGATGCCCGCCTGGA  >  NC_002947/3297164‑3297462

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