BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,922,959	G→A	5.5%	V122I (GTC→ATC) ‡	PP_4328 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,922,959	0	G	A	5.5%	941.1 / 27.9	293	V122I (GTC→ATC) ‡	PP_4328	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (128/148); new base A (11/5); total (139/154)
Fisher's exact test for biased strand distribution p-value = 1.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 301 total aligned reads displayed.

CTTCAGGCCAACGCCAGCCAGCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTC  >  NC_002947/4922812‑4923108

 ttCAGGCCAACGCCAGCCAGCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCg                                                                                                                                                    >  2:787325/1‑150 (MQ=255)

 ttCAGGCCAACGCCAGCCAGCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCg                                                                                                                                                    <  2:1841880/150‑1 (MQ=255)

   cAGGCCAACGCCAGCAAGCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGACGGc                                                                                                                                                  >  1:531387/1‑150 (MQ=255)

    aGGCCAACGCCAGCCAGCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCa                                                                                                                                                 <  1:4212159/150‑1 (MQ=255)

    aGGCCAACGCCAGCCAGCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCa                                                                                                                                                 <  2:2341727/150‑1 (MQ=255)

       ccAACGCCAGCCAGCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCt                                                                                                                                              >  1:2034938/1‑150 (MQ=255)

       ccAACGCCAGCCAGCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCt                                                                                                                                              >  1:614816/1‑150 (MQ=255)

        cAACGCCAGCCAGCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCtt                                                                                                                                             <  1:248562/150‑1 (MQ=255)

        cAACGCCAGCCAGCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGACCGGCACCtt                                                                                                                                             <  2:5359983/150‑1 (MQ=255)

         aaCGCCAGCCAGCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTtg                                                                                                                                            <  2:3491471/150‑1 (MQ=255)

          aCGCCAGCCAGCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTtgc                                                                                                                                           >  2:5089683/1‑150 (MQ=255)

            gccagccaGCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGACCGGCACCTtgctg                                                                                                                                         <  2:3430850/150‑1 (MQ=255)

             ccagccagCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTtgctgc                                                                                                                                        <  1:6843352/150‑1 (MQ=255)

              cagccagcCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTgca                                                                                                                                       <  1:6342861/150‑1 (MQ=255)

                     ccGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTgcaaggcaag                                                                                                                                <  1:7337377/150‑1 (MQ=255)

                     ccGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCACCAAGGTGCCGGTCGGCACCTTGCTgcaaggcaag                                                                                                                                <  2:6807570/150‑1 (MQ=255)

                       gCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGt                                                                                                                              >  2:5532178/1‑150 (MQ=255)

                         tGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTtgctgc                                                                                                                                        <  2:5898862/138‑1 (MQ=255)

                          gCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCt                                                                                                                           <  2:4639568/150‑1 (MQ=255)

                            cccAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGa                                                                                                                         >  1:1214080/1‑150 (MQ=255)

                            cccAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGa                                                                                                                         <  2:3861161/150‑1 (MQ=255)

                                 ggCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGaccaac                                                                                                                    <  1:1340958/150‑1 (MQ=255)

                                    aGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGa                                                                                                                         <  2:4699252/142‑1 (MQ=255)

                                     gCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCACGGTCCTGACCAACCAgc                                                                                                                >  2:4552137/1‑150 (MQ=255)

                                       cAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAgctg                                                                                                              <  1:3908427/150‑1 (MQ=255)

                                         gATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAgctgct                                                                                                            >  1:3026468/1‑150 (MQ=255)

                                          aTCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGctgctg                                                                                                           <  1:4684506/150‑1 (MQ=255)

                                                 tCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAgg                                                                                                    >  2:1528183/1‑150 (MQ=255)

                                                   aGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCa                                                                                                  <  1:6527102/150‑1 (MQ=255)

                                                     ccAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCAcc                                                                                                >  1:1847240/1‑150 (MQ=255)

                                                      cAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCAccc                                                                                               <  1:579659/150‑1 (MQ=255)

                                                      cAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCAccc                                                                                               >  2:7325686/1‑150 (MQ=255)

                                                       aGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCg                                                                                              <  2:2034938/150‑1 (MQ=255)

                                                           cGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCa                                                                                          >  2:6339943/1‑150 (MQ=255)

                                                              aaGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAAcc                                                                                       >  2:6591353/1‑150 (MQ=255)

                                                                  aaTCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGcc                                                                                   <  1:2706049/150‑1 (MQ=255)

                                                                   aTCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCgg                                                                                     >  1:1030230/1‑147 (MQ=255)

                                                                   aTCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCg                                                                                  <  2:7236391/150‑1 (MQ=255)

                                                                    tCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCgt                                                                                 >  1:606100/1‑150 (MQ=255)

                                                                        cTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTAc                                                                             >  2:4833903/1‑150 (MQ=255)

                                                                         tGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTAcc                                                                            <  2:2418053/150‑1 (MQ=255)

                                                                          ggCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCg                                                                           >  1:359705/1‑150 (MQ=255)

                                                                             gATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCTAGGTGCCGGTCGGCACCTTGCTGGAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCATCCGGCCGTGTACCGGgc                                                                        <  1:6753644/150‑1 (MQ=255)

                                                                             gATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGgc                                                                        >  2:1045297/1‑150 (MQ=255)

                                                                               tCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGAACCCGGCCAACCGGCCGTGTACCGGgcgc                                                                      >  2:5219541/1‑150 (MQ=255)

                                                                                cATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCt                                                                     <  2:1355296/150‑1 (MQ=255)

                                                                                cATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCt                                                                     >  2:1286752/1‑150 (MQ=255)

                                                                                   gctgcAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGt                                                                  <  1:2159827/150‑1 (MQ=255)

                                                                                    ctgcAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTc                                                                 >  1:276600/1‑150 (MQ=255)

                                                                                      gcAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCAccc                                                                                               >  1:5194118/1‑118 (MQ=255)

                                                                                      gcAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCAccc                                                                                               <  1:4069750/118‑1 (MQ=255)

                                                                                       cAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGc                                                              <  1:4545539/150‑1 (MQ=255)

                                                                                        aaCAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGcc                                                             >  2:3686363/1‑150 (MQ=255)

                                                                                        aaCAGGCCAGCGCCAGCCAGGTCGCCACACTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCGGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGcc                                                             >  2:5998230/1‑150 (MQ=255)

                                                                                           aGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGc                                                          <  1:6637091/150‑1 (MQ=255)

                                                                                            ggCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCt                                                         <  2:3551293/150‑1 (MQ=255)

                                                                                               cagcgccagcCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCaa                                                      >  1:558342/1‑150 (MQ=255)

                                                                                                agcgccagcCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAAc                                                     >  2:7194626/1‑150 (MQ=255)

                                                                                                 gcgccagcCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACa                                                    <  1:4941108/150‑1 (MQ=255)

                                                                                                   gccagccaGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGc                                                  >  2:6067623/1‑150 (MQ=255)

                                                                                                   gccagccaGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGc                                                  >  1:1546439/1‑150 (MQ=255)

                                                                                                    ccagccagGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCa                                                 <  1:3220843/150‑1 (MQ=255)

                                                                                                    ccagccagGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCa                                                 >  1:976897/1‑150 (MQ=255)

                                                                                                     cagccagGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCaa                                                >  1:764237/1‑150 (MQ=255)

                                                                                                      agccagGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGCTCCTGACCAACCAGCTGCAGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAAc                                               <  2:4308537/150‑1 (MQ=255)

                                                                                                       gccagGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTCCCGGGCGCTGGTCAGCCTGCTCAACAGCAAcc                                              >  2:6275167/1‑150 (MQ=255)

                                                                                                           ggTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGc                                          <  1:2663908/150‑1 (MQ=255)

                                                                                                                      tCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGc                                                              <  2:3148484/119‑1 (MQ=255)

                                                                                                                      tCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCacac                               <  2:64405/150‑1 (MQ=255)

                                                                                                                      tCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCacac                               >  1:3789550/1‑150 (MQ=255)

                                                                                                                       cACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACt                              <  2:3200723/150‑1 (MQ=255)

                                                                                                                       cACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACt                              <  1:2561178/150‑1 (MQ=255)

                                                                                                                        aCCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTg                             <  1:5098612/150‑1 (MQ=255)

                                                                                                                         cccAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGa                            <  2:614816/150‑1 (MQ=255)

                                                                                                                         cccAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGa                            <  1:6615090/150‑1 (MQ=255)

                                                                                                                          ccAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAg                           >  1:2928205/1‑150 (MQ=255)

                                                                                                                              cTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATc                       <  2:2577733/150‑1 (MQ=255)

                                                                                                                               tGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCg                      >  2:6997406/1‑150 (MQ=255)

                                                                                                                                 gACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGAc                    >  1:710398/1‑150 (MQ=255)

                                                                                                                                  acacCAGCATGGTGCCGGTGGTCACCTTGCTGCAAGGCAAGTTCCTGACCAACCAGCTGCTGGCCAAGGCAGCCGGCCAACCGGCGGTGTTCAGGGCGATGGTCAGCCTGCTCATCATCAACCAGTCCTTTGCCAGTCTGAGCATCGACa                   <  2:2099050/150‑1 (MQ=255)

                                                                                                                                  acacCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACa                   <  2:5312149/150‑1 (MQ=255)

                                                                                                                                    acCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGa                            >  1:6956591/1‑139 (MQ=255)

                                                                                                                                    acCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGa                            <  2:6956591/139‑1 (MQ=255)

                                                                                                                                      cAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGc                                                              <  2:4081763/103‑1 (MQ=255)

                                                                                                                                      cAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCACAAAACCAACCAGGCCGGTGCCACACTGAGCATCGACAcccc               >  2:1178994/1‑150 (MQ=255)

                                                                                                                                       aGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCACCCGGCCGTGTACCGGGCGCGGGTCAGCCTGCTCAACCGCAACCAGGCCGGTGCCACACTGAGCATCGTCAgcccg              >  2:4183834/1‑150 (MQ=255)

                                                                                                                                       aGCAAGGTGCCGACCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAgcccg              >  1:1711503/1‑150 (MQ=255)

                                                                                                                                       aGCAAGGTGCCGACCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAgcccg              >  2:558163/1‑150 (MQ=255)

                                                                                                                                        gCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAgcccgc             <  2:7109550/150‑1 (MQ=255)

                                                                                                                                        gCAAGGTGCCGACCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAgcccgc             >  1:5471508/1‑150 (MQ=255)

                                                                                                                                        gCAAGGTGCCGACCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAgcccgc             >  1:5274756/1‑150 (MQ=255)

                                                                                                                                        gCAAGGTGCCGACCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAgcccgc             >  2:4434839/1‑150 (MQ=255)

                                                                                                                                         cAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCa                                                                                                      <  1:1093078/60‑1 (MQ=255)

                                                                                                                                         cAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCa                                                                                                      <  1:4829601/60‑1 (MQ=255)

                                                                                                                                         cAAGGTGCCGACCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAgcccgcg            >  2:3282280/1‑150 (MQ=255)

                                                                                                                                         cAAGGTGCCGACCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCGGGGCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAgcccgcg            >  1:4338137/1‑150 (MQ=255)

                                                                                                                                               gCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCt      <  2:6658945/150‑1 (MQ=255)

                                                                                                                                                ccGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTg     <  2:2632272/150‑1 (MQ=255)

                                                                                                                                                  ggTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTggg   >  2:3009539/1‑149 (MQ=255)

                                                                                                                                                  ggTCGGCACCTTGCTGCAAGGCAAGGTCCAGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGt   >  1:3757653/1‑150 (MQ=255)

CTTCAGGCCAACGCCAGCCAGCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTC  >  NC_002947/4922812‑4923108

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