BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,395,612	C→G	5.0%	A254P (GCC→CCC)	PP_2997 ←	protein CnmA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,395,612	0	C	G	5.0%	612.6 / 17.3	179	A254P (GCC→CCC)	PP_2997	protein CnmA
Reads supporting (aligned to +/- strand): ref base C (92/77); new base G (4/5); total (97/82)
Fisher's exact test for biased strand distribution p-value = 7.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

Only 100 of 181 total aligned reads displayed.

TTGGCTGGCTCACTGGAACAATCCGCGAACGACCTCCCTCAACCATCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAGCATGATCCGCCGTTCGATCCCGGCCCGCGCCAATCCCGCTTCCAGCAGCTTCTGCCCG  >  NC_002947/3395463‑3395760

ttGGCTGGCTCACTGGAACAATCCGCGAACGACCTCCCTCAACCATCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGc                                                                                                                                                      >  1:2801667/1‑150 (MQ=255)

  ggctggctCACTTGAACAATCCGCGAACGACCTCCCTCAACCATCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCa                                                                                                                                                    <  1:68087/150‑1 (MQ=255)

  ggctggctCACTGGAACAATCCGCGAACGACCTCCCTCAACCATCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCa                                                                                                                                                    <  1:4008021/150‑1 (MQ=255)

         tCACTGGAACAATCCGCGAACGACCTCCCTCAACCATCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTg                                                                                                                                             >  2:3952312/1‑150 (MQ=255)

         tCACTGGAACAATCCGCGAACGACCTCCCTCAACCATCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTg                                                                                                                                             >  2:1960876/1‑150 (MQ=255)

         tCACTGGAACAATCCGCGAACGACCTCCCTCAACCATCGATTTCCAGTGTCCTGATGATAGCGAGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCCGTGTCTg                                                                                                                                             >  2:268274/1‑150 (MQ=255)

            cTGGAACAATCCGCGAACGACCTCCCTCAACCATCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGccc                                                                                                                                          >  2:3449345/1‑150 (MQ=255)

            cTGGAACAATCCGCGAACGACCTCCCTCAACCATCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTGGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGccc                                                                                                                                          <  2:1512323/150‑1 (MQ=255)

            cTGGAACAATCCGCGAACGACCTCCCTCAACCATCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCGGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATCTCGGCCAGTGTCTGccc                                                                                                                                          >  2:379353/1‑150 (MQ=255)

             tGGAACAATCCGCGAACGACCTCCCTCAACCATCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCg                                                                                                                                         >  2:409542/1‑150 (MQ=255)

               gAACAATCCGCGAACGACCTCCCTCAACCATCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAgt                                                                                                                                                  <  1:2329255/139‑1 (MQ=255)

               gAACAATCCGCGAACGACCTCCCTCAACCATCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAgt                                                                                                                                                  >  2:2329255/1‑139 (MQ=255)

                aaCAATCCGCGAACGACCTCCCTCAACCATCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGata                                                                                                                                      <  1:3952312/150‑1 (MQ=255)

                aaCAATCCGCGAACGACCTCCCTCAACCATCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGata                                                                                                                                      <  2:1988161/150‑1 (MQ=255)

                     tCCGCGAACGACCTCCCTCAACCATCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGGCAGACACTGGCCGATATGcat                                                                                                                                 <  2:1174829/150‑3 (MQ=255)

                     tCCGCGAACGACCTCCCTCAACCATCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCcg                                                                                                                                 >  1:4104574/1‑150 (MQ=255)

                      ccGCGAACGACCTCCCTCAACCATCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGGCAGACACTGGCCGATATGcatg                                                                                                                                <  1:1617521/150‑4 (MQ=255)

                      ccGCGAACGACCTCCCTCAACCATCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCcgc                                                                                                                                >  2:1704997/1‑150 (MQ=255)

                       cgcgAACGACCTCCCTCAACCATCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCcgcg                                                                                                                               <  1:445253/150‑1 (MQ=255)

                       cgcgAACGACCTCCCTCAACCATCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCcgcg                                                                                                                               <  1:3871876/150‑1 (MQ=255)

                        gcgAACGACCTCCCTCAACCATCGATTTCCCTTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCGGTAGCCCATGCATATCGGCCAGTGTCTGCCCGCTATGCCGCgg                                                                                                                              >  2:4305590/1‑150 (MQ=255)

                          gAACGACCTCCCTCAACCATCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCa                                                                                                                            >  2:4389099/1‑150 (MQ=255)

                           aaCGACCTCCCTCAACCATCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCaa                                                                                                                           <  1:1057091/150‑1 (MQ=255)

                              gACCTCCCTCAACCATCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGt                                                                                                                        <  1:3449399/150‑1 (MQ=255)

                              gACCTCCCTCAACCATCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGt                                                                                                                        <  1:2490180/150‑1 (MQ=255)

                                    cctcAACCATCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAg                                                                                                                  >  2:4217215/1‑150 (MQ=255)

                                     ctcACCCATCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGt                                                                                                                 <  2:2875409/150‑1 (MQ=255)

                                      tcAACCATCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTa                                                                                                                >  1:1253175/1‑150 (MQ=255)

                                       cAACCATCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAg                                                                                                               >  2:4354819/1‑150 (MQ=255)

                                          ccATCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTc                                                                                                            <  1:1716953/150‑1 (MQ=255)

                                          ccATCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTc                                                                                                            <  2:442298/150‑1 (MQ=255)

                                           cATCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCg                                                                                                           <  1:3336465/150‑1 (MQ=255)

                                            aTCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGTCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCgg                                                                                                          >  2:1300073/1‑150 (MQ=255)

                                            aTCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCgg                                                                                                          >  2:1093530/1‑150 (MQ=255)

                                            aTCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCgg                                                                                                          >  2:3662673/1‑150 (MQ=255)

                                            aTCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCgg                                                                                                          >  1:1411698/1‑150 (MQ=255)

                                                    ccATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCg                                                                                                  >  1:4394643/1‑150 (MQ=255)

                                                     cATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGa                                                                                                 <  1:3756975/150‑1 (MQ=255)

                                                     cATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGa                                                                                                 <  1:3018183/150‑1 (MQ=255)

                                                     cATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGa                                                                                                 <  1:1275301/150‑1 (MQ=255)

                                                      aTTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGAt                                                                                                >  2:2832600/1‑150 (MQ=255)

                                                            cTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCtt                                                                                          >  2:828014/1‑150 (MQ=255)

                                                             tgatgaTAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTCGGTCGGTAGTCCCGATTATCTTc                                                                                         >  2:3247870/1‑150 (MQ=255)

                                                             tgatgaTAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTc                                                                                         >  1:941951/1‑150 (MQ=255)

                                                             tgatgaTAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTc                                                                                         >  1:3585814/1‑150 (MQ=255)

                                                             tgatgaTAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTc                                                                                         <  1:751460/150‑1 (MQ=255)

                                                              gatgatAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTcc                                                                                        >  2:921344/1‑150 (MQ=255)

                                                              gatgatAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTcc                                                                                        <  2:3615305/150‑1 (MQ=255)

                                                               atgatAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTccc                                                                                       <  1:2997758/150‑1 (MQ=255)

                                                               atgatAGCGGGCGTGCCAATGCTGGTTGACCGTGAATCCCTCCACTGCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTccc                                                                                       <  1:198237/150‑1 (MQ=255)

                                                                 gatAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCa                                                                                     <  1:2870362/150‑1 (MQ=255)

                                                                 gatAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCa                                                                                     >  1:2587744/1‑150 (MQ=255)

                                                                      cGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCg                                                                                >  1:1279290/1‑150 (MQ=255)

                                                                       gggCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCga                                                                               <  2:2983713/150‑1 (MQ=255)

                                                                       gggCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCga                                                                               <  2:1174713/150‑1 (MQ=255)

                                                                       gggCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCga                                                                               >  2:1262434/1‑150 (MQ=255)

                                                                        ggCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTATAGCCCATACATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCgag                                                                              <  1:1078392/150‑1 (MQ=255)

                                                                          cGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGaa                                                                            <  1:3469691/150‑1 (MQ=255)

                                                                            tGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCTGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAAc                                                                          <  1:4070731/150‑1 (MQ=255)

                                                                            tGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAAc                                                                          <  1:1858653/150‑1 (MQ=255)

                                                                            tGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAAc                                                                          <  1:2060401/150‑1 (MQ=255)

                                                                               ctaTGCTGTTTGACCGTGAATCCCTCCACTTCAAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCt                                                                       <  2:4325757/148‑1 (MQ=255)

                                                                                 aTGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTgg                                                                     <  1:3105428/150‑1 (MQ=255)

                                                                                   gCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCa                                                                   <  1:3206586/150‑1 (MQ=255)

                                                                                     tGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGc                                                                 >  2:4497234/1‑150 (MQ=255)

                                                                                     tGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCACGCATATCGGGCAGACACTGGCCGATATGcat                                                                                                                                 <  1:84594/86‑3 (MQ=25)

                                                                                       tttGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTc                                                               <  1:4088419/150‑1 (MQ=255)

                                                                                          gACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAg                                                            <  2:382739/150‑1 (MQ=255)

                                                                                             cGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAGCAt                                                         >  1:3571210/1‑150 (MQ=255)

                                                                                              gTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAGCATg                                                        >  2:4108578/1‑150 (MQ=255)

                                                                                               tGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAGCATGa                                                       >  1:390230/1‑150 (MQ=255)

                                                                                                      cTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAGCATGATccgccg                                                >  2:3743691/1‑150 (MQ=255)

                                                                                                      cTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAGCATGATccgccg                                                >  2:2254472/1‑150 (MQ=255)

                                                                                                       tCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAGCATGATCCGCCGt                                               >  1:1074048/1‑150 (MQ=255)

                                                                                                          aCTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAGCATGATCCGCCGTTCg                                            >  1:1137399/1‑150 (MQ=255)

                                                                                                            ttCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAGCATGATCCGCCGTTCGAt                                          >  1:277909/1‑150 (MQ=255)

                                                                                                             tCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCGTACCCAGCCCGAGAAACGCTGGCAGCTCGAGCATGATCCGCCGTTCGATc                                         >  1:1180685/1‑150 (MQ=255)

                                                                                                                   gggCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAGCATGATCCGCCGTTCGATCCCGGcc                                   >  2:2066168/1‑150 (MQ=255)

                                                                                                                      cAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAGCATGATCCGCCGTTCGATCCCGGCCcgc                                <  2:115631/150‑1 (MQ=255)

                                                                                                                         tCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAGCATGATCCGCCGTTCGATCCCGGCCCGCGcc                             >  2:2477816/1‑150 (MQ=255)

                                                                                                                         tCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAGCATGATCCGCCGTTCGATCCCGGCCCGCGcc                             >  2:2974417/1‑150 (MQ=255)

                                                                                                                          cATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAGCATGATCCGCCGTTCGATCCCGGCCCGCGCCa                            >  1:3711194/1‑150 (MQ=255)

                                                                                                                           aTGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAGCATGATCCGCCGTTCGATCCCGGCCCGCGCCaa                           >  1:4232960/1‑150 (MQ=255)

                                                                                                                             ggACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATTTGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAGCATGATCCGCCGTTCGATCCCGGCCCGCGCCAATc                         <  2:555906/150‑1 (MQ=255)

                                                                                                                                 cTGTAGCCCATGCAAATCGGCCAGTGTCAGCCCGATATGCCGCGGCAAGGTGCCCTGTACGACGGTAGTCCCGATTATCTACCCCAGCCCGTGAAACGCTGCCTCCTCGAGCATGATCCGCCCTACGAACCCGCCCCGCGCCAATCCCGc                     >  2:366801/1‑150 (MQ=255)

                                                                                                                                  tGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAGCATGATCCGCCGTTCGATCCCGGCCCGCGCCAATCCCGCt                    >  1:3502739/1‑150 (MQ=255)

                                                                                                                                       cccATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAGCATGATCCGCCGTTCGATCCCGGCCCGCGCCAATCCCGCTTCcag               <  2:1997541/150‑1 (MQ=255)

                                                                                                                                       cccATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGACCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAGCATGATCCGCCGTTCGATCCCGGCCCGCGCCAATCCCGCTTCcag               >  2:1749346/1‑150 (MQ=255)

                                                                                                                                         cATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAGCATGATCCGCCGTTCGATCCCGGCCCGCGCCAATCCCGCTTCcagca             <  1:4305590/150‑1 (MQ=255)

                                                                                                                                         cATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAGCATGATCCGCCGTTCGATCCCGGCCCGCGCCAATCCCGCTTCcagca             <  1:625176/150‑1 (MQ=255)

                                                                                                                                          cggCATATCGGGCAGACACTGGCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAGCATGATCCGCCGTTCGATCCCGGCCCGCGCCAATCCCGCTTCcagcag            >  1:2381583/3‑150 (MQ=255)

                                                                                                                                           tGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAGCATGATCCGCCGTTCGATCCCGGCCCGCGCCAATCCCGCTTCCagcagc           >  2:1707148/1‑150 (MQ=255)

                                                                                                                                               tatCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAGCATGATCCGCCGTTCGATCCCGGCCCGCGCCAATCCCGCTTCCAGCAGCTTCt       >  1:1402752/1‑150 (MQ=255)

                                                                                                                                                 tCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAGCATGATCCGCCGTTCGATCCCGGCCCGCGCCAATCCCGCTTCCAGCAGCTTCTGc     <  1:3100865/150‑1 (MQ=255)

                                                                                                                                                  cGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAGCATGATCCGCCGTTCGATCCCGGCCCGCGCCAATCCCGCTTCCAGCAGCTTCTGcc    >  1:218799/1‑150 (MQ=255)

                                                                                                                                                  cGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAGCATGATCCGCCGTTCGATCCCGGCCCGCGCCAATCCCGCTTCCAGCAGCTTCTGcc    >  1:218558/1‑150 (MQ=255)

                                                                                                                                                   ggCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAGCATGATCCGCCGTTCGATCCCGGCCCGCGCCAATCCCGCTTCCAGCAGCTTCTGccc   >  1:3103807/1‑150 (MQ=255)

                                                                                                                                                    gCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAGCATGATCCGCCGTTCGATCCCGGCCCGCGCCAATCCCGCTTCCAGCAGCTTCTGCCCg  >  1:524187/1‑150 (MQ=255)

                                                                                                                                                    gCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAGCATGATCCGCCGTTCGATCCCGGCCCGCGCCAATCCCGCTTCCAGCAGCTTCTGCCCg  <  2:1411698/150‑1 (MQ=255)

                                                                                                                                                    gCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAGCATGATCCGCCGTTCGATCCCGGCCCGCGCCAATCCCGCTTCCAGCAGCTTCTGCCCg  >  2:3254204/1‑150 (MQ=255)

TTGGCTGGCTCACTGGAACAATCCGCGAACGACCTCCCTCAACCATCGATTTCCATTGTCCTGATGATAGCGGGCGTGCCAATGCTGTTTGACCGTGAATCCCTCCACTTCGAACGGGCAGTCATGGACCTGTAGCCCATGCATATCGGCCAGTGTCTGCCCGATATGCCGCGGCAAGGTGCCCAGTAGGTCGGTAGTCCCGATTATCTTCCCCAGCCCGAGAAACGCTGGCAGCTCGAGCATGATCCGCCGTTCGATCCCGGCCCGCGCCAATCCCGCTTCCAGCAGCTTCTGCCCG  >  NC_002947/3395463‑3395760

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