BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,033,816	T→A	9.4%	F145I (TTC→ATC)	PP_0891 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,033,816	0	T	A	9.4%	480.0 / 32.1	160	F145I (TTC→ATC)	PP_0891	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (68/76); new base A (7/8); total (76/84)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01

Only 100 of 161 total aligned reads displayed.

AGTGAAGGCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCTGATCTTCACTGCGCG  >  NC_002947/1033668‑1033961

 gTGAAGGCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTc                                                                                                                                                  >  2:1093086/1‑150 (MQ=255)

  tGAAGGCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCACCTGGAAGCGCCGCGCATcc                                                                                                                                                 <  1:5074152/150‑1 (MQ=255)

   gAAGGCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCa                                                                                                                                                >  2:4283107/1‑150 (MQ=255)

   gAAGGCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGAAGCGCCGCGCATCCa                                                                                                                                                <  1:3777980/150‑1 (MQ=255)

      ggCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAgg‑t                                                                                                                                            <  2:1706990/150‑1 (MQ=255)

       gCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAgg‑tg                                                                                                                                           >  2:1790875/1‑150 (MQ=255)

       gCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAgg‑tg                                                                                                                                           >  1:4120194/1‑150 (MQ=255)

        cacaGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCGGACGGTCACCCTGCATGGCCAGTCGCACAGCGACACTGTCCTGGTGCTGGCCACCCGCCTGGAAGCGCCGCGCATCCAgg‑cgg                                                                                                                                          <  2:1022875/150‑1 (MQ=255)

         acaGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAgg‑tggt                                                                                                                                         >  1:724543/1‑150 (MQ=255)

           aGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGa                                                                                                                                       <  2:2847007/150‑1 (MQ=255)

                 aCCTCCAGCAGATCAATGACATGTCCTACCGTCAACATATCGAATTGCGCCTGCCGCTGACCTTCACACTGCATGCCCAGTCGCACACCAACAAACACCTCAAGAAGTCCACCCGACTGGAAGCGCGGCGCTTCCAGC‑TCGAGACCCTTg                                                                                                                                 >  2:4046420/1‑150 (MQ=255)

                   cTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCCGCGCTTCCAGG‑TGGTGACCCTTGAg                                                                                                                               >  2:844629/1‑150 (MQ=255)

                    tGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGc                                                                                                                              <  2:2633207/150‑1 (MQ=255)

                      cAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCct                                                                                                                            >  2:1499625/1‑150 (MQ=255)

                           gATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGt                                                                                                                       >  1:1400204/1‑150 (MQ=255)

                           gATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGt                                                                                                                       >  1:2552665/1‑150 (MQ=255)

                           gATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGt                                                                                                                       >  1:940738/1‑150 (MQ=255)

                             tCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGtgc                                                                                                                     >  1:1261350/1‑150 (MQ=255)

                                aTGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGtgctgc                                                                                                                  >  1:404377/1‑150 (MQ=255)

                                aTGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGtgctgc                                                                                                                  <  1:4608092/150‑1 (MQ=255)

                                        gCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGaa                                                                                                          <  2:2346922/150‑1 (MQ=255)

                                         ccGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAg                                                                                                         >  2:2674430/1‑150 (MQ=255)

                                         ccGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAg                                                                                                         >  2:2655120/1‑150 (MQ=255)

                                         ccGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCGACAATGTCCGGTTGCTGGCCACCCGCCGGGATGCGCGGCGCGTCCTGG‑TGGTGACCCTGGTGCCTCTGGTGCTGCGCGCCGAAg                                                                                                         >  2:2480712/1‑150 (MQ=255)

                                             cGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGAttt                                                                                                     >  2:4799008/1‑150 (MQ=255)

                                              ggTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGAtttt                                                                                                    <  1:5045038/150‑1 (MQ=255)

                                              ggTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGAtttt                                                                                                    <  1:1301076/150‑1 (MQ=255)

                                               gTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTg                                                                                                   >  1:2776311/1‑150 (MQ=255)

                                               gTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTg                                                                                                   <  1:1597008/150‑1 (MQ=255)

                                                tGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTgg                                                                                                  >  1:1044812/1‑150 (MQ=255)

                                                 gCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGc                                                                                                 >  1:3986602/1‑150 (MQ=255)

                                                   acaGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCtt                                                                                               >  2:725587/1‑150 (MQ=255)

                                                   acaGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCtt                                                                                               <  1:2416264/150‑1 (MQ=255)

                                                      gATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTtgct                                                                                            <  2:3515890/150‑1 (MQ=255)

                                                             ttGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACACTGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCTAAGATTTTGGCTTGCTGCCgggg                                                                                     <  2:2153294/150‑1 (MQ=255)

                                                                cgcCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTg                                                                                  <  1:4725702/150‑1 (MQ=255)

                                                                 gcCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATCTTGGCTTGCTGCCGGGGCTgg                                                                                 <  2:2877097/150‑1 (MQ=255)

                                                                    tGCCGCTGTCGGTCACCCTGCTTGGCCAGTCGCACAGCTACAATGTCGTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGTAGATTTTGGCTTGCTGCCGGGGCTGGaaa                                                                              <  2:4308886/150‑1 (MQ=255)

                                                                    tGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGaaa                                                                              <  1:2821970/150‑1 (MQ=255)

                                                                      ccGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGc                                                                            <  1:2783080/150‑1 (MQ=255)

                                                                        gCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCt                                                                          <  1:2368098/150‑1 (MQ=255)

                                                                              ggTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCaa                                                                    <  2:4567658/150‑1 (MQ=255)

                                                                              ggTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCaa                                                                    >  2:1284839/1‑150 (MQ=255)

                                                                               gTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAg                                                                   <  1:3829306/150‑1 (MQ=255)

                                                                                  aCCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTc                                                                >  1:93581/1‑150 (MQ=255)

                                                                                   cccTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCg                                                               >  2:4681508/1‑150 (MQ=255)

                                                                                    ccTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGc                                                              <  2:3599187/150‑1 (MQ=255)

                                                                                    ccTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGc                                                              >  2:3746194/1‑150 (MQ=255)

                                                                                    ccTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGc                                                              <  1:4745558/150‑1 (MQ=255)

                                                                                      tGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCa                                                            >  1:753132/1‑150 (MQ=255)

                                                                                      tGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCa                                                            >  1:753319/1‑150 (MQ=255)

                                                                                      tGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCa                                                            >  1:2737874/1‑150 (MQ=255)

                                                                                      tGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCACG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCa                                                            <  1:4613449/150‑1 (MQ=255)

                                                                                       gCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAg                                                           >  1:1458293/1‑150 (MQ=255)

                                                                                         aTGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGt                                                         <  2:745354/150‑1 (MQ=255)

                                                                                         aTGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGt                                                         >  1:1482958/1‑150 (MQ=255)

                                                                                            gCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGa                                                      >  1:416588/1‑150 (MQ=255)

                                                                                             ccAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGaa                                                     <  2:3978116/150‑1 (MQ=255)

                                                                                             ccAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGaa                                                     <  2:3978117/150‑1 (MQ=255)

                                                                                             ccAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGaa                                                     <  2:4059855/150‑1 (MQ=255)

                                                                                             ccAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGaa                                                     <  1:1779956/150‑1 (MQ=255)

                                                                                             ccAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGaa                                                     <  1:3833366/150‑1 (MQ=255)

                                                                                             ccAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGaa                                                     <  2:1196968/150‑1 (MQ=255)

                                                                                                gTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAgcc                                                  <  1:1143989/150‑3 (MQ=255)

                                                                                                gTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTc                                                  <  1:1144007/150‑1 (MQ=255)

                                                                                                    cacaGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATc                                              >  1:3730616/1‑150 (MQ=255)

                                                                                                     acaGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCa                                             >  1:554497/1‑150 (MQ=255)

                                                                                                      caGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCaa                                            <  1:2050820/150‑1 (MQ=255)

                                                                                                         cTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAAccc                                         <  1:1340727/150‑1 (MQ=255)

                                                                                                         cTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAAccc                                         <  1:1340832/150‑1 (MQ=255)

                                                                                                         cTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAAccc                                         >  1:1544308/1‑150 (MQ=255)

                                                                                                          tACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCGGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCg                                        >  2:394131/1‑150 (MQ=255)

                                                                                                          tACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCg                                        >  2:4363512/1‑150 (MQ=255)

                                                                                                          tACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCg                                        >  2:1784953/1‑150 (MQ=255)

                                                                                                            cAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTc                                      <  2:4889998/150‑1 (MQ=255)

                                                                                                            cAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTc                                      >  1:649425/1‑150 (MQ=255)

                                                                                                            cAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCAAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTc                                      <  1:2128620/150‑1 (MQ=255)

                                                                                                             aaTGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCg                                     >  1:4066078/1‑150 (MQ=255)

                                                                                                                   cTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCg                               <  2:4778120/150‑1 (MQ=255)

                                                                                                                    tGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCgg                              >  1:4214800/1‑150 (MQ=255)

                                                                                                                    tGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAAACCGTCGGTACCgg                              >  1:4798002/1‑150 (MQ=255)

                                                                                                                     gTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGt                             >  1:1814338/1‑150 (MQ=255)

                                                                                                                      ttGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGtt                            <  2:617566/150‑1 (MQ=255)

                                                                                                                           ggCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAgcgc                       <  2:1997714/150‑1 (MQ=255)

                                                                                                                             gtcaccaccTGGAAGCGCCGCGCATCCAGGCGGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGcgcg                      >  2:1424276/8‑150 (MQ=38)

                                                                                                                               aCCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTg                   <  2:3188992/150‑1 (MQ=255)

                                                                                                                               aCCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTg                   >  1:1622731/1‑150 (MQ=255)

                                                                                                                                 ccGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCt                 >  1:3205011/1‑150 (MQ=255)

                                                                                                                                 ccGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCt                 <  2:130003/150‑1 (MQ=255)

                                                                                                                                  caccTGGAAGCGCCGCGCATCCAGGCGGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCt                 >  2:3289982/3‑150 (MQ=255)

                                                                                                                                   gCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCTGa               <  2:554497/150‑1 (MQ=255)

                                                                                                                                      tGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCTGATCt            <  2:1442422/150‑1 (MQ=255)

                                                                                                                                      tGGAAGCGCCGCGCATCCAGGCGGGTGGCCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGACAGCCTGCGCACGTTCGCCCGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCTGATc             >  1:1089866/1‑150 (MQ=255)

                                                                                                                                      tGGAAGCGCCGCGCATCCAGGCGGGTGGCCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCTGATc             >  1:87656/1‑150 (MQ=255)

                                                                                                                                       ggATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCTGATCtt           >  2:1414143/1‑150 (MQ=255)

                                                                                                                                         aTGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCTGATCTTCa         <  2:220867/150‑1 (MQ=255)

                                                                                                                                            cgccgcGCATCCAGGCGGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCTGATCTTCACt       <  1:1424601/150‑1 (MQ=255)

                                                                                                                                             gcggcgCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCTGATCTTCACTgc     >  2:5086328/1‑150 (MQ=255)

                                                                                                                                             gcggcgCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCTGATCTTCACTgc     >  2:1499099/1‑150 (MQ=255)

                                                                                                                                                gcgcTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCTGATCTTCACTgcgcg  >  1:2295456/1‑150 (MQ=255)

AGTGAAGGCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGG‑TGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCTGATCTTCACTGCGCG  >  NC_002947/1033668‑1033961

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