BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,436,263	T→G	5.8%	L347W (TTG→TGG) ‡	PP_3925 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,436,263	0	T	G	5.8%	1198.5 / 42.2	353	L347W (TTG→TGG) ‡	PP_3925	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (194/137); new base G (9/12); total (204/149)
Fisher's exact test for biased strand distribution p-value = 1.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 356 total aligned reads displayed.

GACAACTTCGATTACTACTATTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAGGCCATCGCC  >  NC_002947/4436115‑4436411

  cAACTTCGATTACTACTATTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCa                                                                                                                                                   <  2:1649173/150‑1 (MQ=255)

   aaCTTCGATTACTACTATTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAg                                                                                                                                                  <  1:4486278/150‑1 (MQ=255)

     cTTCGATTACTACTATTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGt                                                                                                                                                >  2:882102/1‑150 (MQ=255)

      ttCGATTACTACTATTGCTATGGCCTGTTCCGCCTGGCGGGCATCGTTCAGCAGATCTCCTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTAACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTc                                                                                                                                               >  2:6540548/1‑150 (MQ=255)

      ttCGATTACTACTATTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTc                                                                                                                                               >  1:6308627/1‑150 (MQ=255)

          aTTACTACTATTGCTATGGCCTGTTCCGCCTGGCATGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTGCAAGGTCAttt                                                                                                                                           <  1:5238571/150‑1 (MQ=255)

          aTTACTACTATTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCTCTTGCTGGAACAAATGACCTGCAAGGTCAttt                                                                                                                                           <  2:6246942/150‑1 (MQ=255)

          aTTACTACTATTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTGCAAGGTCAttt                                                                                                                                           <  1:3507093/150‑1 (MQ=255)

              ctactaTTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCaa                                                                                                                                       <  2:534664/150‑1 (MQ=255)

              ctactaTTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCaa                                                                                                                                       <  2:4392072/150‑1 (MQ=255)

              ctactaTTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTGCAAGGTCAtttgttc                                                                                                                                       <  2:2874006/150‑5 (MQ=255)

               tactaTTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAg                                                                                                                                      >  1:819704/1‑150 (MQ=255)

               tactaTTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTAGCAGGTCATATCCAAg                                                                                                                                      >  2:3230566/1‑150 (MQ=255)

                actaTTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTGCAAGGTCAtttgttcca                                                                                                                                     <  1:1355803/150‑7 (MQ=255)

                  taTTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGAACTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATATCCAAGTcc                                                                                                                                   >  1:1101948/1‑150 (MQ=255)

                   aTTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCa                                                                                                                                  <  2:266917/150‑1 (MQ=255)

                    ttGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAg                                                                                                                                 >  2:2724372/1‑150 (MQ=255)

                     tGCTATGGCCTGTTCCGCCTGGCAGGCCTCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGc                                                                                                                                <  1:2662428/150‑1 (MQ=255)

                       cTATGTCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCct                                                                                                                              >  2:6434253/1‑150 (MQ=255)

                        tATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCctc                                                                                                                             >  1:5723891/1‑150 (MQ=255)

                         aTGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCctct                                                                                                                            <  1:4312813/150‑1 (MQ=255)

                         aTGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCctct                                                                                                                            <  2:3529411/150‑1 (MQ=255)

                          tGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTg                                                                                                                           >  2:4652973/1‑150 (MQ=255)

                             ccTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGc                                                                                                                        <  2:6308246/150‑1 (MQ=255)

                             ccTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGc                                                                                                                        >  1:406959/1‑150 (MQ=255)

                              cTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGca                                                                                                                       >  1:2268085/1‑150 (MQ=255)

                              cTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGca                                                                                                                       <  1:1980788/150‑1 (MQ=255)

                              cTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGca                                                                                                                       >  2:6045440/1‑150 (MQ=255)

                               tGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGcac                                                                                                                      >  1:5979269/1‑150 (MQ=255)

                               tGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGcac                                                                                                                      >  2:2438533/1‑150 (MQ=255)

                               tGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGcac                                                                                                                      >  2:3709850/1‑150 (MQ=255)

                                      ccTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCac                                                                                                               >  2:5882128/1‑150 (MQ=255)

                                        tGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCacga                                                                                                             >  2:5799065/1‑150 (MQ=255)

                                           caggcaTCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCacgacga                                                                                                          >  1:6567332/1‑150 (MQ=255)

                                            aggcaTCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCacgacgac                                                                                                         <  2:3318200/150‑1 (MQ=255)

                                            aggcaTCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCacgacgac                                                                                                         >  2:1021865/1‑150 (MQ=255)

                                                aTCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGAt                                                                                                     <  2:6717820/150‑1 (MQ=255)

                                                aTCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGAt                                                                                                     <  1:3386431/150‑1 (MQ=255)

                                                 tCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATa                                                                                                    <  1:2902119/150‑1 (MQ=255)

                                                  cGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATaa                                                                                                   <  1:6447611/150‑1 (MQ=255)

                                                   gTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAAc                                                                                                  >  2:5482751/1‑150 (MQ=255)

                                                    ttCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACa                                                                                                 >  2:3064000/1‑150 (MQ=255)

                                                     tCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACaa                                                                                                <  2:1386238/150‑1 (MQ=255)

                                                      cagcagATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAg                                                                                               <  2:2095578/150‑1 (MQ=255)

                                                        gcagATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGa                                                                                             >  2:3336037/1‑150 (MQ=255)

                                                        gcagATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGa                                                                                             <  1:1761571/150‑1 (MQ=255)

                                                         cagATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGaa                                                                                            <  2:2308165/150‑1 (MQ=255)

                                                         cagATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGaa                                                                                            <  1:5344231/150‑1 (MQ=255)

                                                          agATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGaaa                                                                                           >  2:5379210/1‑150 (MQ=255)

                                                           gATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCTCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGaaaa                                                                                          >  1:3247932/1‑150 (MQ=255)

                                                            aTCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAAc                                                                                         <  1:4619070/150‑1 (MQ=255)

                                                            aTCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAAc                                                                                         <  1:681951/150‑1 (MQ=255)

                                                              ctactaCCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAg                                                                                       <  2:5154818/150‑1 (MQ=255)

                                                                 ctacCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCAt                                                                                    >  1:2566357/1‑150 (MQ=255)

                                                                     cGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTcc                                                                                <  2:4334020/150‑1 (MQ=255)

                                                                      gCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCa                                                                               >  1:6656306/1‑150 (MQ=255)

                                                                           tacCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGAcc                                                                          >  2:6866968/1‑150 (MQ=255)

                                                                             ccATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCa                                                                        >  2:3597008/1‑150 (MQ=255)

                                                                             ccATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCa                                                                        >  2:6080484/1‑150 (MQ=255)

                                                                                   ccAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTcc                                                                                                                                         <  2:4174686/79‑1 (MQ=255)

                                                                                   ccAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCATGACCCACCTGtt                                                                  >  2:3411160/1‑150 (MQ=255)

                                                                                     aGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCg                                                                >  1:3147343/1‑150 (MQ=255)

                                                                                         ccAGGACAAGCGTTTCTCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCt                                                            >  2:468925/1‑150 (MQ=255)

                                                                                         ccAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCt                                                            >  1:1009040/1‑150 (MQ=255)

                                                                                         ccAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCt                                                            >  2:2522844/1‑150 (MQ=255)

                                                                                         ccAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCt                                                            >  2:1758898/1‑150 (MQ=255)

                                                                                         ccAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCGGTTCGACCt                                                            >  1:4910053/1‑150 (MQ=255)

                                                                                           aGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCg                                                          >  2:5512404/1‑150 (MQ=255)

                                                                                                aaGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGc                                                     >  2:2174603/1‑150 (MQ=255)

                                                                                                aaGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGc                                                     >  2:3917197/1‑150 (MQ=255)

                                                                                                     tttCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGAt                                                >  1:4689640/1‑150 (MQ=255)

                                                                                                        cGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGc                                             <  2:4991552/150‑1 (MQ=255)

                                                                                                         gCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGcc                                            >  2:6160465/1‑150 (MQ=255)

                                                                                                            cAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCttt                                         <  2:4549403/150‑1 (MQ=255)

                                                                                                                 cattcaCATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTc                                    <  1:5427941/150‑1 (MQ=255)

                                                                                                                  attcaCATGAACCGCTTGCTTGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTcc                                   <  1:4416517/150‑1 (MQ=255)

                                                                                                                   ttcaCATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCg                                  <  1:4810822/150‑1 (MQ=255)

                                                                                                                   ttcaCATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCg                                  <  1:1909713/150‑1 (MQ=255)

                                                                                                                    tcaCATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCgg                                 >  1:5823726/1‑150 (MQ=255)

                                                                                                                    tcaCATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCgg                                 >  2:524942/1‑150 (MQ=255)

                                                                                                                       caTGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGgcgc                              >  1:1512671/1‑150 (MQ=255)

                                                                                                                       caTGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGgcgc                              >  1:3438532/1‑150 (MQ=255)

                                                                                                                         tGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCa                            >  2:4290319/1‑150 (MQ=255)

                                                                                                                            aCCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGcc                         <  1:3978966/150‑1 (MQ=255)

                                                                                                                             ccGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCg                        >  1:5501456/1‑150 (MQ=255)

                                                                                                                              cGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGt                       <  1:2438533/150‑1 (MQ=255)

                                                                                                                              cGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGt                       <  2:1473034/150‑1 (MQ=255)

                                                                                                                                 aggCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGc                    >  1:5041199/3‑150 (MQ=255)

                                                                                                                                     tGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCg                >  1:2492732/1‑150 (MQ=255)

                                                                                                                                       ttggaaaTGACCTGCAAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGc              >  2:3968192/5‑150 (MQ=255)

                                                                                                                                       ttggaaaTGACCTGCAAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGc              >  1:2304418/5‑150 (MQ=255)

                                                                                                                                       gAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGc              >  1:6206879/1‑150 (MQ=255)

                                                                                                                                       gAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGc              >  2:1914569/1‑150 (MQ=255)

                                                                                                                                         aCAAATGAACTTGCAGGTCATTTCCTAGTCCAGCCTCTGAGCACAGCCGACGACGACGGATATCAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCTGCGa            <  2:1101948/150‑1 (MQ=255)

                                                                                                                                          gaaaTGACCTGCAAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAg           >  1:3072956/2‑150 (MQ=255)

                                                                                                                                           aaaTGACCTGCAAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAgg          >  1:4283062/1‑150 (MQ=255)

                                                                                                                                             aTGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAGGcc        >  1:1290775/1‑150 (MQ=255)

                                                                                                                                                aCCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAGGCCATc     <  1:5882128/150‑1 (MQ=255)

                                                                                                                                                aCCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAGGCCATc     <  1:3709850/150‑1 (MQ=255)

                                                                                                                                                   ttGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAGGCCATCGcc  >  1:3635954/1‑150 (MQ=255)

GACAACTTCGATTACTACTATTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAGGCCATCGCC  >  NC_002947/4436115‑4436411

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