BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,193,345	A→C	6.0%	E83A (GAG→GCG) ‡	PP_3684 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,193,345	0	A	C	6.0%	934.6 / 33.6	268	E83A (GAG→GCG) ‡	PP_3684	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (126/126); new base C (7/9); total (133/135)
Fisher's exact test for biased strand distribution p-value = 7.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

Only 100 of 277 total aligned reads displayed.

CTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTCCGTATAGACAAG  >  NC_002947/4193196‑4193494

cTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGa                                                                                                                                                       <  2:3732878/150‑1 (MQ=255)

cTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGa                                                                                                                                                       <  2:473538/150‑1 (MQ=255)

 ttCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAg                                                                                                                                                      >  1:993075/1‑150 (MQ=255)

    ccAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGct                                                                                                                                                   >  1:1535823/1‑150 (MQ=255)

    ccAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGct                                                                                                                                                   >  1:6832281/1‑150 (MQ=255)

     cAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGctc                                                                                                                                                  >  1:5208710/1‑150 (MQ=255)

     cAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGctc                                                                                                                                                  >  1:919731/1‑150 (MQ=255)

      aGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGctct                                                                                                                                                 >  1:2522825/1‑150 (MQ=255)

        aCCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCa                                                                                                                                               >  2:4287839/1‑150 (MQ=255)

         ccGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCaa                                                                                                                                              <  2:389703/150‑1 (MQ=255)

          cGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCTCACGGCTCTCAAg                                                                                                                                             <  1:2912028/150‑1 (MQ=255)

          cGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAg                                                                                                                                             <  2:2759446/150‑1 (MQ=255)

          cGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAg                                                                                                                                             >  2:5884213/1‑150 (MQ=255)

           gAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGc                                                                                                                                            >  2:6801613/1‑150 (MQ=255)

              gACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACt                                                                                                                                         <  1:6576646/150‑1 (MQ=255)

               aCCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACtt                                                                                                                                        <  2:5433438/150‑1 (MQ=255)

                ccAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTg                                                                                                                                       <  1:2299092/150‑1 (MQ=255)

                 cAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTgg                                                                                                                                      >  2:5201397/1‑150 (MQ=255)

                 cAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTgg                                                                                                                                      >  2:3304666/1‑150 (MQ=255)

                 cAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTgg                                                                                                                                      <  1:6777348/150‑1 (MQ=255)

                 cAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTgg                                                                                                                                      >  2:977819/1‑150 (MQ=255)

                  aGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGt                                                                                                                                     >  1:799832/1‑150 (MQ=255)

                  aGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGt                                                                                                                                     <  1:3687173/150‑1 (MQ=255)

                   gAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTc                                                                                                                                    <  1:4028987/150‑1 (MQ=255)

                   gAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTc                                                                                                                                    >  1:6523457/1‑150 (MQ=255)

                      gCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGAc                                                                                                                                 >  1:3494657/1‑150 (MQ=255)

                          agaTAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCgg                                                                                                                             <  1:3054571/150‑1 (MQ=255)

                            aTAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGgcg                                                                                                                           <  2:4789034/150‑1 (MQ=255)

                            aTAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGgcg                                                                                                                           <  2:1430718/150‑1 (MQ=255)

                              aaGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCa                                                                                                                         <  2:3512001/150‑1 (MQ=255)

                                    tGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGaa                                                                                                                   >  1:468172/1‑150 (MQ=255)

                                       gcgcCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGc                                                                                                                >  2:240300/1‑150 (MQ=255)

                                       gcgcCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGc                                                                                                                >  1:3308016/1‑150 (MQ=255)

                                        cgcCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCg                                                                                                               <  1:2337786/150‑1 (MQ=255)

                                         gcCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCTCACGGCTCTCAAg                                                                                                                                             >  1:4149143/1‑119 (MQ=255)

                                         gcCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCTCACGGCTCTCAAg                                                                                                                                             <  2:4149143/119‑1 (MQ=255)

                                               tCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACaaaa                                                                                                        <  2:3360781/150‑1 (MQ=255)

                                                cATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAAc                                                                                                       <  2:5758293/150‑1 (MQ=255)

                                                  tCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCt                                                                                                     >  1:1966742/1‑150 (MQ=255)

                                                  tCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCt                                                                                                     <  1:4489292/150‑1 (MQ=255)

                                                    gAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGACCAAAACCTGa                                                                                                   >  2:4560264/1‑150 (MQ=255)

                                                    gAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGa                                                                                                   >  1:6344141/1‑150 (MQ=255)

                                                    gAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGa                                                                                                   <  2:6894875/150‑1 (MQ=255)

                                                    gAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGa                                                                                                   >  1:185089/1‑150 (MQ=255)

                                                       gCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAgcg                                                                                                >  1:5923918/1‑150 (MQ=255)

                                                             aaCAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCtt                                                                                          <  1:5031122/150‑1 (MQ=255)

                                                                aGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTtgtg                                                                                       <  2:4871458/150‑1 (MQ=255)

                                                                aGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTtgtg                                                                                       <  1:1022999/150‑1 (MQ=255)

                                                                   cTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGtt                                                                                    >  2:4092957/1‑150 (MQ=255)

                                                                   cTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGtt                                                                                    <  1:6686282/150‑1 (MQ=255)

                                                                    tCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCTCACGGCTCTCAAg                                                                                                                                             >  1:2493274/1‑92 (MQ=39)

                                                                    tCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCTCACGGCTCTCAAg                                                                                                                                             <  2:2493274/92‑1 (MQ=39)

                                                                    tCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGttt                                                                                   >  1:7027916/1‑150 (MQ=255)

                                                                    tCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGttt                                                                                   >  1:5373195/1‑150 (MQ=255)

                                                                    tCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCGTGTGGttt                                                                                   >  2:3292130/1‑150 (MQ=255)

                                                                     cAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTc                                                                                  >  1:5809692/1‑150 (MQ=255)

                                                                      aaaaCCTTGAGTGTGACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATGCTTGAGAGCCGTGAGGCTCTATAGGACTTGGTCGACCCGGCGCATCGGAAAGCGTACAAAACCTGAGCGCTGCTTGTGTTTTCa                                                                                 <  1:5650545/150‑1 (MQ=255)

                                                                      aaaaCCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCa                                                                                 >  1:1192420/1‑150 (MQ=255)

                                                                       aaaCCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAg                                                                                <  2:4024018/150‑1 (MQ=255)

                                                                       aaaCCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAg                                                                                <  2:4759230/150‑1 (MQ=255)

                                                                       aaaCCTAGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGAGAGGCTCACAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTATTGGTTTCAg                                                                                >  1:6193810/1‑150 (MQ=255)

                                                                         cccTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGAt                                                                              <  2:6747851/149‑1 (MQ=255)

                                                                             tGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCa                                                                          <  2:559940/150‑1 (MQ=255)

                                                                              gtgggTCACTGACCGAGTCGGGATGTCGATCGACCCTGAAGAGTTACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCaa                                                                         <  2:1224592/148‑1 (MQ=255)

                                                                              gAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCaa                                                                         <  2:7310405/150‑1 (MQ=255)

                                                                              gAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCaa                                                                         <  1:2835392/150‑1 (MQ=255)

                                                                               aGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAg                                                                        <  2:2041920/150‑1 (MQ=255)

                                                                                gggTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGc                                                                       >  1:551856/1‑150 (MQ=255)

                                                                                 ggTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGcc                                                                      <  2:2986203/150‑1 (MQ=255)

                                                                                 ggTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGcc                                                                      <  2:6636747/150‑1 (MQ=255)

                                                                                   tCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGgc                                                                                                                            >  1:4796540/1‑94 (MQ=255)

                                                                                     aCTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCctct                                                                  >  1:3413194/1‑150 (MQ=255)

                                                                                     aCTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCctct                                                                  <  2:6054695/150‑1 (MQ=255)

                                                                                         aCCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAAc                                                              >  1:7240866/1‑150 (MQ=255)

                                                                                         aCCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAAc                                                              >  1:364342/1‑150 (MQ=255)

                                                                                           cGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTg                                                            >  2:1174397/1‑150 (MQ=255)

                                                                                             aGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCa                                                          <  2:2555820/150‑1 (MQ=255)

                                                                                             aGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCa                                                          <  2:7391711/150‑1 (MQ=255)

                                                                                             aGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCa                                                          <  1:1575840/150‑1 (MQ=255)

                                                                                                cGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTTAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCg                                                       >  1:6871419/1‑150 (MQ=255)

                                                                                                   gATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACt                                                    >  1:4777271/1‑150 (MQ=255)

                                                                                                         gatcgaCCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACt                                              >  2:810408/1‑150 (MQ=255)

                                                                                                          atcgaCCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTg                                             >  2:5104574/1‑150 (MQ=255)

                                                                                                                ccTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACAGCGTAtt                                       >  2:7178716/1‑150 (MQ=255)

                                                                                                                     agagATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTtgtg                                  <  1:1760677/150‑1 (MQ=255)

                                                                                                                       agaTACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTTAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGcc                                <  2:6871419/150‑1 (MQ=255)

                                                                                                                        gaTACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGccc                               <  1:4880508/150‑1 (MQ=255)

                                                                                                                            cGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTg                           <  1:7133131/150‑1 (MQ=255)

                                                                                                                             gtgatgtgatCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTgc                          >  2:5422490/1‑150 (MQ=255)

                                                                                                                                 tgtgatCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTc                      <  2:2962306/150‑1 (MQ=255)

                                                                                                                                 tgtgatCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTc                      <  2:551856/150‑1 (MQ=255)

                                                                                                                                 tgtgatCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTc                      <  1:6964814/150‑1 (MQ=255)

                                                                                                                                 tgtgatCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTc                      <  1:410672/150‑1 (MQ=255)

                                                                                                                                        cTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACt               >  1:434086/1‑150 (MQ=255)

                                                                                                                                        cTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACt               >  1:544234/1‑150 (MQ=255)

                                                                                                                                          tGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTcc             <  1:2054038/150‑1 (MQ=255)

                                                                                                                                           gagaGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTCCg            >  2:1059738/1‑150 (MQ=255)

                                                                                                                                              agCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTCCGtat         >  2:1781709/1‑150 (MQ=255)

                                                                                                                                                 cGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTCCGTATAGa      >  1:6926005/1‑150 (MQ=255)

                                                                                                                                                    gAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTCCGTATAGACaa   <  2:6398844/150‑1 (MQ=255)

CTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTCCGTATAGACAAG  >  NC_002947/4193196‑4193494

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