BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,755,000	A→T	8.8%	V330V (GTA→GTT)	glnA →	glutamine synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,755,000	0	A	T	8.8%	449.2 / 27.0	147	V330V (GTA→GTT)	glnA	glutamine synthetase
Reads supporting (aligned to +/- strand): ref base A (79/55); new base T (6/7); total (85/62)
Fisher's exact test for biased strand distribution p-value = 3.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.43e-01

Only 100 of 148 total aligned reads displayed.

ACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCGGACCCATCGGCCAACCCGTACCTGGCGTTCGCGGCCCTGCTGATGGCTGGCCTG  >  NC_002947/5754851‑5755147

acCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTa                                                                                                                                                     >  1:4928409/1‑150 (MQ=255)

    tctcGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGc                                                                                                                                                 >  2:3810017/1‑150 (MQ=255)

      tcGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTg                                                                                                                                               <  1:3692842/150‑1 (MQ=255)

      tcGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTg                                                                                                                                               <  1:2845382/150‑1 (MQ=255)

               gaatATGCCGGCCTGTCGGAAAACTGCCTGTTCTTCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCTATAAGCGTCTGGTCCCAGGTTATGAAGCTCCGGTTCTGCTGGCCTATTCg                                                                                                                                      >  1:2706268/4‑150 (MQ=255)

                gCTATGCCGTCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCgg                                                                                                                                     >  1:4682483/1‑150 (MQ=255)

                gCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCgg                                                                                                                                     >  2:3142784/1‑150 (MQ=255)

                  tATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCt                                                                                                                                   >  1:4118091/1‑150 (MQ=255)

                   aTGCCGGCCTGTCGGAAAACTGCCTGTTCTTCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCTATAAGCGTCTGGTCCCAGGTTATGAAGCTCCGGTTCTGCTGGCCTATTCGGCTc                                                                                                                                  <  2:4368945/150‑1 (MQ=255)

                   aTGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTc                                                                                                                                  <  1:111699/150‑1 (MQ=255)

                   aTGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTc                                                                                                                                  <  2:3212508/150‑1 (MQ=255)

                   aTGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTc                                                                                                                                  <  1:5096372/150‑1 (MQ=255)

                   aTGCCGCCCTGTCCGACACCGGCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTc                                                                                                                                  <  1:223332/150‑1 (MQ=255)

                       cGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCaa                                                                                                                              >  1:4893180/1‑150 (MQ=255)

                            tGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGtt                                                                                                                         <  2:1350387/150‑1 (MQ=255)

                             gTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTc                                                                                                                        <  1:3086335/150‑1 (MQ=255)

                             gTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGATGGCCTACTCGGCTCGCAACCGTTc                                                                                                                        <  1:2328743/150‑1 (MQ=255)

                              tCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTcc                                                                                                                       >  1:596726/1‑150 (MQ=255)

                                     ccGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGa                                                                                                                <  2:2922177/150‑1 (MQ=255)

                                      cGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGAt                                                                                                               >  2:927809/1‑150 (MQ=255)

                                         ccTGTTCTTCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCTATAAGCGTCTGGTCCCAGGTTATGAAGCTCCGGTTCTGCTGGCCTATTCGGCTCGTAACCGTTCCGCTTcctgccg                                                                                                            <  1:3597516/150‑7 (MQ=255)

                                          cTGTTCTTCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCTATAAGCGTCTGGTCCCAGGTTATGAAGCTCCGGTTCTGCTGGCCTATTCGGCTCGTAACCGTTCCGCTTcctgccgt                                                                                                           >  2:1799262/1‑143 (MQ=255)

                                            gTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTAt                                                                                                         <  2:2922998/150‑1 (MQ=255)

                                              tcTTCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCTATAAGCGTCTGGTCCCAGGTTATGAAGCTCCGGTTCTGCTGGCCTATTCGGCTCGTAACCGTTCCGCTTCCTGCCGTATTc                                                                                                       <  1:1479215/149‑1 (MQ=255)

                                                ttCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCt                                                                                                     >  1:2382176/1‑150 (MQ=255)

                                                ttCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCTATAAGCGTCTGGTCCCAGGTTATGAAGCTCCGGTTCTGCTGGCCTATTCGGCTCGTAACCGTTCCGCTTCCTGCCGTATTcca                                                                                                     >  1:4849105/1‑149 (MQ=255)

                                                    tCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACg                                                                                                 >  2:1246615/1‑150 (MQ=255)

                                                    tCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACg                                                                                                 >  1:3728937/1‑150 (MQ=255)

                                                    tCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCTATAAGCGTCTGGTCCCAGGTTATGAAGCTCCGGTTCTGCTGGCCTATTCGGCTCGTAACCGTTCCGCTTCCTGCCGTATTccattcg                                                                                                 <  1:2726341/150‑6 (MQ=255)

                                                      ggcggcGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCTCTAAGCGTCTGGTCCCAGGTTATGAAGCTCCGGTTCTGCTGGCCTATTCGGCTCGTAACCGTTCCGCTTCCTGCCGTATTccattcggc                                                                                               <  2:2643808/150‑8 (MQ=255)

                                                      ggcggcGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCTATAAGCGTCTGGTCCCAGGTTATGAAGCTCCGGTTCTGCTGGCCTATTCGGCTCGTAACCGTTCCGCTTCCTGCCGTATTccattcggc                                                                                               >  2:3120032/1‑143 (MQ=255)

                                                      ggcggTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGtt                                                                                               <  2:2610682/150‑1 (MQ=255)

                                                              catcaAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGccc                                                                                       <  2:1747332/150‑1 (MQ=255)

                                                               atcaAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCg                                                                                      >  1:794497/1‑150 (MQ=255)

                                                               atcaAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCg                                                                                      >  1:1948090/1‑150 (MQ=255)

                                                                tcaAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGa                                                                                     >  1:2584009/1‑150 (MQ=255)

                                                                tcaAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGa                                                                                     >  1:2964169/1‑150 (MQ=255)

                                                                  aaGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCTCGaaa                                                                                   <  1:4331590/150‑1 (MQ=255)

                                                                  aaGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGaaa                                                                                   <  1:1550705/150‑1 (MQ=255)

                                                                    gCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTACTCCTTACGTTGGCAGCCCGAAAGc                                                                                 >  1:3360346/1‑150 (MQ=255)

                                                                      aCGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGccc                                                                               <  1:4958914/150‑1 (MQ=255)

                                                                        ggTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGa                                                                             <  2:1220317/150‑1 (MQ=255)

                                                                        ggTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGa                                                                             >  2:33574/1‑150 (MQ=255)

                                                                         gTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGAc                                                                            <  2:2870125/150‑1 (MQ=255)

                                                                         gTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGAc                                                                            <  1:3194737/150‑1 (MQ=255)

                                                                          tAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACg                                                                           >  1:1838797/1‑150 (MQ=255)

                                                                          tAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACg                                                                           >  1:1838456/1‑150 (MQ=255)

                                                                              gCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATc                                                                       >  2:3537630/1‑150 (MQ=255)

                                                                                ccTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGa                                                                     >  1:140142/1‑150 (MQ=255)

                                                                                 cTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGaa                                                                    <  1:1738925/150‑1 (MQ=255)

                                                                                 cTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGaa                                                                    <  1:696580/150‑1 (MQ=255)

                                                                                  tGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAg                                                                   <  2:1331829/150‑1 (MQ=255)

                                                                                     aCGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCAc                                                                >  1:1934198/1‑150 (MQ=255)

                                                                                     aCGGCTTCACCAACCCGTCGACCAACTCGGACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCAc                                                                <  2:1876748/150‑1 (MQ=255)

                                                                                        gCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCt                                                             >  1:4405296/1‑150 (MQ=255)

                                                                                         cTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCtt                                                            >  2:272607/1‑150 (MQ=255)

                                                                                         cTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCtt                                                            >  1:3260861/1‑150 (MQ=255)

                                                                                          ttCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTc                                                           >  1:4563676/1‑150 (MQ=255)

                                                                                          ttCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCATTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTc                                                           >  1:194592/1‑150 (MQ=255)

                                                                                             aCCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCg                                                        >  2:123767/1‑150 (MQ=255)

                                                                                              ccAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCgg                                                       >  2:4905795/1‑150 (MQ=255)

                                                                                                 aCCCGTCGACCCACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCGGAcc                                                    >  1:3746977/1‑150 (MQ=255)

                                                                                                   ccGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGAAGCCCGAAAGCCCGAAGTATCGAAGCACGCTTCCCGGACCCa                                                  >  1:209191/1‑150 (MQ=255)

                                                                                                    cGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCGGACCCAt                                                 >  1:2259672/1‑150 (MQ=255)

                                                                                                        gACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCGGACCCATCGGc                                             <  1:913797/150‑1 (MQ=255)

                                                                                                           cAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGtt                                                                                               <  2:1838712/97‑1 (MQ=255)

                                                                                                           cAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCGGACCCATCGGCCaa                                          <  2:4139359/150‑1 (MQ=255)

                                                                                                           cAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCGGACCCATCGGCCaa                                          >  2:1645727/1‑150 (MQ=255)

                                                                                                            aaCTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCGGACCCATCGGCCAAc                                         >  2:2649486/1‑150 (MQ=255)

                                                                                                                    cAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAGGCACGCTTCCCGGACCCATCGGCCAACCCGTACCt                                 >  2:2600660/1‑150 (MQ=255)

                                                                                                                    cAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGCTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACCTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCGGACCCATCGGCCAACCCGTACCt                                 <  2:4931313/150‑1 (MQ=255)

                                                                                                                     aaGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCGGACCCATCGGCCAACCCGTACCTg                                >  2:1690223/1‑150 (MQ=255)

                                                                                                                     aaGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCGGACCCATCGGCCAACCCGTACCTg                                >  2:326515/1‑150 (MQ=255)

                                                                                                                     aaGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCGGACCCATCGGCCAACCCGTACCTg                                >  1:585825/1‑150 (MQ=255)

                                                                                                                       gCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCGGACCCATCGGCCAACCCGTACCTGGc                              >  1:3461652/1‑150 (MQ=255)

                                                                                                                           cTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCGGACCCATCGGCCAACCCGTACCTGGCGTTc                          <  2:2382176/150‑1 (MQ=255)

                                                                                                                              gTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCGGACCCATCGGCCAACCCGTACCTGGCGTTcgcg                       <  2:482366/150‑1 (MQ=255)

                                                                                                                                 ccAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCGGACCCATCGGCCAACCCGTACCTGGCGTTCGCGGcc                    >  2:2407715/1‑150 (MQ=255)

                                                                                                                                 ccAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCGGACCCATCGGCCAACCCGTACCTGGCGTTCGCGGcc                    >  1:4751154/1‑150 (MQ=255)

                                                                                                                                 ccAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCGGACCCATCGGCCAACCCGTACCTGGCGTTCGCGGcc                    >  1:4338098/1‑150 (MQ=255)

                                                                                                                                   aGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCGGACCCATCGGCCAACCCGTACCTGGCGTTCGCGGCCct                  <  1:4290315/150‑1 (MQ=255)

                                                                                                                                   aGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCGGACCCATCGGCCAACCCGTACCTGGCGTTCGCGGCCct                  <  1:4290363/150‑1 (MQ=255)

                                                                                                                                     gCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCGGACCCATCGGCCAACCCGTACCTGGCGTTCGCGGCCctgc                <  2:2611590/150‑1 (MQ=255)

                                                                                                                                      cTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCGGACCCATCGGCCAACCCGTACCTGGCGTTCGCGGCCctgct               <  2:596726/150‑1 (MQ=255)

                                                                                                                                       ttCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCGGACCCATCGGCCAACCCGTACCTGGCGTTCGCGGCCctgctg              >  1:1376496/1‑150 (MQ=255)

                                                                                                                                       ttCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCGGACCCATCGGCCAACCCGTACCTGGCGTTCGCGGCCctgctg              >  2:4657090/1‑150 (MQ=255)

                                                                                                                                        tCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCGGACCCATCGGCCAACCCGTACCTGGCGTTCGCGGCCCTGCTGa             <  1:3241057/150‑1 (MQ=255)

                                                                                                                                        tCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCGGACCCATCGGCCAACCCGTACCTGGCGTTCGCGGCCCTGCTGa             >  2:3425766/1‑150 (MQ=255)

                                                                                                                                        tCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCGGACCCATCGGCCAACCCGTACCTGGCGTTCGCGGCCCTGCTGa             <  1:448582/150‑1 (MQ=255)

                                                                                                                                        tCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTAGTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCGGACCCATCGGCCAACCCGTACCTGGCGTTCGCGGCCCTGCTGa             <  1:448607/150‑1 (MQ=255)

                                                                                                                                         cGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCGGACCCATCGGCCAACCCGTACCTGGCGTTCGCGGCCCTGCTGAt            >  1:2064057/1‑150 (MQ=255)

                                                                                                                                          gAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCGGACCCATCGGCCAACCCGTACCTGGCGTTCGCGGCCCTGCTGAtg           <  1:416701/150‑1 (MQ=255)

                                                                                                                                              ccccGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCGGACCCATCGGCCAACCCGTACCTGGCGTTCGCGGCCCTGCTGAtggctg       >  1:4583648/1‑150 (MQ=255)

                                                                                                                                               cccGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCGGACCCATCGGCCAACCCGTACCTGGCGTTCGCGGCCCTGCTGAtggctgg      <  1:1101259/150‑1 (MQ=255)

                                                                                                                                               cccGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCGGACCCATCGGCCAACCCGTACCTGGCGTTCGCGGCCCTGCTGAtggctgg      >  1:446033/1‑150 (MQ=255)

                                                                                                                                               cccGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCGGACCCATCGGCCAACCCGTACCTGGCGTTCGCGGCCCTGCTGAtggctgg      >  2:5011208/1‑150 (MQ=255)

                                                                                                                                                ccGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCGGACCCATCGGCCAACCCGTACCTGGCGTTCGCGGCCCTGCTGAtggctggc     >  1:1277803/1‑150 (MQ=255)

                                                                                                                                                ccGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCGGACCCATCGGCCAACCCGTACCTGGCGTTCGCGGCCCTGCTGAtggctggc     >  2:3346402/1‑150 (MQ=255)

                                                                                                                                                   gTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCGGACCCATCGGCCAACCCGTACCTGGCGTTCGCGGCCCTGCTGATGGCTGGCCTg  >  2:1069541/1‑150 (MQ=255)

                                                                                                                                                   gTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCGGACCCATCGGCCAACCCGTACCTGGCGTTCGCGGCCCTGCTGATGGCTGGCCTg  >  2:1113627/1‑150 (MQ=255)

ACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCGGACCCATCGGCCAACCCGTACCTGGCGTTCGCGGCCCTGCTGATGGCTGGCCTG  >  NC_002947/5754851‑5755147

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