BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,000,788	G→T	5.7%	A112A (GCC→GCA)	ycaC‑IV ←	putative hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,000,788	0	G	T	5.7%	1177.5 / 38.9	351	A112A (GCC→GCA)	ycaC‑IV	putative hydrolase
Reads supporting (aligned to +/- strand): ref base G (164/166); new base T (9/12); total (173/178)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 355 total aligned reads displayed.

TGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTGGCACCAGTGGGCCGTTGGGGCCTTG  >  NC_002947/6000639‑6000937

tGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGACAGCTACGCACACTTCAGTCACAACGCCg                                                                                                                                                       >  1:8737650/1‑150 (MQ=255)

     aCCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCgat                                                                                                                                                  >  1:5536943/1‑150 (MQ=255)

       ctgcctgGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCgatga                                                                                                                                                <  2:7212260/150‑1 (MQ=255)

         gcctgGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATc                                                                                                                                              >  2:1037219/1‑150 (MQ=255)

            tgGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGt                                                                                                                                           >  2:5602576/1‑150 (MQ=255)

                tCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGct                                                                                                                                       >  2:6837281/1‑150 (MQ=255)

                tCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGct                                                                                                                                       >  2:8025814/1‑150 (MQ=255)

                tCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGct                                                                                                                                       >  1:6016204/1‑150 (MQ=255)

                  aTGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGcttc                                                                                                                                     >  1:4731175/1‑150 (MQ=255)

                    gCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGcttctt                                                                                                                                   <  1:6371181/150‑1 (MQ=255)

                      ggTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGc                                                                                                                                 >  2:2209353/1‑150 (MQ=255)

                        tCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCa                                                                                                                               >  2:447252/1‑150 (MQ=255)

                           tGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTg                                                                                                                            <  2:1569896/150‑1 (MQ=255)

                           tGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGAACAGTTGCTTCTTGCCAGTg                                                                                                                            <  2:3955410/150‑1 (MQ=255)

                                   gTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCAc                                                                                                                    >  2:141800/1‑150 (MQ=255)

                                    tCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACt                                                                                                                   <  1:2504142/150‑1 (MQ=255)

                                         ggTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCtt                                                                                                              <  2:847225/150‑1 (MQ=255)

                                            cATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGAc                                                                                                           <  1:7567882/150‑1 (MQ=255)

                                              tGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGa                                                                                                         >  1:5004311/1‑150 (MQ=255)

                                               ttcattGAAGGTTCCCGACGCATCTGTGACGACAAAGACATCGAAACCGGCTTCGATAGCGGACAATGCCGGGAAAGCCACGCAGACTTCGGTCACAACACCTGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGaa                                                                                                        >  1:1030459/5‑150 (MQ=255)

                                               ggCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGaa                                                                                                        <  1:2050050/150‑1 (MQ=255)

                                                  attGAAGGTTCCCGACGCATCTGTGACGACAAAGACATCGAAACCGGCTTCGATAGCGGACAATGCCGGGAAAGCCACGCAGACTTCGGTCACAACACCTGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGAATTc                                                                                                     >  1:2147562/2‑150 (MQ=255)

                                                  attGAAGGTTCCCGACGCATCTGTGACGACAAAGACATCGAAACCGGCTTCGATAGCGGACAATGCCGGGAAAGCCACGCAGACTTCGGTCACAACACCTGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGAATTc                                                                                                     <  2:6315253/149‑1 (MQ=255)

                                                    tGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCtt                                                                                                   <  2:1844390/150‑1 (MQ=255)

                                                    tGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCtt                                                                                                   <  2:5458944/150‑1 (MQ=255)

                                                     gAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTc                                                                                                  >  2:4220440/1‑150 (MQ=255)

                                                     gAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTc                                                                                                  <  1:3212456/150‑1 (MQ=255)

                                                     gAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTc                                                                                                  <  2:9110949/150‑1 (MQ=255)

                                                      aaGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCg                                                                                                 <  2:2645127/150‑1 (MQ=255)

                                                       aGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGt                                                                                                >  2:3572904/1‑150 (MQ=255)

                                                          tACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGt                                                                                             <  2:2796374/150‑1 (MQ=255)

                                                            ccGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTcc                                                                                           >  2:8955348/1‑150 (MQ=255)

                                                            ccGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTcc                                                                                           >  2:8298013/1‑150 (MQ=255)

                                                             ccgacgCATCTGTGACGACAAAGACATCGAAACCGGCTTCGATAGCGGACAATGCCGGGAAAGCAACGCAGACTTCGGTCACAACACATGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGAATTCCTCATTGTcac                                                                                          >  2:4744878/6‑148 (MQ=255)

                                                              cgacgCATCTGTGACGACAAAGACATCGAAACCGGCTTCGATAGCGGACAATGCCGGGAAAGCCACGCAGACTTCGGTCACAACACCTGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGAATTCCTCATTGTCCCa                                                                                         <  1:6581095/146‑1 (MQ=255)

                                                                aGGCATCAGTGACCACGAACACGTCGAACACTTCTTCCAGGGCCGACAGGGCCGGGAAAGCAACGCACACATCAGTCACAACGCCGGCGAGGATCAGTAGCTTCTAGCCAGTGGCCAACACTGCCTAGACGAAGTCTTCGATGTCCCAgg                                                                                       >  2:4946191/1‑150 (MQ=255)

                                                                 ggCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGc                                                                                      <  2:407993/150‑1 (MQ=255)

                                                                   cATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGt                                                                                    <  1:3654974/150‑1 (MQ=255)

                                                                    aTCAGTGACCACGAACACGTCGAACTCTTGTTCGAGGGCCGAGAGGGACTTTAAAGCTACGCTCCCTTGCGTCACAACTCCGGCGATGATGAGGTGGGTCTTGCGACGTGACTTCTCTTCCTTGACGTAGTCTTCGTTGTCCCAGGCGtt                                                                                   <  2:7644962/150‑1 (MQ=255)

                                                                    aTCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGtt                                                                                   >  1:7303127/1‑150 (MQ=255)

                                                                    aTCAGGGACCAAGAACACGTCGAACTCTTCTTCCAGGGCCGACCGGGCCGGGACAGCTAAGCACACATCAGTCACAACGCCGGCGATGATCAGTGGCTTCTTGCCAGTGGCCTGCACTGCCTGGACGACGTCTTCGTTGTCCCAGGCGtt                                                                                   >  2:6366672/1‑150 (MQ=255)

                                                                       aGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGAt                                                                                <  1:1429672/150‑1 (MQ=255)

                                                                          gACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTg                                                                             <  1:3470463/150‑1 (MQ=255)

                                                                           aCCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTgg                                                                            <  2:3715692/150‑1 (MQ=255)

                                                                             cACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGTCCTTCACTGCCTTGACGAAGTCTTCGTTGTCACAGTCGTTGATCTGGcc                                                                          >  2:6540841/1‑150 (MQ=255)

                                                                             cACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGcc                                                                          >  2:3179039/1‑150 (MQ=255)

                                                                              aCGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCa                                                                         >  1:1516623/1‑150 (MQ=255)

                                                                               cGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAg                                                                        >  1:5922893/1‑150 (MQ=255)

                                                                                 aaCACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAggg                                                                      <  1:8071427/150‑1 (MQ=255)

                                                                                 aaCACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAggg                                                                      <  1:1570470/150‑1 (MQ=255)

                                                                                  acacGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAgggc                                                                     >  1:5746350/1‑150 (MQ=255)

                                                                                   cacGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAgggcg                                                                    >  1:8187957/1‑150 (MQ=255)

                                                                                      gTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAgggcgggc                                                                 >  2:1296680/1‑150 (MQ=255)

                                                                                       tCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTc                                                                                                  >  1:1134809/1‑116 (MQ=255)

                                                                                       tCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGAGGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTc                                                                                                  <  2:1134809/116‑1 (MQ=255)

                                                                                         gAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAAt                                                              <  2:2489731/150‑1 (MQ=255)

                                                                                          aaCTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATg                                                             >  2:6238858/1‑150 (MQ=255)

                                                                                           aCTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGt                                                            <  1:4081094/150‑1 (MQ=255)

                                                                                            ctctTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTa                                                           <  1:9133913/150‑1 (MQ=255)

                                                                                             tcttctTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGt                                                                                                      <  1:5492511/106‑1 (MQ=255)

                                                                                             tcttctTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTAc                                                          <  1:5899663/150‑1 (MQ=255)

                                                                                             tcttctTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTAc                                                          >  1:6190484/1‑150 (MQ=255)

                                                                                             acttcttCCATGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTAc                                                          <  2:9145419/149‑1 (MQ=255)

                                                                                              cttcttCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACg                                                         <  2:2142047/150‑1 (MQ=255)

                                                                                              cttcttCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACg                                                         >  1:1741333/1‑150 (MQ=255)

                                                                                                  ttcCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGc                                                     >  1:309745/1‑150 (MQ=255)

                                                                                                  ttcCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGc                                                     >  2:57749/1‑150 (MQ=255)

                                                                                                     cAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTc                                                  >  2:4589065/1‑150 (MQ=255)

                                                                                                        ggCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTc                                                                                                  >  2:4984602/1‑99 (MQ=255)

                                                                                                            gACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAAc                                           >  1:4634116/1‑150 (MQ=255)

                                                                                                             tcAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCATCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGTGGCGTCCTGGAACa                                          <  2:8878854/149‑1 (MQ=255)

                                                                                                              cAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAg                                         <  2:167789/150‑1 (MQ=255)

                                                                                                              cAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAg                                         <  2:5854523/150‑1 (MQ=255)

                                                                                                              cAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAg                                         <  2:3147515/150‑1 (MQ=255)

                                                                                                                gggCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGTTCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAggg                                       <  1:7702963/150‑1 (MQ=255)

                                                                                                                gggCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAggg                                       <  1:8320398/150‑1 (MQ=255)

                                                                                                                 ggCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGc                                      <  1:2023206/150‑1 (MQ=255)

                                                                                                                 ggCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGc                                      <  1:398151/150‑1 (MQ=255)

                                                                                                                  gCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCt                                     >  1:4783930/1‑150 (MQ=255)

                                                                                                                    cGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCttt                                   <  2:5161732/150‑1 (MQ=255)

                                                                                                                      ggAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCa                                 <  1:853338/150‑1 (MQ=255)

                                                                                                                      ggAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCa                                 <  1:6992849/150‑1 (MQ=255)

                                                                                                                           gCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCa                                                                                                                     <  2:8804389/61‑1 (MQ=255)

                                                                                                                           gCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCa                                                                                                                     >  1:8804389/1‑61 (MQ=255)

                                                                                                                           gCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCt                            <  1:1218996/150‑1 (MQ=255)

                                                                                                                            cTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTg                           <  2:3280389/150‑1 (MQ=255)

                                                                                                                            cTACGCACACTTCAGTCACAACGCCGGCGATGATCAGGTGCTTCGTGCCAGTGGCCTGCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTg                           <  2:7022847/150‑1 (MQ=255)

                                                                                                                             tACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTgg                          >  1:5787113/1‑150 (MQ=255)

                                                                                                                             tACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTgg                          >  1:4418827/1‑150 (MQ=255)

                                                                                                                             tACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTgg                          >  1:4335519/1‑150 (MQ=255)

                                                                                                                             tACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTgg                          <  2:7320625/150‑1 (MQ=255)

                                                                                                                             tACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTgg                          >  1:6877475/1‑150 (MQ=255)

                                                                                                                                    acTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTGGCACCAGt                   >  1:9071771/1‑150 (MQ=255)

                                                                                                                                      ttCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGaa                                                                                                        <  2:3435190/63‑1 (MQ=255)

                                                                                                                                      ttCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTGGCACCAGTgg                 <  2:4390050/150‑1 (MQ=255)

                                                                                                                                      ttCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTGGCACCAGTgg                 <  1:4936672/150‑1 (MQ=255)

                                                                                                                                       tCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTGGCACCAGTggg                >  2:2901073/1‑150 (MQ=255)

                                                                                                                                              caACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTGGCACCAGTGGGCCGTTgg         >  2:2450207/1‑150 (MQ=255)

                                                                                                                                                  gCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTGGCACCAGTGGGCCGTTGGGGcc     >  1:798168/1‑150 (MQ=255)

                                                                                                                                                     ggCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTGGCACCAGTGGGCCGTTGGGGCCTTg  >  2:612299/1‑150 (MQ=255)

TGCGCACCTGCCTGGCTCATGCGGTCGTGGGCGGCGTCGCGGGTCATGGCGTTGAAGGTACCGGAGGCATCAGTGACCACGAACACGTCGAACTCTTCTTCCAGGGCCGACAGGGCCGGGAAAGCTACGCACACTTCAGTCACAACGCCGGCGATGATCAGTTGCTTCTTGCCAGTGGCCTTCACTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGGGCGTCCGGGAACAGGGCTTTCAGTTCTGGCACCAGTGGGCCGTTGGGGCCTTG  >  NC_002947/6000639‑6000937

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