BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,357,078	T→A	5.5%	E33D (GAA→GAT)	PP_2950 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,357,078	0	T	A	5.5%	1111.6 / 35.5	326	E33D (GAA→GAT)	PP_2950	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (179/129); new base A (15/3); total (194/132)
Fisher's exact test for biased strand distribution p-value = 4.61e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.85e-01

Only 100 of 329 total aligned reads displayed.

AAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATGGTGGTGCCTCCAATCTTTTGGATGGTTCAGTTGGGTGGGTATTTCTTTAA  >  NC_002947/3356929‑3357227

     ccAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTcgcg                                                                                                                                                  >  1:7934174/1‑150 (MQ=255)

     ccAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTcgcg                                                                                                                                                  >  1:101337/1‑150 (MQ=255)

     ccAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTcgcg                                                                                                                                                  >  2:4892600/1‑150 (MQ=255)

     ccAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTcgcg                                                                                                                                                  >  2:1855301/1‑150 (MQ=255)

     ccAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTcgcg                                                                                                                                                  >  2:420494/1‑150 (MQ=255)

     ccAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTcgcg                                                                                                                                                  <  2:237990/150‑1 (MQ=255)

        gagCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCTGGTTCGCGGGt                                                                                                                                               >  1:6982953/1‑150 (MQ=255)

          gCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTAt                                                                                                                                             <  2:6900060/150‑1 (MQ=255)

            tctTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTc                                                                                                                                           >  1:7437472/1‑150 (MQ=255)

            tctTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTc                                                                                                                                           <  2:7648900/150‑1 (MQ=255)

              tttGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCa                                                                                                                                         >  1:5676858/1‑150 (MQ=255)

              tttGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCa                                                                                                                                         >  1:8149253/1‑150 (MQ=255)

               ttGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAg                                                                                                                                        >  1:7247732/1‑150 (MQ=255)

                tGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGc                                                                                                                                       >  2:2270646/1‑150 (MQ=255)

                tGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGc                                                                                                                                       <  1:7580482/150‑1 (MQ=255)

                 gATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGcc                                                                                                                                      <  1:1230019/150‑1 (MQ=255)

                   tCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCAt                                                                                                                                    >  1:5255898/1‑150 (MQ=255)

                     aaCCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATcc                                                                                                                                  >  1:2875830/1‑150 (MQ=255)

                     aaCCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATcc                                                                                                                                  >  2:8329773/1‑150 (MQ=255)

                      aCCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCt                                                                                                                                 >  2:2091720/1‑150 (MQ=255)

                      aCCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCATTTCGGGTTCGCGGGTATTCCAGCCATCCt                                                                                                                                 >  2:6152350/1‑150 (MQ=255)

                         gATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAg                                                                                                                              >  2:8904025/1‑150 (MQ=255)

                             gTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAg                                                                                                                          >  2:4383047/1‑150 (MQ=255)

                              tCCAGTTCCTTATTCCACTTGCGGGTCAGGAAGGCTTCGGCTTCGGCACGGTTTGTAACGAATTCAACGCGGTTGCGCCACCGCGTGTCCAAGGTATAGGCGAGAGCGACTTTCGCGGGATCGCGGGTATTCCAGCCATCTTCCGCAAgg                                                                                                                         >  1:534974/1‑149 (MQ=255)

                               ccAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTc                                                                                                                        >  2:5882960/1‑150 (MQ=255)

                                cAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCg                                                                                                                       >  1:884670/1‑150 (MQ=255)

                                cAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCg                                                                                                                       >  1:1934175/1‑150 (MQ=255)

                                      tttcttCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCtt                                                                                                                 >  2:8635875/1‑150 (MQ=255)

                                      cttattcCACTTGCGGGTCAGGAAGGCTTCGGCTTCGGCACGGTTTGTAACGAATTCAACGCGGTTGCGCCACCGCGTGTCCAAGGTATAGGCGAGAGCGACTTTCGCGGGATCGCGGGTATTCCAGCCATCTTCCGCAAGGCGGACCtt                                                                                                                 >  2:1425408/5‑150 (MQ=255)

                                      cttattcCACTTGCGGGTCAGGAAGGCTTCGGCTTCGGCACGGTTTGTAACGAATTCAACGCGGTTGCGCCACCGCGTGTCCAAGGTATAGGCGAGAGCGACTTTCGCGGGATCGCGGGTATTCCAGCCATCTTCCGCAAGGCGGACCtt                                                                                                                 >  2:1619703/5‑150 (MQ=255)

                                       ttcttcCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCttt                                                                                                                <  2:2830679/150‑1 (MQ=255)

                                       ttcttcCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCttt                                                                                                                <  1:2664431/150‑1 (MQ=255)

                                       ttattcCACTTGCGGGTCAGGAAGGCTTCGGCTTCGGCACGGTTTGTAACGAATTCAACGCGGTTGCGCCACCGCGTGTCCAAGGTATAGGCGAGAGCGACTTTCGCGGGATCGCGGGTATTCCAGCCATCTTCCGCAAGGCGGACCttt                                                                                                                >  2:3987359/4‑150 (MQ=255)

                                          ttcCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCa                                                                                                             >  2:9088814/1‑150 (MQ=255)

                                          ttcCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCa                                                                                                             >  2:3215354/1‑150 (MQ=255)

                                                ttGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGcc                                                                                                       >  1:4100912/1‑150 (MQ=255)

                                                  gCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGa                                                                                                     >  2:1512059/1‑150 (MQ=255)

                                                   cGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGAt                                                                                                    <  2:5268225/150‑1 (MQ=255)

                                                    gAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGAtt                                                                                                   >  1:139307/1‑150 (MQ=255)

                                                    gAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGAtt                                                                                                   >  1:5477866/1‑150 (MQ=255)

                                                       ttAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTcgc                                                                                                <  2:914336/150‑1 (MQ=255)

                                                       ttAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTcgc                                                                                                >  2:5006492/1‑150 (MQ=255)

                                                       ttAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTcgc                                                                                                <  2:6832411/150‑1 (MQ=255)

                                                       ttAGGAAGCCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTcgc                                                                                                >  1:3164551/1‑150 (MQ=255)

                                                           gAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGt                                                                                            <  1:5171014/150‑1 (MQ=255)

                                                           gAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGt                                                                                            <  1:2279675/150‑1 (MQ=255)

                                                           gAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGt                                                                                            >  2:6693899/1‑150 (MQ=255)

                                                           gAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACATTTTCAATCGCCGATTCGCGAGt                                                                                            >  1:4596436/1‑150 (MQ=255)

                                                               tCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGaaa                                                                                        <  1:1712646/150‑1 (MQ=255)

                                                               tCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGaaa                                                                                        >  2:5455230/1‑150 (MQ=255)

                                                               tCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGaaa                                                                                        <  1:2801528/150‑1 (MQ=255)

                                                               tCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGaaa                                                                                        <  2:4560067/150‑1 (MQ=255)

                                                               tCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGaaa                                                                                        <  2:8740834/150‑1 (MQ=255)

                                                                cccGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAag                                                                                       >  2:6421660/1‑150 (MQ=255)

                                                                cccGAGCTTCTTCACGGTTGTTCGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAag                                                                                       <  2:7902719/150‑1 (MQ=255)

                                                                 ccGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTAt                                                                                                                                             <  1:182504/95‑1 (MQ=255)

                                                                  cGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAagga                                                                                     <  1:4677061/150‑1 (MQ=255)

                                                                  cGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAagga                                                                                     <  1:8494567/150‑1 (MQ=255)

                                                                  cGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAagga                                                                                     <  1:6885362/150‑1 (MQ=255)

                                                                   gAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAaggag                                                                                    <  1:8635875/150‑1 (MQ=255)

                                                                          ttcACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTgg                                                                             >  1:719614/1‑150 (MQ=255)

                                                                          ttcACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTgg                                                                             >  2:3784077/1‑150 (MQ=255)

                                                                           tcACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGc                                                                            >  2:2772944/1‑150 (MQ=255)

                                                                           tcACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGc                                                                            >  2:3310167/1‑150 (MQ=255)

                                                                            cACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCg                                                                           <  2:4599515/150‑1 (MQ=255)

                                                                               ggTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGACGTGGCGGAc                                                                        >  1:4433971/1‑150 (MQ=255)

                                                                               ggTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGAc                                                                        >  1:9035653/1‑150 (MQ=255)

                                                                                     ttGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAgtgt                                                                  <  2:5543512/150‑1 (MQ=255)

                                                                                         aatcTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGTGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTc                                                              <  1:8770772/147‑1 (MQ=255)

                                                                                         aaaCTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTc                                                              >  1:2926160/1‑150 (MQ=255)

                                                                                             tCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTg                                                          >  1:7679383/1‑150 (MQ=255)

                                                                                                gCACGGTTGCGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGAc                                                       >  1:6602155/1‑150 (MQ=255)

                                                                                                      ttACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACAtggtg                                                 <  2:1354169/150‑1 (MQ=255)

                                                                                                      ttACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACAtggtg                                                 <  2:2659993/150‑1 (MQ=255)

                                                                                                             aTTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATGGTGGTGCCTc                                          <  2:7552385/150‑1 (MQ=255)

                                                                                                             aTTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATGGTGGTGCCTc                                          >  2:2898825/1‑150 (MQ=255)

                                                                                                             aTTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATGGTGGTGCCTc                                          <  2:8218471/150‑1 (MQ=255)

                                                                                                                gggTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATGGTGGTGCCTCCaa                                       <  1:4000554/150‑1 (MQ=255)

                                                                                                                 ggTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATGGTGGTGCCTCCAAt                                      >  2:1016465/1‑150 (MQ=255)

                                                                                                                 ggTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCATTGGACATggtggt                                               >  1:8022949/1‑141 (MQ=255)

                                                                                                                     tCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATGGTGGTGCCTCCAATCttt                                  >  2:4791107/1‑150 (MQ=255)

                                                                                                                     tCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATGGTGGTGCCTCCAATCttt                                  >  1:3450560/1‑150 (MQ=255)

                                                                                                                       cACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAggagga                                                                                  >  2:5969025/1‑100 (MQ=255)

                                                                                                                       cACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATGGTGGTGCCTCCAATCTTTTg                                <  1:5017834/150‑1 (MQ=255)

                                                                                                                       cACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATGGTGGTGCCTCCAATCTTTTg                                >  1:487193/1‑150 (MQ=255)

                                                                                                                       cACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATGGTGGTGCCTCCAATCTTTTg                                >  2:7814518/1‑150 (MQ=255)

                                                                                                                       cACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATGGTGGTGCCTCCAATCTTTTg                                >  2:7804418/1‑150 (MQ=255)

                                                                                                                          cGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATGGTGGTGCCTCCAATCTTTTGGAt                             <  2:1906621/150‑1 (MQ=255)

                                                                                                                                ggCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATGGTGGTGCCTCCAATCTTTTGGATGGTTCa                       <  2:1812276/150‑1 (MQ=255)

                                                                                                                                   cAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATGGTGGTGCCTCCAATCTTTTGGATGGTTCAGtt                    <  1:6439108/150‑1 (MQ=255)

                                                                                                                                      ggCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATGGTGGTGCCTCCAATCTTTTGGATGGTTCAGTtggg                 <  1:9153127/150‑1 (MQ=255)

                                                                                                                                       gCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATGGTGGTGCCTCCAATCTTTTGGATGGTTCAGTtgggt                >  2:7273729/1‑150 (MQ=255)

                                                                                                                                       gCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATGGTGGTGCCTCCAATCTTTTGGATGGTTCAGTtgggt                <  1:3539453/150‑1 (MQ=255)

                                                                                                                                       gCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATGGTGGTGCCTCCAATCTTTTGGATGGTTCAGTtgggt                >  2:3437113/1‑150 (MQ=255)

                                                                                                                                             ttttCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATGGTGGTGCCTCCAATCTTTTGGATGGTTCAGTTGGGTGGGTAt          >  1:3346532/1‑150 (MQ=255)

                                                                                                                                                 cGGGTTCGCGGTTATTCCAGCCATCCTCAGCGAGTCTGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATGGTGGTGCCTCCAATCTTTTGGATGGTTCAGTTGGGTGGGTATTTCt      <  2:6982953/150‑1 (MQ=255)

                                                                                                                                                 cGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATGGTGGTGCCTCCAATCTTTTGGATGGTTCAGTTGGGTGGGTATTTCt      <  1:5575767/150‑1 (MQ=255)

                                                                                                                                                    gTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATGGTGGTGCCTCCAATCTTTTGGATGGTTCAGTTGGGTGGGTATTTCTTTa   >  2:4570342/1‑150 (MQ=255)

                                                                                                                                                     ttCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATGGTGGTGCCTCCAATCTTTTGGATGGTTCAGTTGGGTGGGTATTTCTTTaa  <  1:6152350/150‑1 (MQ=255)

                                                                                                                                                     ttCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATGGTGGTGCCTCCAATCTTTTGGATGGTTCAGTTGGGTGGGTATTTCTTTaa  <  1:6299684/150‑1 (MQ=255)

AAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATGGTGGTGCCTCCAATCTTTTGGATGGTTCAGTTGGGTGGGTATTTCTTTAA  >  NC_002947/3356929‑3357227

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