BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,190,316	T→C	6.5%	R82R (CGT→CGC) ‡	PP_3684 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,190,316	0	T	C	6.5%	824.9 / 26.8	263	R82R (CGT→CGC) ‡	PP_3684	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (117/127); new base C (8/9); total (126/137)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 269 total aligned reads displayed.

GTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTCCGTATAGACA  >  NC_002947/4190167‑4190465

gTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGt                                                                                                                                                       >  1:1148320/1‑150 (MQ=255)

gTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGt                                                                                                                                                       <  2:6664462/150‑1 (MQ=255)

gTCTTCCCAGACAGAAGACCAGAAGAAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGt                                                                                                                                                       <  2:8423303/150‑1 (MQ=255)

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

      ccAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGAACGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGct                                                                                                                                                 >  2:5761356/1‑150 (MQ=255)

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

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

       cAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGctc                                                                                                                                                <  2:7800830/150‑1 (MQ=255)

         gACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCtctc                                                                                                                                              <  1:4056879/150‑1 (MQ=255)

           ccGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCTCACGGCTCTCaa                                                                                                                                            <  1:5901528/150‑1 (MQ=255)

           ccGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGACGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCaa                                                                                                                                            <  2:3150061/150‑1 (MQ=255)

            cGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCCTCGAGGCGGTGAACAGCCTCAACACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAg                                                                                                                                           <  2:7407844/150‑1 (MQ=255)

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

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

             gAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCTCACGGCTCTCAAgg                                                                                                                                          <  1:8038428/150‑2 (MQ=255)

              aaGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCa                                                                                                                                         >  2:161127/1‑150 (MQ=255)

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

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

                       aGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGa                                                                                                                                <  2:8573179/150‑1 (MQ=255)

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

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

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

                         cAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGAcc                                                                                                                              <  2:715671/150‑1 (MQ=255)

                               tAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGgcgc                                                                                                                        <  1:8851994/150‑1 (MQ=255)

                               tAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGgcgc                                                                                                                        <  1:8851560/150‑1 (MQ=255)

                                  ggCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATc                                                                                                                     >  2:6804744/1‑150 (MQ=255)

                                  ggCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATc                                                                                                                     >  1:7499071/1‑150 (MQ=255)

                                  ggCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATc                                                                                                                     >  1:6815923/1‑150 (MQ=255)

                                   gCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCg                                                                                                                    >  2:3738389/1‑150 (MQ=255)

                                       gAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGaaa                                                                                                                <  2:6077377/150‑1 (MQ=255)

                                       gAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGaaa                                                                                                                <  1:163608/150‑1 (MQ=255)

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

                                          cgcCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCg                                                                                                             >  2:7601441/1‑150 (MQ=255)

                                            ccAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGaa                                                                                                           <  1:2201147/150‑1 (MQ=255)

                                            ccAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGaa                                                                                                           <  1:1910176/150‑1 (MQ=255)

                                                ggttttgAGGCTGTTCACCGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACaaa                                                                                                       >  2:4855074/7‑150 (MQ=255)

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

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

                                                   tttgAGGCTGTTCACCGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAAcc                                                                                                    >  2:5489581/4‑150 (MQ=255)

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

                                                      gAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACGTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCAGa                                                                                                 >  1:4821592/1‑150 (MQ=255)

                                                      gAGGCGGTGAACAGCCTCAAAACCTGGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGa                                                                                                 <  2:5625150/150‑1 (MQ=255)

                                                       aGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAg                                                                                                <  1:3481096/150‑1 (MQ=255)

                                                            gTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCgctg                                                                                           >  2:1085329/1‑150 (MQ=255)

                                                             tGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCgctgc                                                                                          <  1:6437005/150‑1 (MQ=255)

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

                                                                  aGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTtgtg                                                                                     >  2:1384837/1‑150 (MQ=255)

                                                                   gCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTgg                                                                                    <  1:2376909/150‑1 (MQ=255)

                                                                   gCCTCAAAACCTGGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTgg                                                                                    <  1:6994877/150‑1 (MQ=255)

                                                                    ccTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGt                                                                                   >  1:2643/1‑150 (MQ=255)

                                                                    ccTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGt                                                                                   >  1:2663/1‑150 (MQ=255)

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

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

                                                                      tCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGgc                                                                                                                          <  1:2859903/109‑1 (MQ=255)

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

                                                                        aaaaCCTTGAGGGTCACTGACCGAGGCGGGTTGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCa                                                                               <  1:4528508/150‑1 (MQ=255)

                                                                          aaCCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGa                                                                             <  1:2479993/150‑1 (MQ=255)

                                                                              ttGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGAtctc                                                                         >  1:4333083/1‑150 (MQ=255)

                                                                              ttGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGAtctc                                                                         >  1:5244777/1‑150 (MQ=255)

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

                                                                                     tCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCct                                                                  <  2:8743761/150‑1 (MQ=255)

                                                                                            ccGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACt                                                           >  2:4706010/1‑150 (MQ=255)

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

                                                                                                 gCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGc                                                      >  2:4120819/1‑150 (MQ=255)

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

                                                                                                              cgaCCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCg                                         <  2:780570/150‑1 (MQ=255)

                                                                                                                aCCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTa                                       <  2:1148320/150‑1 (MQ=255)

                                                                                                                    tGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGc                                          <  2:392046/143‑1 (MQ=255)

                                                                                                                    tGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGc                                          >  1:392046/1‑143 (MQ=255)

                                                                                                                    tGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTAtttt                                   >  2:870680/1‑150 (MQ=255)

                                                                                                                    tGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTAtttt                                   <  1:8158790/150‑1 (MQ=255)

                                                                                                                     gAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGTAGCGACTTCAACTGCGTATTTtg                                  <  2:8815292/150‑1 (MQ=255)

                                                                                                                      aaGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTtgt                                 <  2:7039919/150‑1 (MQ=255)

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

                                                                                                                        gagaTACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCa                                                                               >  1:8162522/1‑102 (MQ=255)

                                                                                                                        gagaTACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCa                                                                               <  2:8162522/102‑1 (MQ=255)

                                                                                                                         agaTACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTACTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGcc                              <  2:5765511/150‑1 (MQ=255)

                                                                                                                          gaTACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGccc                             >  1:3571879/1‑150 (MQ=255)

                                                                                                                           attcGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCTAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCt                            <  2:4417255/147‑1 (MQ=255)

                                                                                                                            tACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCtt                           >  1:7078759/1‑150 (MQ=255)

                                                                                                                             ccGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCttt                          <  2:8229314/149‑1 (MQ=255)

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

                                                                                                                                   tgtgatCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTc                    >  2:4328058/1‑150 (MQ=255)

                                                                                                                                     tgatCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGa                  >  2:9007265/1‑150 (MQ=255)

                                                                                                                                        tgcTTGAGAGCCTCACGGCTCTCAAGGACTTGGTCGACCCGGCGCACCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACaa               >  2:4501863/3‑150 (MQ=255)

                                                                                                                                        tCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACaa               >  1:7171515/1‑150 (MQ=255)

                                                                                                                                        tCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACaa               >  2:2153519/1‑150 (MQ=255)

                                                                                                                                           ttGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTc            >  1:1714692/1‑150 (MQ=255)

                                                                                                                                           ttGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTc            >  1:5026429/1‑150 (MQ=255)

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

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

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

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

                                                                                                                                             gagaGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTCCg          <  1:3002883/150‑1 (MQ=255)

                                                                                                                                              agagCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTCCGt         >  1:4703718/1‑150 (MQ=255)

                                                                                                                                               gagCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTCCGta        <  1:5489581/150‑1 (MQ=255)

                                                                                                                                                 gCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTCCGtata      <  1:5988200/150‑1 (MQ=255)

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

                                                                                                                                                   cGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTCCGTATAGa    <  2:3543556/150‑1 (MQ=255)

                                                                                                                                                    gTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTCCGTATAGAc   >  1:9085425/1‑150 (MQ=255)

GTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTCCGTATAGACA  >  NC_002947/4190167‑4190465

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