BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,143,579	T→G	6.1%	S29A (TCG→GCG) ‡	PP_3643 →	oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,143,579	0	T	G	6.1%	758.6 / 26.5	224	S29A (TCG→GCG) ‡	PP_3643	oxidoreductase
Reads supporting (aligned to +/- strand): ref base T (100/109); new base G (4/10); total (105/119)
Fisher's exact test for biased strand distribution p-value = 1.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 227 total aligned reads displayed.

TATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCACCTGCTGTCCGGTGAACAGCAG  >  NC_002947/4143431‑4143728

tataCCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTc                                                                                                                                                      >  1:1691441/1‑150 (MQ=255)

 ataCCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCg                                                                                                                                                     <  2:2583274/150‑1 (MQ=255)

  taCCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCgg                                                                                                                                                    >  2:1817834/1‑150 (MQ=255)

    ccATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGgtc                                                                                                                                                  <  1:2995820/150‑1 (MQ=255)

     cATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGgtcg                                                                                                                                                 >  2:4366407/1‑150 (MQ=255)

     cATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGgtcg                                                                                                                                                 >  1:5103809/1‑150 (MQ=255)

      aTCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGgtcgt                                                                                                                                                >  2:3855514/1‑150 (MQ=255)

       tCAGACACATCAACCATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGgtcgtc                                                                                                                                               <  2:2247102/150‑1 (MQ=255)

        cAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGCCGACGACCGAGCCGtcgtcg                                                                                                                                              <  1:488877/150‑1 (MQ=255)

         aGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCAGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCgg                                                                                                                                             >  1:5275007/1‑150 (MQ=255)

          gACACATCAACGATAACAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACCGAGCCGtcgtcgtc                                                                                                                                            <  2:5568847/150‑3 (MQ=255)

          gACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGCCGACGACCGAGCCGtcgtcgtc                                                                                                                                            <  1:4786756/150‑3 (MQ=255)

          gACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGc                                                                                                                                            <  2:3205141/150‑1 (MQ=255)

          gACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGc                                                                                                                                            >  1:2736177/1‑150 (MQ=255)

          gACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGc                                                                                                                                            <  1:3112629/150‑1 (MQ=255)

           acacATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGCCGACGACCGAGCCGTCGTCGTcc                                                                                                                                           <  2:5427490/150‑1 (MQ=255)

           acacATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGCCGACGACCGAGCCGTCGTCGTcc                                                                                                                                           <  2:4824861/150‑1 (MQ=255)

           acacATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGcc                                                                                                                                           <  1:853322/150‑1 (MQ=255)

            cacaTCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGCCGACGACCGAGCCGTCGTCGTCCt                                                                                                                                          <  1:4508601/150‑1 (MQ=255)

            cacaTCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGCCGACGACCGAGCCGTCGTCGTCCt                                                                                                                                          <  2:575425/150‑1 (MQ=255)

              caTCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGa                                                                                                                                        <  1:3672531/150‑1 (MQ=255)

               aTCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGAc                                                                                                                                       >  2:2938432/1‑150 (MQ=255)

               aTCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGAc                                                                                                                                       >  1:2407487/1‑150 (MQ=255)

                tCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGAcc                                                                                                                                      >  1:723408/1‑150 (MQ=255)

                   aCGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGCCGACGACCGAGCCGTCGTCGTCCTgggtggt                                                                                                                                   <  1:3884397/150‑7 (MQ=38)

                      aTAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAgc                                                                                                                                <  1:167228/150‑1 (MQ=255)

                            aCCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTtc                                                                                                                          >  2:2716746/1‑150 (MQ=255)

                            aCCGCGGACCACCGCAGCCCTTCAGAGGTACACCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCCGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTtc                                                                                                                          <  1:945965/150‑1 (MQ=255)

                               gcgGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCAGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTctcc                                                                                                                       <  2:5275007/150‑1 (MQ=255)

                                 ggACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTctcctc                                                                                                                     <  1:5611568/150‑1 (MQ=255)

                                 gccccaccGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCCGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTctcctc                                                                                                                     <  1:2345282/147‑1 (MQ=255)

                                  gACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCCGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCtcctcc                                                                                                                    <  2:1111791/150‑1 (MQ=255)

                                          cAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGAt                                                                                                            >  2:4748004/1‑150 (MQ=255)

                                          cAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGAt                                                                                                            <  2:2639316/150‑1 (MQ=255)

                                          cAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGAt                                                                                                            <  2:625763/150‑1 (MQ=255)

                                           aGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGTGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATg                                                                                                           >  2:1311814/1‑150 (MQ=255)

                                             cccTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGa                                                                                                         <  2:3066234/150‑1 (MQ=255)

                                             cccTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGa                                                                                                         >  1:4197296/1‑150 (MQ=255)

                                              ccTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGAc                                                                                                        <  2:4911710/150‑1 (MQ=255)

                                               cTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGAcc                                                                                                       >  1:2704947/1‑150 (MQ=255)

                                               cTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGAcc                                                                                                       >  1:3251856/1‑150 (MQ=255)

                                                ttCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGAccc                                                                                                      >  2:4414574/1‑150 (MQ=255)

                                                gtCAGAGGTACTCCCGATGATCGACGCAACTGTCTACACGACCGTTCGGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGCCGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGAccc                                                                                                      <  1:1042905/149‑1 (MQ=255)

                                                   agagGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACccgcc                                                                                                   >  2:4569615/1‑150 (MQ=255)

                                                    gagGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACccgccg                                                                                                  <  1:3634616/150‑1 (MQ=255)

                                                       gTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTg                                                                                               <  2:3386779/150‑1 (MQ=255)

                                                        tACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTgg                                                                                              >  1:4110696/1‑150 (MQ=255)

                                                          cTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTg                                                                                            <  1:920459/150‑1 (MQ=255)

                                                           tCCCGATGATCGACGCAACTGTCTACAAGAACGATCTGGGCTCCTTCCCGTCCGGCGTTACCATCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGc                                                                                           >  2:3377270/1‑150 (MQ=255)

                                                               gatgatCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGtt                                                                                       >  2:3520424/1‑150 (MQ=255)

                                                                 tgatCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCt                                                                                     <  1:4665359/150‑1 (MQ=255)

                                                                   ctCGACGCCACGGTCTTCTAGCACGTTCTGGGCTCCTTCCCGTCCGGCGTGACCGTCATCACCACCGAGGACGAGGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCGGTCGATGGCCCCGCCGCTGGTGCTGTTCTGc                                                                                   <  2:1163253/149‑1 (MQ=255)

                                                                     cGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGccc                                                                                 >  2:1039176/1‑150 (MQ=255)

                                                                      gACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGcccc                                                                                <  1:1828795/150‑1 (MQ=255)

                                                                       aCGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCa                                                                               >  2:1539380/1‑150 (MQ=255)

                                                                       aCGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCa                                                                               >  1:896487/1‑150 (MQ=255)

                                                                        cGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCaa                                                                              <  2:1790179/150‑1 (MQ=255)

                                                                         gCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAAc                                                                             <  1:3335572/150‑1 (MQ=255)

                                                                           aaCTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTa                                                                           >  2:2024472/1‑150 (MQ=255)

                                                                            aCTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTac                                                                          <  1:600103/150‑1 (MQ=255)

                                                                              tGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCGTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGTCCCGCCGCTGGTGCTGTTCTGCCCCTACTacac                                                                        >  2:3306006/1‑150 (MQ=255)

                                                                                  tACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTcc                                                                    >  2:172471/1‑150 (MQ=255)

                                                                                  tACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTcc                                                                    <  2:2736201/150‑1 (MQ=255)

                                                                                     aaGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGAc                                                                 >  2:548952/1‑150 (MQ=255)

                                                                                     aaGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGAc                                                                 >  2:1518986/1‑150 (MQ=255)

                                                                                        aaCGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTcc                                                              >  2:3723536/1‑150 (MQ=255)

                                                                                           gTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTAc                                                           >  2:1931044/1‑150 (MQ=255)

                                                                                            ttCGGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTAcc                                                          <  1:3033497/150‑1 (MQ=255)

                                                                                             ctcaggCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCGCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTcccc                                                         >  2:402642/5‑150 (MQ=255)

                                                                                                 ggCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTa                                                     <  2:821521/150‑1 (MQ=255)

                                                                                                     ccTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGa                                                 <  2:2688049/150‑1 (MQ=255)

                                                                                                       ttCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATc                                               <  2:4871216/150‑1 (MQ=255)

                                                                                                       ttCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATc                                               >  2:3489960/1‑150 (MQ=255)

                                                                                                         cccGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCaa                                             <  2:4308789/150‑1 (MQ=255)

                                                                                                         cccGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCaa                                             <  2:785456/150‑1 (MQ=255)

                                                                                                           cGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGc                                           >  1:21190/1‑150 (MQ=255)

                                                                                                               cGGCGTTACCGTCATCACCACCCACCACGACGACTGCTCTGTCTTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGACGCTGTTGATCTTCTGCCCCAACTACACCTCAGACTCCTACCCGCTACTGAACAACCACaa                                       <  2:1688767/150‑1 (MQ=255)

                                                                                                                      aCCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCttt                                <  2:4081671/150‑1 (MQ=255)

                                                                                                                      aCCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCttt                                <  1:1313266/150‑1 (MQ=255)

                                                                                                                      aCCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCttt                                <  2:4299091/150‑1 (MQ=255)

                                                                                                                       ccatcatcaCCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTg                               <  1:3377270/147‑1 (MQ=255)

                                                                                                                              caccacCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCa                        <  1:5306267/150‑1 (MQ=255)

                                                                                                                               accaccCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCAc                       <  1:5366987/150‑1 (MQ=255)

                                                                                                                                ccaccCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGa                                                                  >  2:4266614/1‑106 (MQ=255)

                                                                                                                                ccaccCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGa                                                                  <  1:4266614/106‑1 (MQ=255)

                                                                                                                                 caccCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCACct                     >  2:4070726/1‑150 (MQ=255)

                                                                                                                                 caccCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCACct                     <  1:1215746/150‑1 (MQ=255)

                                                                                                                                   cccAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCACctgc                   >  1:5234573/1‑150 (MQ=255)

                                                                                                                                    ccAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCACctgct                  <  2:4955962/150‑1 (MQ=255)

                                                                                                                                     cAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCACctgctg                 <  2:2762192/150‑1 (MQ=255)

                                                                                                                                       ggACGACGTCGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCACCTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCACCTGCTGTc               <  2:1741261/150‑1 (MQ=255)

                                                                                                                                        gccgacgaccGAGCCGTCGTCGTCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCACCTGCTGTcc              >  2:2381552/3‑150 (MQ=37)

                                                                                                                                           gacgacgGCTCGGTCGTCGGCCTGACCGCCGGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCACCTGCTGTCCGGt           >  1:2434486/1‑150 (MQ=255)

                                                                                                                                           gacgacgGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCACCTGCTGTCCGGt           <  2:3324268/150‑1 (MQ=255)

                                                                                                                                           gacgacgGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCACCTGCTGTCCGGt           >  1:2511034/1‑150 (MQ=255)

                                                                                                                                           gacgacgGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCACCTGCTGTCCGGt           >  2:3254746/1‑150 (MQ=255)

                                                                                                                                             cgacgGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCACCTGCTGTCCGGTGa         <  1:2750102/150‑1 (MQ=255)

                                                                                                                                               acgGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCACCTGCTGTCCGGTGAAc       <  1:1039123/150‑1 (MQ=255)

                                                                                                                                                   cTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCACCTGCTGTCCGGTGAAcagca   >  1:3479135/1‑150 (MQ=255)

                                                                                                                                                   cTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCACCTGCTGTCCGGTGAAcagca   >  1:3707719/1‑150 (MQ=255)

TATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCACCTGCTGTCCGGTGAACAGCAG  >  NC_002947/4143431‑4143728

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 11 ≤` `ATCG/ATCG` `< 25 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`