BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,046,140	T→A	6.2%	I15F (ATC→TTC) ‡	PP_5293 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,046,140	0	T	A	6.2%	1142.2 / 42.1	348	I15F (ATC→TTC) ‡	PP_5293	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (162/164); new base A (7/15); total (169/179)
Fisher's exact test for biased strand distribution p-value = 1.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 355 total aligned reads displayed.

AAAGTTCAGCGCCTCTTTGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCAAAGCGAGCATGGCGTTGAAATCTTC  >  NC_002947/6045991‑6046289

   gTTCAGCGCCTCTTTGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTc                                                                                                                                                    >  2:6272272/1‑150 (MQ=255)

   gTTCAGCGCCTCTTTGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTc                                                                                                                                                    <  1:2023486/150‑1 (MQ=255)

     tCAGCGCCTCTTTGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCgg                                                                                                                                                  >  1:549186/1‑150 (MQ=255)

        gcgcCTCTTTGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCAt                                                                                                                                               <  1:1308353/150‑1 (MQ=255)

        gcgcCTCTTTGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCAATGCCGAAATCCGGCaa                                                                                                                                               >  2:3460042/1‑149 (MQ=255)

         cgcCTCTTTGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCAtt                                                                                                                                              >  2:5581768/1‑150 (MQ=255)

          gcCTCTTTGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCGGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTg                                                                                                                                             >  2:5673008/1‑150 (MQ=255)

          gcCTCTTTGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTg                                                                                                                                             >  2:5950405/1‑150 (MQ=255)

          gcCTCTTTGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTg                                                                                                                                             >  2:2357092/1‑150 (MQ=255)

            ctctTTGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCTGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCAATGCCGATATCCGGCAATggg                                                                                                                                           <  2:1465973/150‑1 (MQ=255)

            ctctTTGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCAATGCCGAAATCCGGCAATggg                                                                                                                                           <  2:5245676/150‑1 (MQ=255)

            ctctTTGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCAATGCCGAAATCCGGCAATggg                                                                                                                                           <  1:1663752/150‑1 (MQ=255)

            ctctTTGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCAATGCCGAAATCCGGCAATggg                                                                                                                                           <  1:5128685/150‑1 (MQ=255)

             tctTTGCCCTGGGCGTCGATAACAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGc                                                                                                                                          <  2:3794739/150‑1 (MQ=255)

             tctTTGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGc                                                                                                                                          >  1:3057189/1‑150 (MQ=255)

                ttGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCggg                                                                                                                                       <  1:2055811/150‑1 (MQ=255)

                ttGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCggg                                                                                                                                       <  2:3600183/150‑1 (MQ=255)

                  gCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGt                                                                                                                                     >  2:3755286/1‑150 (MQ=255)

                    ccTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGa                                                                                                                                   <  2:3253718/150‑1 (MQ=255)

                     cTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGAt                                                                                                                                  >  2:3899480/1‑150 (MQ=255)

                      tGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATc                                                                                                                                 <  1:1205948/150‑1 (MQ=255)

                      tGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATc                                                                                                                                 >  2:3991750/1‑150 (MQ=255)

                         gCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGAc                                                                                                                              <  1:4627803/150‑1 (MQ=255)

                         gCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGAc                                                                                                                              <  1:6418507/150‑1 (MQ=255)

                          cGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACa                                                                                                                             >  2:6772164/1‑150 (MQ=255)

                           gTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAg                                                                                                                            <  2:2560262/150‑1 (MQ=255)

                             cGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGt                                                                                                                          <  1:5306096/150‑1 (MQ=255)

                              gATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGtt                                                                                                                         >  1:1041091/1‑150 (MQ=255)

                              gATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGtt                                                                                                                         >  1:5630606/1‑150 (MQ=255)

                                tAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCg                                                                                                                       >  1:3337082/1‑150 (MQ=255)

                                  aaaGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGaa                                                                                                                     >  1:1784622/1‑150 (MQ=255)

                                  aaaGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGaa                                                                                                                     <  1:6495177/150‑1 (MQ=255)

                                        tCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTc                                                                                                               >  2:3309579/1‑150 (MQ=255)

                                             cTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGcc                                                                                                          >  2:3334438/1‑150 (MQ=255)

                                              tCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGACTCGCTCATGccc                                                                                                         >  1:4131377/1‑150 (MQ=255)

                                              tCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGccc                                                                                                         >  2:5308080/1‑150 (MQ=255)

                                              tCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGccc                                                                                                         >  1:2390449/1‑150 (MQ=255)

                                                gcgcTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCg                                                                                                       <  1:3865186/150‑1 (MQ=255)

                                                gcgcTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCg                                                                                                       >  1:3667080/1‑150 (MQ=255)

                                                   cTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCAc                                                                                                    <  2:1807833/150‑1 (MQ=255)

                                                       tggTGCTATGGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCAcctcc                                                                                                >  1:2880390/2‑150 (MQ=255)

                                                       aGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCAcctcc                                                                                                <  2:314748/150‑1 (MQ=255)

                                                          tGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTg                                                                                             >  2:4281732/1‑150 (MQ=255)

                                                          tGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTg                                                                                             <  2:5874561/150‑1 (MQ=255)

                                                           gCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTgg                                                                                            >  1:2162517/1‑150 (MQ=255)

                                                           gCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTgg                                                                                            >  1:5654435/1‑150 (MQ=255)

                                                             cATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGcc                                                                                          <  1:6222843/150‑1 (MQ=255)

                                                              aTAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCg                                                                                         >  2:6851468/1‑150 (MQ=255)

                                                                aGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGt                                                                                       <  2:4129791/150‑1 (MQ=255)

                                                                  cACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGttt                                                                                     >  1:3264584/1‑150 (MQ=255)

                                                                   aCCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTtc                                                                                    >  1:3289364/1‑150 (MQ=255)

                                                                    ccAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTtca                                                                                   <  2:4211457/150‑1 (MQ=255)

                                                                    ccAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTtca                                                                                   <  2:2215476/150‑1 (MQ=255)

                                                                    ccAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTtca                                                                                   <  1:5198378/150‑1 (MQ=255)

                                                                     cAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTtcag                                                                                  >  1:6520424/1‑150 (MQ=255)

                                                                     cAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGACATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTtcag                                                                                  <  2:6320409/150‑1 (MQ=255)

                                                                           gCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAg                                                                                                                            <  2:2844897/102‑1 (MQ=255)

                                                                            cccAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCa                                                                           <  2:460966/150‑1 (MQ=255)

                                                                            cccAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCa                                                                           <  1:588933/150‑1 (MQ=255)

                                                                              cagcagGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGcctta                                                                         >  1:4313094/1‑147 (MQ=255)

                                                                              cagcagGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGt                                                                         >  1:4359972/1‑150 (MQ=255)

                                                                                gcagGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGtgc                                                                       >  1:2783068/1‑150 (MQ=255)

                                                                                   ggTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGcgctgc                                                                    <  2:886609/150‑1 (MQ=255)

                                                                                    gTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGttt                                                                                     >  2:2968302/1‑132 (MQ=255)

                                                                                      gCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGctgcctg                                                                 >  2:1611103/1‑150 (MQ=255)

                                                                                       cccAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGctgcctgc                                                                >  2:5633376/1‑150 (MQ=255)

                                                                                         cAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCctgctg                                                              >  1:3446718/1‑150 (MQ=255)

                                                                                         cAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTCTCAGTCAGCCAGTGCTGCctgctg                                                              >  2:4577529/1‑150 (MQ=255)

                                                                                          aaGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCtgctgc                                                             <  1:6636960/150‑1 (MQ=255)

                                                                                                  cTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGaa                                                     <  1:1484987/150‑1 (MQ=255)

                                                                                                      cAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAAtc                                                 >  1:927475/1‑150 (MQ=255)

                                                                                                      cAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAAtc                                                 >  2:755488/1‑150 (MQ=255)

                                                                                                        ggCCCAACAGCGCGGAGAGGTGGCAGGACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAAtctc                                               <  1:6723554/150‑1 (MQ=255)

                                                                                                        ggCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAAtctc                                               <  1:6135863/150‑1 (MQ=255)

                                                                                                         gCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCg                                              <  1:2553724/150‑1 (MQ=255)

                                                                                                          cccAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGt                                             <  2:2440104/150‑1 (MQ=255)

                                                                                                          cccAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGt                                             >  2:206168/1‑150 (MQ=255)

                                                                                                             aaCAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCa                                          <  2:5452802/150‑1 (MQ=255)

                                                                                                                  cgcgGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGccc                                     <  1:6772164/150‑1 (MQ=255)

                                                                                                                     ggAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTc                                  <  2:1069465/150‑1 (MQ=255)

                                                                                                                          ggTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTgcgc                             >  2:6577591/1‑150 (MQ=255)

                                                                                                                           gTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGcc                            >  2:6185170/1‑150 (MQ=255)

                                                                                                                              gCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCaaa                         >  2:6376480/1‑150 (MQ=255)

                                                                                                                                  aaCATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCAAAGCGa                     <  2:4383094/150‑1 (MQ=255)

                                                                                                                                  aaCATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCAAAGCGa                     >  2:4354055/1‑150 (MQ=255)

                                                                                                                                  aaCATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCAAAGCGa                     <  2:13308/150‑1 (MQ=255)

                                                                                                                                     cccGCCCAATGCCGAAATCCGGCAATGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCAAAGCGAGCa                  >  2:3672738/3‑150 (MQ=255)

                                                                                                                                        gCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCAAAGCGAGCATgg               >  1:68704/1‑150 (MQ=255)

                                                                                                                                          ccATAGCCGGAATTCGCCATTGGTCGGGGTGATCTACAGGTTCGAATCGATCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCAGCCTGCTGCAGCGCGAAAATCTCGTTCAAGCCCTTCTGCGACAAAGCGAGCATGGCg             >  2:1846883/1‑150 (MQ=255)

                                                                                                                                           cATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCAAAGCGAGCATGGCGt            >  2:641041/1‑150 (MQ=255)

                                                                                                                                              tGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCAAAGCGAGCATGGCGTTGa         <  2:2735077/150‑1 (MQ=255)

                                                                                                                                               gCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCAAAGCGAGCATGGCGTTGaa        <  2:5097154/150‑1 (MQ=255)

                                                                                                                                                ccGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCAAAGCGAGCATGGCGTTGaaa       <  1:4778482/150‑1 (MQ=255)

                                                                                                                                                ccGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCAAAGCGAGCATGGCGTTGaaa       >  2:6440641/1‑150 (MQ=255)

                                                                                                                                                ccGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCAAAGCGAGCAGGGCGTTGaaa       <  1:1962891/150‑1 (MQ=255)

                                                                                                                                                 cGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCAAAGCGAGCATGGCGTTGAAAt      >  1:5430573/1‑150 (MQ=255)

                                                                                                                                                  ggATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCAAAGCGAGCATGGCGTTGAAATc     <  1:5850172/150‑1 (MQ=255)

                                                                                                                                                   gATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCAAAGCGAGCATGGCGTTGAAATCt    <  1:5950405/150‑1 (MQ=255)

                                                                                                                                                     tttCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCAAAGCGAGCATGGCGTTGAAATCTTc  <  2:6686072/150‑1 (MQ=255)

                                                                                                                                                     tttCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCAAAGCGAGCATGGCGTTGAAATCTTc  <  1:791340/150‑1 (MQ=255)

AAAGTTCAGCGCCTCTTTGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCAAAGCGAGCATGGCGTTGAAATCTTC  >  NC_002947/6045991‑6046289

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