BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,365,368	A→T	13.6%	S261T (TCA→ACA)	PP_2075 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,365,368	0	A	T	13.6%	495.0 / 11.2	204	S261T (TCA→ACA)	PP_2075	membrane protein
Reads supporting (aligned to +/- strand): ref base A (119/33); new base T (2/27); total (122/82)
Fisher's exact test for biased strand distribution p-value = 2.08e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 204 total aligned reads displayed.

CAGCACAGCAGCAGGCTGAGCCAGTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCCAGCAGGGCCATTTCCATGTGAGGCAGGCCGGT  >  NC_002947/2365220‑2365515

cagcacagcaGCAGGCTGAGCCAGTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAg                                                                                                                                                    >  2:2677159/1‑150 (MQ=255)

 agcacagcaGCAGGCTGAGCCAGTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAgg                                                                                                                                                   <  2:1958839/150‑1 (MQ=255)

 agcacagcaGCAGGCTGAGCCAGTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAgg                                                                                                                                                   <  2:3355470/150‑1 (MQ=255)

 agcacagcaGCAGGCTGAGCCAGTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCGTCGAACTGGGCGCTGAATGCTTGGATGGGTGAgg                                                                                                                                                   >  1:5223297/1‑150 (MQ=255)

     cagcagcagGCTGAGCCAGTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAg                                                                                                                                               >  1:4810121/1‑150 (MQ=255)

      agcagcagGCTGAGCCAGTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGc                                                                                                                                              >  2:1026319/1‑150 (MQ=255)

        cagcagGCTGAGCCAGTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTg                                                                                                                                            <  2:658477/150‑1 (MQ=255)

         agcagGCTGAGCCAGTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTgg                                                                                                                                           <  2:2887152/150‑1 (MQ=255)

          gcagGCTGAGCCAGTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTggg                                                                                                                                          >  1:3172961/1‑150 (MQ=255)

                tGAGCCAGTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCa                                                                                                                                    >  1:3078236/1‑150 (MQ=255)

                  aGCCAGTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAgg                                                                                                                                  >  2:1695991/1‑150 (MQ=255)

                  aGCCAGTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAgg                                                                                                                                  >  2:2761879/1‑150 (MQ=255)

                     cAGTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAATGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAggcgg                                                                                                                               <  2:90769/150‑1 (MQ=255)

                       gTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAggcgggg                                                                                                                             >  1:2396145/1‑150 (MQ=255)

                       gTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAggcgggg                                                                                                                             <  2:4729090/150‑1 (MQ=255)

                           gCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAggcgggggcgg                                                                                                                         <  1:2072933/150‑1 (MQ=255)

                           gCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAggcgggggcgg                                                                                                                         <  2:3127527/150‑1 (MQ=255)

                                 ccGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACg                                                                                                                   >  2:759686/1‑150 (MQ=255)

                                  cGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACgg                                                                                                                  >  2:2951781/1‑150 (MQ=255)

                                   gatgatGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGc                                                                                                                 >  2:5416692/1‑150 (MQ=255)

                                         gacgacCGCGGGCCAGACAAACAGCTGTGAAACTTTGAACAGCCAGGTGCTTTTGCGTAATAAACGCCAAAGGCCTTCGAACTGGGCGCTGATTGGTTGTTTGGGTGTGGGTAGCTGGGTGCTCTGGGGGGGGGGGGGGACGGCGTGTgg                                                                                                           <  1:2906756/150‑1 (MQ=14)

                                         gaCGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTgg                                                                                                           >  1:2423340/1‑150 (MQ=255)

                                         gaCGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTgg                                                                                                           >  2:4589673/1‑150 (MQ=255)

                                         gaCGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTgg                                                                                                           >  1:4043975/1‑150 (MQ=255)

                                           ctaaggcgGCCCATAAAAACAGCAGTGAAACTTTGAACCGCCTGCTGCTTTTGGGGAAGAACCGCCCACGGCCTTCGAAGTGGGCGGTGAATGGTTGGTTGGGTGTGGGTTGGTGGGTGCCCCGGGGGGGGCGGGGGACGGCGTGTGGGc                                                                                                         <  2:3736172/145‑1 (MQ=255)

                                               cgggggCAAGACAAACAGCAGTGCAACTTTGAACAGCCTGGTGCTTTTGCGGAATAAACGCCCACGGCCTTCTAACTGGGCGCTTAATGCTTGGTTGGGTGAGGGTAGCTGGGTGCTCTGGGGGGGGGGGGGGACGGCGTGTGGGCcttc                                                                                                     <  1:3293390/147‑1 (MQ=14)

                                                 ctggCCAGACAAACAGCATTCAAAATTTGAACAGCCAGCTTCTTTTGCTCAAAAAACCCCAAAGCCCTTCGAACTGGGCGCTGAATTCTTGGTTGGTTTTGTGTAGATGGGTGATAAGGCGGGGGCGGGGCACGGCGTGTGGGCcttctt                                                                                                   <  2:2930593/148‑1 (MQ=9)

                                                  gggCCAGACAACCCGCCGTGCAAGTTTGAACCGCCCGGTGGTTTTGCGGAAGAACCGCCCACGGCCTTCGAACTGGGCGGTGACTGGTGGGTGGGGGGTGGGTCGCTGGGTGCTCTGGGGGGGGGGGGGGACGGCGTGTGGGCCttcttc                                                                                                  <  1:4123491/150‑1 (MQ=255)

                                                    gCCTGACAAACAGCTTTGAAACTTTTAACTGCCTGCTGCTTTTGCGGAATAAACGCCCAAGGCATTCGACCTGGGCGGCGACTGCTTGGTTGGGGGCGGGTAGCTGGTTGCTCTGGGGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCa                                                                                                <  1:1644848/150‑1 (MQ=255)

                                                    gCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCa                                                                                                >  2:3145396/1‑150 (MQ=255)

                                                      cAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAAGTGGGCGCTTAATGCTTGGATGGGTGCGGGTAGGTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCata                                                                                              <  1:3392289/150‑1 (MQ=255)

                                                      cAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCata                                                                                              >  1:1453645/1‑150 (MQ=255)

                                                      cAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCata                                                                                              >  1:2358950/1‑150 (MQ=255)

                                                       aGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCatat                                                                                             >  1:2927279/1‑150 (MQ=255)

                                                           aaaCAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAacg                                                                                         >  2:3688032/1‑150 (MQ=255)

                                                              cagctttGCAATTTTGAACAGCCTGCTTCTTTTGCGGAAGAAACGCCCAAGGCCTTCTAACTGGGCGCTTTATGGTTGGTTGGGTGTGGGTAGCTGGGTGCTCTGGGGGGGGGGGGGGACGGCGTGTGGGCCTTCTTCCATATAacgacg                                                                                      <  1:4047071/144‑1 (MQ=14)

                                                               agcagTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACgg                                                                                     >  1:4003479/1‑150 (MQ=255)

                                                                  agTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATg                                                                                  >  1:4879314/1‑150 (MQ=255)

                                                                    tGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGaa                                                                                >  1:709238/1‑150 (MQ=255)

                                                                      cAACTTTGAACAGCCTGGTGCTTTTGCGTATGAATCGCCAATGGCCTTCTATGTGGGCGCTTATTGCTTGGTTGGGTGTGGGTAGGTGGGTGCTCTGGGGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAAtt                                                                              <  2:5019012/150‑1 (MQ=14)

                                                                      cAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTGGAAGGGGGCGCTTAAGGGTGGGAGGGTTGGGGTTCGCTGGTTCAACTGGCGGGGCCGGGGGAGGCCGTGTGGGCGTTGTACATTTAAACGAGGGATGAAta                                                                              >  1:4180093/1‑149 (MQ=9)

                                                                      cAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAAtt                                                                              >  1:5613962/1‑150 (MQ=255)

                                                                         cTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTccc                                                                           >  1:2218384/1‑150 (MQ=255)

                                                                         cTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTccc                                                                           >  1:1459973/1‑150 (MQ=255)

                                                                          ttttAACAGCCTGCTGTTTTTGCGTATGAAACGCCCACGGCCTTCGAACGGGGGGCTGAATGCTTGGTTGGGTGTGGGTCGCGGGGTGCTCAGGGGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCa                                                                          <  1:3761854/150‑1 (MQ=14)

                                                                          tttGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCa                                                                          >  1:5261365/1‑150 (MQ=255)

                                                                          tttGAACACCCAGCTGCTTTTGCGGAATAAACGCCCAAGGCCTTCGATCTGGGCGCTGAATGCTTGGTTGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGACGAATTCCCa                                                                          <  2:4593555/150‑1 (MQ=14)

                                                                           tttaaCAGCCCGCTGCTTTTGCGGACGAACCGCCCACGGCCTCCGAAGGGGGCGCTGATTGCTTGGTTGGGGGTGGGTTGGTGGGTGCTCTGGGGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCAc                                                                         <  1:899827/147‑1 (MQ=255)

                                                                           ttGAACAGCCCGCTGCTTTTGCGGAAGAAACGCCCACGGCCTTCGAACTGGGCGCTGATTGCTTGGATGGGTGTGGGTAGCTGGGTGCTCAGGGGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCAc                                                                         <  2:4414872/150‑1 (MQ=14)

                                                                             gAACAGCCCGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGATGTGGGCGCTGAATGCTTGGTTGGGTGTGGGTCGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACtc                                                                       <  1:1702893/150‑1 (MQ=14)

                                                                                cagccagcTGCTTTTGCGTAAGAAACGCCCATGGCATTCGAACTGGGCGCTGATTGGTTGGTTGGGTGAGGGTAGGTGGGTGCTCAGGGGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACtcatc                                                                    <  1:528611/150‑1 (MQ=14)

                                                                                  gcccgcTGCTTTTGCTGAAGAACCGCCCAAGGCCTTCGAACTGGGCGCTGAATGGTTGGTTGGGTGCGGGTAGCTGGGTGCTCTGGGGGGGGGGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCAc                                                                  <  2:3374320/150‑1 (MQ=14)

                                                                                   ccagcTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACg                                                                 >  1:1131960/1‑150 (MQ=255)

                                                                                       gtgctTTGGCGTAAGAACCGCCCAGGGCCTTCTAACTGGGGGATGTATGCTTGTAGGGGTGTGGGTCGGTGGGTGCGCCGGGGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGt                                                             <  2:2526068/149‑1 (MQ=14)

                                                                                           ttttGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGc                                                         >  2:5119131/1‑150 (MQ=40)

                                                                                            tttGCGGAAGAAACGCCCACGGCCTCCGAACTGGGCGCTGATTGCTTGGATGGGTGAGGTTTGCTGGGTGCTCCGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCt                                                        <  2:5223300/150‑1 (MQ=14)

                                                                                                 ggAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTgccgc                                                   >  1:735380/1‑150 (MQ=39)

                                                                                                 ggAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTgccgc                                                   >  2:3887457/1‑150 (MQ=39)

                                                                                                 ggAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTgccgc                                                   >  2:3946304/1‑150 (MQ=39)

                                                                                                 ggAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTgccgc                                                   >  2:4091480/1‑150 (MQ=39)

                                                                                                  gaagaaACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTgccgcc                                                  >  1:1346869/1‑150 (MQ=39)

                                                                                                  gaagaaACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTgccgcc                                                  >  1:158735/1‑150 (MQ=39)

                                                                                                      aaaCGCCAAAGGCCTTCGAAGTGGGCGCTGATTGCTTGGTTGGGTGCGGGTCGCGGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCgccttgg                                              <  1:4551641/150‑1 (MQ=14)

                                                                                                       aaCGCCCCCGGCCGGCGACCGGGGCGCGGAAGGCTTGGCTGGGTGCGGGTAGGTGGGGGCTCTGGGGGGGGGGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCgccttggc                                             <  2:4483748/150‑1 (MQ=14)

                                                                                                          gCCCTAGGCCTTCGAAGTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGGGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTAGCTGCCggcttggcgtt                                          >  2:4304872/1‑150 (MQ=14)

                                                                                                          gCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCgccttggcctt                                          >  2:3798012/1‑150 (MQ=39)

                                                                                                           cccAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGTGGGTAGGTGGGTGCTCTGGGGGGGGGGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCgccttggccttg                                         <  2:3574818/150‑1 (MQ=14)

                                                                                                            ccAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCCTATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGt                                        >  1:1732975/1‑150 (MQ=255)

                                                                                                             cAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGtt                                       >  1:38502/1‑150 (MQ=39)

                                                                                                              aaGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGGGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGttt                                      <  2:3459380/150‑1 (MQ=255)

                                                                                                              aaGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGttt                                      >  1:2877104/1‑150 (MQ=39)

                                                                                                              aaGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGttt                                      >  2:5577249/1‑150 (MQ=39)

                                                                                                              aaGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGttt                                      >  1:2899481/1‑150 (MQ=39)

                                                                                                               cggCCTTCGACCGGGGCGCTTATTGCTTGGTTGGGGGCGGGTCGCGGGGTGCTCTGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTg                                     <  2:4731492/149‑1 (MQ=14)

                                                                                                                    ttCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCcag                                >  2:777574/1‑150 (MQ=255)

                                                                                                                    ttCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCcag                                >  2:1539727/1‑150 (MQ=255)

                                                                                                                      cGACCTGGGCGCTGCCTGCTTGGAGGGGTGCGGGTAGCTGGGTGCTCAGGGGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCcagca                              <  2:824639/150‑1 (MQ=14)

                                                                                                                         aCTGGGCGCTGAATGCTTGGTGGGGTGCGGGTCGCTGGGTGCTCAGGGGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCCAGCAggg                           <  1:2605556/150‑1 (MQ=14)

                                                                                                                           tGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCCAGCAGGGcc                         >  1:14879/1‑150 (MQ=255)

                                                                                                                           tGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCCAGCAGGGcc                         >  2:1288197/1‑150 (MQ=255)

                                                                                                                           tGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCCAGCAGGGcc                         >  1:2003936/1‑150 (MQ=255)

                                                                                                                              gcgcTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCCAGCAGGGCCAtt                      >  1:4722127/1‑150 (MQ=255)

                                                                                                                               cgcTGAATGCTTGGTGGGGTGCGGGGCGCGGGGTGCTCCGGGGGGGGGGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCCAGCAGGGCCAttt                     <  1:2951781/150‑1 (MQ=14)

                                                                                                                                   gAATGCTTGGCTGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCCAGCAGGGCCATTTCCAt                 <  2:2435841/150‑1 (MQ=255)

                                                                                                                                    aaTGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCCAGCAGGGCCATTTCCAtg                >  2:3188873/1‑150 (MQ=255)

                                                                                                                                     aTGGTTGGTTGGGTGTGGGTTGGTGGGTGCTCCGGGGGGGGGGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCCAGCAGGGCCATTTCCAtgt               <  2:2396145/150‑1 (MQ=14)

                                                                                                                                     aTGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCCAGCAGGGCCATTTCCAtgt               <  2:5198254/150‑1 (MQ=255)

                                                                                                                                     aTGCTTGCATGGGTGAGGGTAGCTGGGTGCACAGGCGGGGGCGGGCGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCCAGCAGGGCCATTTCCAtgt               <  2:4043975/150‑1 (MQ=255)

                                                                                                                                       gCTTGGCTGGGTGCGGGGAGCTGGGGGCTCTGGGGGGGGGGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCCAGCAGGGCCATTTCCATGTGa             <  1:74554/150‑1 (MQ=14)

                                                                                                                                       gCTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCCAGCAGGGCCATTTCCATGTGa             <  2:1381160/150‑1 (MQ=255)

                                                                                                                                       gCTTGGATGGGTGAGGGTAGCTGGGTGCGCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCCAGCAGGGCCATTTCCATGTGa             <  2:1179058/150‑1 (MQ=255)

                                                                                                                                       gCTGGGTTGGGTGCGGGTCGCGGGGTGCTCTGGGGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCCAGCAGGGCCATTTCCATGTGa             <  1:4169278/150‑1 (MQ=14)

                                                                                                                                        cTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCCAGCAGGGCCATTTCCATGTGag            >  1:695199/1‑150 (MQ=255)

                                                                                                                                         ttGGATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCCAGCAGGGCCATTTCCATGTGagg           >  2:477706/1‑150 (MQ=255)

                                                                                                                                           ggATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCCAGCAGGGCCATTTCCATGTGaggca         >  1:4485498/1‑150 (MQ=255)

                                                                                                                                           ggATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCCAGCAGGGCCATTTCCATGTGaggca         >  1:3567981/1‑150 (MQ=255)

                                                                                                                                           ggATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCCAGCAGGGCCATTTCCATGTGaggca         >  1:307474/1‑150 (MQ=255)

                                                                                                                                            gtggggTGCGGGGCGCGGGGTGCCCCGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCCAGCAGGGCCATTTCCATGTGaggcag        <  1:598673/147‑1 (MQ=14)

                                                                                                                                            gATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCCAGCAGGGCCATTTCCATGTGaggcag        >  1:3417506/1‑150 (MQ=255)

                                                                                                                                               gggTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCCAGCAGGGCCATTTCCATGTGAGGCAGGcc     >  2:5059776/1‑150 (MQ=255)

CAGCACAGCAGCAGGCTGAGCCAGTGGGCGGTGCCGATGATGACGAACGCGGGCCAGACAAACAGCAGTGCAACTTTGAACAGCCAGCTGCTTTTGCGGAAGAAACGCCCAAGGCCTTCGAACTGGGCGCTGAATGCTTGGATGGGTGAGGGTAGCTGGGTGCTCAGGCGGGGGCGGGGGACGGCGTGTGGGCCTTCTTCCATATAACGACGGATGAATTCCCACTCATCACGCAGTTCGCTGCCGCCTTGGCCTTGTTTGCCCAGCAGGGCCATTTCCATGTGAGGCAGGCCGGT  >  NC_002947/2365220‑2365515

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